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STABILITY AND PERFORMANCE OF MANNED CONTROL SYSTEMS 


by 
Alfred Skolnick 


Adviser: John G. Truxal 


Submitted in partial fulfillment of the requirements 


for the degree of Doctor of Philosophy (Systems) 


The human-operator transform is usually specified as a function of the 
controlled-element characteristics and type of input power spectrum. Thus, any 
single description for the human operator really is valid only for a set of limited 
conditions. This study demonstrates the human-operator transfer function ap- 
proach may be generalized; from a variety of published experimental data, ''capa- 
bility bounds" upon the transfer function parameters are formed. With such ranges 
defined, these parameters, forming atype of variable structure, can be used as a 
mapping function to display in the complex plane, the boundaries of human- adaptive 
Capacity. These boundaries contain a collection of ''critical points", similar to the 
-1+j0 point on the usual Nyquist diagram; when the plant Nyquist contour is plot- 


ted in conjunction with these bounded areas (at particular frequencies), interpreta- 
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tion of system stability characteristics is possible. Several plots of such bounda- 
ries are presented for important frequencies in the human range. A man-moni- 
tored system configuration also appears to be susceptible to such stzbility analy- 
sis, 

Experimental verification of the theory is obtained by analyzing certain 
systems which have been studied empirically in the past and for which results 
have been published in the literature. To this end, transparent overlays are pro- 
vided which facilitate stability determination. The theoretical predictions parallel 
the experimental results remarkably well; one of the sets of data relates to a sys- 
tem which is open-loop unstable. 

A discussion of scalar performance indices and some of their limitations 
is followed by a laconic survey of methods for solving the "best'' performance 
problem, especially for varying-parameter plants. The notion of a ''vector" per- 
formance index is introduced and a technique for establishing one response as 
"better'' than another is presented. This technique rests intimately upon the 
formulation of an ''artificial'' function which, when reduced to zero, identically 
satisfies the upper bounds on system requirements. These requirements are set 
forth in the ''vector'' index. The design for a ''best-overall'' compensator to meet 
specifications upon a varying-parameter plant is obtained by using a ''direct 
search'' method (especially suited to the electronic computer). This method in- 
corporates a ''pattern strategy'' to aid in finding the next set of trial values for the 
compensator; the strategy affords rapid convergence to a solution, as is demon- 
strated by an example. 

Finally, the example itself displays the versatility of the variable-para- 
meter mappings presented previously. Setting the human-operator transfer func- 


t 


tion equal to unity shows that the mappings have use for "inanimate" system design. 
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"Instead, we shall regard all that has been presented so far in this volume, 
and all that will follow, as merely a set of conceptual and mathematical tricks, a 
bit of scientific legerdemain, which, for reasons unknown, but for which we are 
profoundly grateful, prove to be startlingly more successful in many applications 


than anyone who examines the basic assumptions would have any right to expect." 
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CHAPTER ONE 


STABILITY AND PERFORMANCE OF MANNED CONTROL SYSTEMS 


ial introduction and Historical Background 


The unique capabilities of man, while seeming to pale in the face of 
increasingly sophisticated and automatic equipment, still offer characteristics 
either wholly unattainable by machines or, at the least, constrained in their 
attainment by weight, reliability, or economic considerations. Ultimately, the 
purpose of any good system, e.g., economic, political or control, is to ease 
the burdens of man, either directly (automatically serving his convenience, 
transporting him to some preferred location), or indirectly (aiding him in 
attaining his desires, whether the production of goods or the retention of vehicle 
control). Itis not surprising, therefore, to find that a good deal of study and 
experiment have been devoted to various aspects of human behavior. These 
range from investigations of a psychological variety to those of a less ambitious 
nature involving the characterization of human dynamics. Even in the latter, 
less extensive domain, any attempt at describing the rich diversity of human 
dexterity is an ambitious one. To permit mathematical description of human 
response to stimuli, the scope of the problem must be reduced (to a tracking 
task with visual inputs) and the input-signal characteristics constrained to cer- 
tain classes. Using such techniques, Tustin, working on the problem of man- 
operated-gun-fire-control equipment, first attempted to describe the human 


(1) (2) 


controller in servomechanism terms. About that time, Sobczyk and 


(1) 


Tustin, A., An Investigation of the Operator's Response in Manual Control of 
a Power-Driven Gun, C.S. Memo 169, Metropolitan-Vickers Co., Ltd., 
Sheffield, England, 1944. 


2 
( Pabezyk, A., Aided Tracking, MIT Rad. Lab. Rpt. 430, Sept., 1943. 





at 
1” 


Phillips’ performed similar work on "aided" gun-laying systems. . It was 
Tustin, however, who first arbitrarily separated the input-output characteristics 
of man into a linear operator term anda "remnant'' term, to account for that 
part of the human-output response not attributable to the linear operation on the 
maput, Fig. 1-1. Itis clear also from Fig. 1-1 that the system under considera- 
tion restricts the human controller to processing the error signal; this configura- 
tion is called "compensatory tracking", a term assigned by the engineering 
psychologists. Certainly, not all human control can be classified as tracking, 
nor does the tracking task imply human control. A sufficiently large group of 
situations is covered by the compensatory tracking mode, however, to make 
its study very worthwhile; even in those situations where, due to other complica- 
tions, this mode does not provide an accurate description of system conditions, 
an understanding of such manual-tracking capabilities can provide insight into the 
more complex situation. 

It is evident that a mechanism so adaptable that it can ''shoot marbles" 
or pilot high performance aircraft does not lend itself to description by a single 
transfer function. After World War II, several now-classic studies of the 
human controller were performed. eye) From these and the comprehensive 
work of McRuer and Krendel, there began to evolve a general form of human- 
operator-transfer-characteristic in which some confidence could be placed. (4) 
It became apparent that, in general, the particular transfer description was a 
function of both the controlled element and the input power spectrum. Thus, 
tn aided tracking, the human controller's output signal directly affects both 


the position and rate of the gun; thus, it is possible to simultaneously correct 
errors in gun-mount displacement and speed. 

1 
( Phillips, R.S., H.M. James, and N.B. Nichols, ‘"Iheory of Servomechanisms | 
McGraw-Hill, New York, 1947. 

(2) 


Russell, L., Characteristics of the Human as a Linear Servo Llement, 
hie. Ihesis, MIT, May, 1951. 


2) 
( elkind, Jele, Characteristics of Simple Manual ControMsyctenis, ao 


Thesis, MIT, April, 1956; also issued as MIT Lincoln Lab. Rep. No. I1l. 


4 
( McRuer, D.T. and E.S. Krendel, Dynamic Response of Human Operators, 


WADC Technical Rept. 56-524, Oct., 1957. 
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combinations of useful sets of plant dynamics and input power spectra were 
presented ina table with certain ''adjustment rules" which allowed the designer 
to 'set'' the parameters in the human-operator transfer function in accordance 
with empirical evidence. As a result of the environmental dependence of the 
human-operator transfer function and the occasional lack of point-to-point 


correspondence between the model and actual tracking records, the term 


"quasi-linear''is used to qualify the transfer function description; thus the 


model portrays a smoothed, time-averaged type of "linearity"'. Since the 


publication of the McRuer compendium, additional studies have demonstrated 


(1)(2) 


further the utility of the transfer function approach. The assumption of 
human operator stationarity | is inherent in most of the studies involving the 
transfer function approach, although non-steady behavior is, to some extent, 


Inveasured by the remnant term; work has been conducted to determine time- 


(3) (4) 


varying descriptions of the human controller. Although many elegant 


models have been presented in the last few years, none has been tested and 
established to ee extent of the model contained in the studies of McRuer and 
Krendel. This quasi-linear mathematical description, in its most general form, 
appears as 


Ke (itjeT |) 
H (jw) oat o.. <n (1-1) 
(14jwT (14527) 


L 


Thatis, during any given measurement run, one assumes that human response 
ieemot time dependent. 


] 
( all, I.A.M., Effects of Controlled Element on the Human Pilot, Princeton 


Warvy., WADC Tech. Rept. No. 57-509, Aug., 1958. 


Z 
( ISmith, Pete wOn tie limits of Manual Control, IDEE Urans. on. 


MaretOrs in Electronics, Sept., 1963. 


(3) 


Sheridan, T.B., Studies of Adaptive Characteristics of the Human Controller 
Tech. Rept. No. ESD-TDR 62-351, L.G. Hanscom Field, Bedford, Mass., 
Wec., 1962. 


+ 
ijer ville, W.W., G.A. Gagne, A Theory for the Optimal Determi mice 
Characterization of the Time-Varying Dynamics of the Human Operator, 
@ernell Aero. Lab., NASA Contractor Rept. CR-170, Feb., 1965: 





It is noted that the description contains a gain term, a pure lag 
(frequently associated with human reaction time), a pole with time constant Ty 
(taken to represent neuromuscular lag), and an equalization characteristic with 


time constants TL and T- This equalization capability represents the human 


operator's capacity to "'tune"! 


his response in much the same manner a control 
engineer compensates a system by introducing lead or lag networks. Considerable 
experiment has shown that descriptions more complicated than (1-1) offer little 
added accuracy in ''goodness of fit'' to the data; frequently, however, the form 
of (1-1) can be simplified by eliminating some of the time constants whenever 
they are not required to obtain a good-fit to actual human-operator response 
data. In part, this represents the analytical analog of the human-operator 
adaptive capacity. Equation (1-1) is obtained by analysis of power spectra 
records; the expression is derived under the assumptions of time-invariance 
and "'trained'' operators, i.e., the human-controller transfer characteristic 
represents post-learning period response-capability. This assumed familiarity 
with the task introduces sufficient response-consistency to permit mathematical 


description of the data. 


1.2 Human-Operator Adaptive Capacity 


The transfer function approach continues to offer the most promise to 


engineers involved with human-controlled systems. Studies designed to increase 


1)-(5 
the utility of this technique are still being Seemed. | )-(>) This thesis presents 


De iickler. F.A. and R.W. Obermayer, The Use of Man in Booster Guidance and 


Control, NASA CR-81, July, 1964. 


2 
( Wemeller, F.A. and R.W. Obermayer, Control System Lags and Man-Machine 


System Performance, NASA CR-83, July, 1964. 


3 
( ers. L.G., and K. Ziedman, A Study of Manual Control Methodology with 


mAmnotated Bibliography, NASA CR-125, Nov., 1964. 
(4) 


Mereeron, H.P., J.K. Kincaid, and J.J. Adams, Measured Human Transfer 
Functions in Simulated Single-Degree-of-Freedom Nonlinear Control Systems, 
NASA TN D-2569, Jan., 1965. 


5 
( Geer mayer, R.W. and F.A. Muckler, On the Inverse Optimal Control Problem 


in Manual Control Systems, NASA CR-208, April, 1965 
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a different approach to meeting the stability and performance requirements of 
certain man-machine systems. In the many studies conducted, there appear 
many different values for the parameters K,,T, I,,T_and Ty of (1-1); these 


h al 


‘character of the human-operator 


distinct values indicate the ‘'quasi-linear' 
transfer function; also they simultaneously imply the adaptive capacity of the 
human controller. From these empirically derived values, ‘''capability bounds" 
upon the human operator's adaptive range can be formed. Certain mathematical 
manipulations and a suitable method for displaying such bounds, when used in 
conjunction with plots of the controlled-element dynamics, permit stability 
estimates to be made for the combined closed-loop system. The assumption 
that the human operator regards system stability as a primary criterion is 
intrinsically contained in these estimates; it has been established in experiments 
that this is a valid supposition. Experiments also have shown that the human 
controller is most satisfied with the dynamics of a system when he is able to 


(1)(2) 


participate (and maintain control) as a simple amplifier. This property 


' remaining can serve as a measure 


implies that any ''margin of human capacity' 
of human satisfaction with the system ("handling qualities''); these margins 
appear as a part of the ''capability bounds" display. The bounds formulated and 
used in this thesis cannot be regarded as definitive parameter ranges for large 
populations since such do not now exist. Indeed, if past experiments share any 
common feature, it is the limited number of subjects used to provide data. 
Usually, the small range of variation among well-trained subjects is cited as 


justification for limited subject sample size. The cross-sectional aspect of this 


thesis, i.e., its use of the data from many experiments, to some extent 


Seeemningham, H.P. and F.V. Taylor, A Human Engineering Approach to the 


Design of Man-Operated Continuous Control Systems, NRL Rept. 43333, 
April, 1954. 


BiieRuer, D, 1. and E.S. Krendel, The Human Operator as a Servo System 


Element, J. Franklin Inst., Vol. 267, Nos. 5 and 6, May and June, 1959. 








meenuates the effect of reduced sample size. The ranges established foruihe 


parameters of (1-1) are as follows: 


KL >: "02, Gea aia0 (1-2a) 
7 : O.1 - 0.4 (1 -2b) 
Th 2 Ures5 26 (l=2c) 
T, 2 0 = 25 (1-2d) 
Ty : O-0.7 (1-2e) 


The units for K depend upon the units of the input-signal to the human 
controller, the others have units of seconds. It should be noted that the ranges 
on KL and tT do not restrict the analysis in any way; as is seen below, these two 
parameters enter into the method in sucha way that any range desired can be 


accommodated quickly. 
meome tine ‘Difference Vector" Plot 


Por the system of Fig. |!-1 the closed-loop system function is” 


Vv gee © 
Me OEP (1-3) 
R ae 

O 


The system is stable if the denominator contains no zeros in the right-half plane. 
It is noteworthy that in a monitoring configuration, _ Fig. 1-2, the closed-loop 
function is 

Pret Grae x 
- O 
oe (1-4) 
Ry ts GG 

a ee 


[= 


“Equation (1-3) does not consider the effect of n(t).. The stability analysis which 
follows does not depend upon this exclusion. 


In a monitoring configuration, the human operator observes system error and 
Sweets the loop as a ''signal processor’ only after he decides his contributions 
can improve system response. To include the "no entry" condition the lower 


bound on KL (l-2a) is reduced to zero. 
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Fig.1-2 : The Human Controller as a Monitor 





Again, the system is stable only if the denominator contains no right-half plane 


zeros. The denominator may be written 
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In (1-7), the term containing the variable-parameter portion of the system is 
i Ge (1-8) 
o Pp 


The term 


1+G G (1-9) 
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in (1-7) presents no serious problem in stability analysis since it has a fixed- 
parameter structure. It should be noted that (1-8) has the form of the denomi- 
nator of (1-3); of course, any effects from the poles of (1-6) remain to be 
considered. Expression (1-8), set equal to zero, is the characteristic equation 
for the system; for stability this equation must have no roots in the right-half 


plane. 
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Substituting (1-1) into (1-10) and manipulating, one obtains 
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In Eq. (1-11), if the first term is called, 9, and the second, L then (1-11) 


becomes 
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By an extension of the Nyquist criterion, ie the encirelements that ee makes 
of the ) region provide stability information about the system. 

Usually, ina fixed-parameter problem the system structure depends for 
its stability upon the nature of encirclements of the ''minus one'! point by its 
Nyquist contour; with a varying parameter structure, this ''critical point" 
generalizes into a ''region of critical points'', each of which is obtained from a 


specific set of values assigned to the triple (T,, T., Ty): Unfortunately in order 


ae | 


to display this "region of variation" for the ranges on T,,T_,,T 


i aN 
(1-12) must be plotted for particular frequencies y = W. Thus, to interpret 


the term ( of 


"encirclements'' and infer stability characteristics of (1-12), several sets of 

such ''regions of variation'’ must be plotted for pertinent frequencies. Fortunately, 
in the human controller problem only a very few frequencies are required as a 
consequence of the rapid deterioration, at frequencies above w20 radians per 
second, of both the human operator and the accuracy of the mathematical model. 

It should be noted that (1-12) 1s written such that represents the VARYING 
portion of the system and Le the FIXED portion. In plotting these to determine 
stability (discussed below), the effects of varying-human-operator gain is 

readily observed on the Nyquist curve for -L the human reaction time delay is 
handled by selecting pertinent values for T and viewing their overall effects. Any 


number of such T values can be included quickly and easily, since an added phase 


l 
( Fiche. Were iscussion, Joeppl Mech... pa aig, Mec., 1940" 
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shift to “Lo is all that is required. 

Rather than plot the sum (1-12) and observe any rotations about the 
origin, the following technique is used. The negative of Le iseplottec laces 
-L. as well as the region of variation 0; Fig. 1-3(a) displays a hypothetical 
Situation. If a vector is drawn from the ~L Gu) loeus to the @ “locus, each ena 


"is actually 


of the vector being at the same value of w, that ''difference vector 
the function : +L and the rotations it experiences give stability information for 
that function. Thus, stability is determined from the nature of the origin 
encirclements, if any, that occur in the "difference vector" plot, Fig. 1-3(b). 
The vectors shown in Fig. 1-3(a) are for specific values of the varying para- 
meters as well as for particular frequencies. Thus, the locus of (1-12) may be 
visualized without carrying out the summation. The replacement of Q by unity 
results in an ordinary Nyquist plot with the head end of the vector fixed at 

fer yo and the other tracing out the -L (jw) locus. Thus, as a dividend the plot 
of -L_ allows immediate observation of the plant-stability characteristics in 


those cases whenG =0; the critical point from which stability is determined in 
e 


S@eminstances is at +1 +)0. 


1.4 The Mapping Function and the Stability Criterion 
With the human operator's gain-varying capability and delay time lumped 
into the fixed-parameter structure, the varying portion of the characteristic 


equation is given by 


Cla a Peayel a yt —) 
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It is shown in Chapter Four that for fixed w, Typ and Th variations in TT, map 


into circular arc segments, while with the other parameters fixed, variations in 


T, or Ty map into straight lines. In this way, the bounded region of variation 


for the ranges selected (1-2, c through e), Fig. 1-4, may be plotted. A typical 
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plot appears in Fig. 1-5; this plot is a mapping through () of the boundaries (end 
points) of Fig. 1-4. Frequently, and in the case of (1-13), the "end point" 
mapping is not "complete". The property of completeness signifies that every 
point of the mapping Q, as the space of points in Fig. 1-4 is scanned, is also 

a point of the "end point'' mapping. This property is required if rigorous 
statements regarding stability criteria are to be made; i.e., completeness 
permits the establishment of necessary and sufficient conditions for stability. 
For human-controller studies, it may be shown that the mapping of Fig. 1-5 is 
"almost complete" and its application to engineering design studies is acceptable. 
Plots of Q for the frequencies w=0.5, 1.0, 2.0 and 3.0 are given in Figs. 1-6 
through 1-9; these plots are "“universal'' in the sense that they may be used with 
any system of interest. 

With the critical points governing stability defined by a region Q), 
estimates of human-operator handling capability are made in the following way: 
if there is some area of 2 unenclosed by the contour, “Lo then presumably the 
human operator can assume such parameter values and stabilize the system. 
The parameter values necessary to accomplish this appear to give insight into 
the amount of effort the human operator must expend to retain control. oe If the 
entire () region is enclosed, the human operator cannot control under any cir- 
cumstances (for the capability ranges assumed and plotted). If no portion of 
the () region is enclosed, there is stability no matter what the actions of the 


human operator (within the sense of the transfer function description). 


The ultimate justification for this approach is provided by predicting 
whether the human operator can control (stabilize) certain systems. The 


systems used in this thesis were analytically equivalent to those reported on in 


] 

( IMieruer, D.R., I. L. Ashkenas and C.L. Guerre, A Systems Analysis View of 
Longitudinal Flying Qualities, WADD 60-43, Wright Air Development Division, 
tan., 1960. 
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Fig. 1-6 : Variable Parameter Mapping 
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Fip.1-7 : Variable Parameter Mapping 
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Fig. 1-8 : Variable Parameter Mapping 
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the literature. The published results were determined experimentally; their 
classifications were obtained from human operator evaluations (pilot opinion 

of handling qualities). The a priori predictions parallel the empirical results 
very well, including some for systems which are open-loop unstable. 

Qualitative information concerning handling qualities and engineering insight into 
how the systems may be modified to improve their stability and the human 


operator's satisfaction with their characteristics is also provided. 


feo eipplications of the Theory 


In a simulator study by Ian Hall at Princeton University the effects upon 


(1) 


a pilot of various controlled-element longitudinal dynamics were investigated. 


McRuer et al performed a detailed analysis on the results of the Hall report in 


Z 
an attempt at developing a ‘‘flying qualities" ‘neo The McRuer approach 


follows the ''adjustment rules'' procedure described in Section 1.1. Another 
experimental study conducted at the University of Cincinnati by Ralph Smith 


investigated the maximum capabilities of a human in controlling certain open- 


(3) 


loop-unstable systems. The results of both these studies furnish experimental 


justification for the theoretical predictions in this thesis. 


As an indication of the analysis technique, a single example from the Hall 
study is offered in this chapter, other applications appear in Chapter Five. The 


experiment was conducted in the compensatory arrangement as shown in Fig. 1-1 


l 
( all, I.A.M., Effects of Controlled Element on the Human Pilot, Princeton 


University, WADC Tech. Rept. No. 57-509, Aug., 1958. 


2 
( ie Ruer, D.T., I. L. Ashkenas and C.L. Guerre, A Systems Analysis View of 


Longitudinal Flying Qualities, WADD 60-43, Wright Air Development Division, 
fan. , 1960. 


3 
( (ornith, R.H., An Experimental Determination of the Limits of Pilot Controll- 


ability for Unstable, Oscillatory Second Order Systems, Thesis for degree of 
Aeronautical Engineer, Univ. of Cincinnati, 1961 (also appeared as a paper, 
mathe Limits of Manual Control, IEEE Trans. on Human Factors in 
Piectronics, Sept., 1963). 
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with the plant (aircraft) dynamics having the generic form (short-period 
approximation): 
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Hall used twenty combinations of € and WO to test pilot response. The plot for 
“Lo = 1 from (1-14) for ¢€ =0.2, w= 1.57 1s shown in Fig. 1-10; this configura- 
tion was labelled C21 in the Hall report and is so labelled in Fig. 1-10. 

Vellum (or transparent) overlays of figures such as Figs. 1-6 through 
1-9 can be made and used to good advantage with system plots as presented in 
Fig. 1-10. Such overlays facilitate envisionment of whether the ''difference 
vector'’ between “Lo and () undergoes some type of rotation. In the case pre- 
sented, a glance at Fig. 1-10 and the general trend of the 2 mappings for the 
frequencies given (w=0.5,1.0,2.0 and 3.0), indicate that system stability is 
Metetimined by the relation between the “Lo plot and the mappings for w=2.0 and 
w=3.0. This situation results since the -LT plot leaves the ) region of 
criticality for wx2; this implies that a condition of instability may pertain. It is 
not yet clear how to take the difference vector between -L G2) and @({2); pie 
conceivably, it could be connected from -L (32) to any point within or on ()(j2); 
for instability, it must be demonstrated that a clockwise encirclement occurs as 
woo irrespective of the selected orientation of this vector (all the Hall configura- 
tions were open-loop stable). 

The time invariance hypothesis implies, however, that once a settled 
condition of operation has been adopted by the human controller, any ''selected 
orientation" of the difference vector defines the sub-region in space to be used for 
the remainder of the stability analysis. For example, if the human operator 
assumes a set of parameters near point B’ of Fig. 1-8 (corresponding toa 


Memenborhood of max T., max T,, max T,.), there results, from Fig. ]-10, an 


L I N 
encirclement by the difference vector and instability. If the point F’ is assumed, 
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Fig.1-10: Plot of Configuration C2) 
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the human can “hold out" longer before enclosure, but eventually an encirclement 
again occurs. Unfortunately, the prospects of being at F’ are not encouraging 
since, although the point corresponds to max T, (lead) which is obviously needed, 
the location F’ also implies Ty =0 which, if the physiological correlations are 
credible, corresponds to no neuromuscular lag. If there is any hope for system 
stabilization or controllability, the human operator must assume parameter 
values somewhere between arc G’O’ and arc A’D’ and in the vicinity of the 
origin of Fig. 1-8, also Fig. 1-5 (detail). Particularly, is this the case when 
the human time-delay is imposed upon the system. Thus, the characteristic of 
C21 projecting beyond the ( region is further aggravated by the introduction of 
the time lag, i.e., the point -L (je) located at approximately 5.0 on the real 
axis rotates, with the distance from the origin to 5+j0 as radius arm, a number 
of degrees proportional to the time delay. 

Since C21 represents a set of dynamics for which empirical evidence of 
human dynamic response exists, it is possible to check the inferences above. A 
Migest fat’ Keene operator description was calculated by Hall and a plot of this 
fienetion Can be made. In Section 5.3 it is shown that the human controller did, 
indeed, assume those characteristics previously implied by the analysis; further- 
more, itis obvious that his capacity is being strained severely just to maintain 
control of the C21 system. That is to say, few points in the variable mapping 
space Q are not enclosed by the plot of Fig. 1-10, particularly after introduction 
of time-delay effects. 

It is interesting that the "subjective evaluation" of this system by two 
experienced Navy pilots was a rating of 'unacceptable-quite unsuitable. '' Good 
correlation between the analytical predictions and the pilot ratings also occur for 
other sets of dynamics (Chapter Five). As implied already, correspondences 
between boundary lines of the mapping and the physical capacity of a human 


Gperator seem to be possible. 
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1.6 "Best'' Performance and ''Best-Overall'' Com ensation of Varying- 


Parameter Plants 


As in all engineering design, several performance measures have been 
formulated for man-machine systems to provide a so-called figure of merit or 
index of quality; typically, deficiencies have been found to exist with all these 
scalar indices. The difficulty in attempting to quantify the performance ofa 
complex system by a single number is that satisfaction of a requirement on this 
unique scalar may still result ina system with other unsatisfactory character- 
istics, irrespective of formulation and weighting techniques. A ''vector"' 
performance index might offer significant advantages. Such an approach to the 
performance problem is pursued in this thesis. Hopefully, it may be suggestive 
of a standard technique, or algorithm, for handling far more severely complex 
situations than can be exhibited here. Also, one feels intuitively that design 
compromise on an adequate compensator should be possible so as to provide 
"best'' performance even over some range of plant-parameter variation. A 
simplified example of a method for satisfying such demands is given below. 

The goal is a ''best-overall'' compensator for any given system with a 
human controller, Fig. 1-11. The approach may be used for a plant with variable 
parameters; the solution to be presented, however, depends intimately upon an 
initial stability study conducted in accordance with the procedure previously out- 
lined. Thus, in Fig. 1-11, the result desired is values for K a and @ which 
provide good performance for all possible values of ¢. This portion of the 
thesis departs from stochastic-input-signal considerations and studies 
deterministic-input-signal response. Thus, performance may be measured by 
the classical indices te on and t. (response rise-time, overshoot and settling - 
time resulting from a unit step input); these follow exactly from some arbitrary 
definition and may be regarded as a vector quantity in three space. Asa result 


of a stability analysis as outlined above, approximate values for Koo a, and 8 
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conducive to stability are determined. In addition, the designer attempts to use 
the freedom in selecting a and 8 (for this imposed class of compensator) so that 
the plant's Nyquist curve passes through that part of the Q region co~responding 
to reduced task requirements on the human operator; i.e., one establishes as 
large a "margin of human capacity''as possible. With the ranges of q and 8 that 
provide acceptable stability available, nominal values for the human controller 
can be read from the plot using the time-invariance assumption. Simultaneously, 
these values must retain the "high-frequency'' end of the compensated Nyquist 
contour in a location that continues to satisfy the stability criterion, as indicated 


by the corresponding (© region. * Thus, the triple (TD > T TN) must have value 
such that it lies to the left of the plant's Nyquist curve (traversed in the direction 
of increasing w) and still be within the range of human adaptivity. Also, 

Th» 7. and Ty should be as small as possible (the human prefers to operate as a 


simple amplifier), particularly T since this parameter is almost a direct 


L’ 
adaptation. 


maeweation of ''comfortable'' 


The above analysis, conducted for a nominal €, is repeated for several 
values covering the range of €. The greatest common range of a and §, obtained 
from the various ranges resulting from the repetitive analysis, then provide 
satisfactory stability over the domain of varying ¢. From these final ranges of 
acceptable q and 8B, the selection is made of those values of a and 8 which produce 
the best performance. There is no difficulty in determining te O ay for any 
given set of KO a and 8; the problem is how to select subsequent trials of q and 8 
from the final range established, so that each iteration is at least as good as the 
last. 

the function 


f=(t -T eto -O oe -T i (1-15) 
Psa ss ss 
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This statement assumes the human operator uses an equivalent performance- 


criterion. 
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is formed, where 


That is, as Oe and a are the maximum acceptable bounds on a 0. and te The 
function is an artifice created specifically to facilitate finding some point on the 
surface of acceptable performance; at f=0, the specifications are met identically 
at their upper bounds. A minimization process conducted on the function (which 
contains time-domain parameters) produces useful information regarding 
acceptable Ko a and 8 (which represent frequency-domain parameters). 
Appropriate reductions in TO ee and T. at timely points in the search procedure 
provide resolution of corresponding boundaries of acceptable-performance in 
both (t > 00 i) and (K a, &) space. 

A very straightforward search technique is used to find a compensator 
satisfying the requirements. (1) This method of "direct search" with "pattern 
strategy’ is especially suitable to minimization of a sum of squares, but its 
overriding virtues are its straightforward simplicity and compatibility with the 
lecrmonic computer. ''Direct some Be means a step-by-step examination of trial 
solutions and a method of comparing eachto determine which is "best'' in some 
sense. Additionally, a method is afforded for choosing the next trial point. This 
last aspect forms the strategy part of the search procedure; a ''pattern'' approach 
is used and from this is derived the strategy's name. 


The simplified example used to demonstrate this search technique 1s 


(1) 


Hooke, R. and T.A. Jeeves, Direct Search Solution of Numerical and 


Statistical Problems, J. Assoc. Comp. Mach., April, 1961. 
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shown in Fig. 1-12; the human-operator description has been set at unity since 
its presence, while having no effect upon the search procedure, complicates the 
explanation. It is assumed that the ranges of q and 8 have been esteblished which 
offer good stability; specific values within that range are now desired which meet 
the step-response requirements on rise time, overshoot and settling time. For 
economy of notation K. is redesignated K. The initial "base point" of the search 
(Ky, ay» 8) is chosen arbitrarily or the designer may exercise a "hunch". Again, 


the function, f, (1-15) is formed; the point (KX, Cas 8.) is said to be "better" than 


Z 


(K),a),8)) if 


f(K »05,8,) < (K,,a),8)) (1-16) 


Z 
By f(K., aL, » 8.) is meant the value of (1-15) found from the substitution of the te 
oe and “a associated with the system step-response when the compensator has 
gain and time constants, K., Cl. » and 8. respectively. If (1-16) holds, a "nothing 
succeeds like success'"' philosophy dictates a continuation of the search in the 
same ''direction''; therefore, exploration is continued by perturbing the inde- 
pendent variables with increments whose proportions are based on the frequency 
of success achieved along each dimension, K,a,8, in the recent past. 

In the system of Fig. 1-12, the specifications are chosen to be 


T == seconds, O, = 15 percent, T = 4 seconds (1-17) 


also, the plant variations possible are given to be 


Onl 


aN 


=a (1-18) 
Exploration is conducted using step sizes as follows: 


For K_ : take step 5, = 0eec (1-19a) 
Cc 
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Fig. 1-12 : Compensation of a Varying- Parameter System 
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For q : take step 5, = 0.1 (1-19b) 


For & : take step 6. = 0.05 (1-19c) 


The initial base point is chosen as 


K, =o (1-20a) 
a = 2 (1-20b) 
8, = O25 (1-20c) 


In (1-19a,b,c) the subscripts 1, 2,3 are assigned to represent step 
increments along the K a and 8 dimensional axes, respectively. Usually, the 
problem and design experience suggest step sizes for each independent variable. 
During the termination procedure, the nominal step sizes are reduced until 
impractical as a consequence of solution indifference or equipment limitations. 

The search is started at the varying plant's ''sluggish end" (€ =1.0). The 
initial base point for the searchis called, by (1-20a,b,c); the horizontal bar 


signifies the vector (multi-dimensional) nature of the point, i.e., 
b, = (1,2, 0.5) ean 


The step response of the system with a compensator which possesses such 


@haracteristics is 


(t ,o ,t ) = (4.8,0, 9) (1-22) 
| & S Ss 
Prom (1-15) and (1-17), 
2 2 2 
file e0s 5r= L422) 410-15) £lSe4) =25G 51 (1222) 


This value for f must now be compared with the value at the next exploration 


point. The next observation is performed at the point b, +8, where 


o| 


een (1. 2,2, 0.5) (1-24) 








a0 


The response is found to have performance measures 


(t_,0.,t.) = (3.5, 0, 6.8) (1-25) 
Therefore 
2 2 2 
£22, 0.5) = (3.5-2] POSS) lon s-4)) eee (1-26) 
Thus, 
(ee as a5 (V2) 


and b, a becomes the new temporary base point, te from which exploratory 
moves along the q axis are conducted. 

The process is continued in this way, the best of three exploratory points 
(only two are shown above, generally, fis also evaluated at b, -5,, also) always 
leading to a new temporary base point. The procedure is performed along each of 
the axes (K_, a, 8) and yields a new base point, Cee after the completion of sucha 
searchcycle. With a successful ''direction" of search established (b,-b,), the 
next temporary base point used for exploration purposes is found by doubling the 
magnitude b,-b, and proceeding in the same direction. This technique forms the 
"pattern strategy" portion of the search procedure since future points of 
exploration are determined as a function of the frequency of success achieved 
along each dimension in the perturbation steps such as (1-24). Great strides 
can be made in pattern moves after very few iterations. This example was solved, 


effectively, after six iterations; the geometric aspects of the first two iterations 


are shown in Fig. 1-13. It is found that the triple 
(K_ 5a, B) =7(5, 170570205) (1-28) 


satisfies the requirements 


T =1.0, 6 =10, T =2 (1-29) 











Pig. 1-13 : "Direct Search" with "Pattern Strategy'' 


Only the first two iterations are shown 
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over the range of 0.01 <C¢ < 1.0; this is superior performance to that originally 
specified. 
Empirical justification of the above approach (with a human operator in 


the loop) can be offered only after additional research is conducted for deter- 


ale 
~ 


ministic inputs driving controlled elements other than pure gains. Itis 
recognized that the bulk of the evidence justifying the form of human-operator 
transfer function used in this study has been accumulated for systems driven by 
stochastic input signals. There is no reason to suppose significant changes in 
Meetorm for deterministic inputs; the initial range selections for K, ,7,7_,, 17 


h L I 
and T,. were quite broad, but modifications, if needed, are easily incorporated. 


N 
Other search techniques might also be used to identify those values of 

Ko a and 8 yielding ''best-overall'' compensation. Any method employed, 
however, must be suitable to finding the optimum of functions which are unknown 
to the observer; i.e., the analytic form of the function is unknown or too 
complicated to be used. This was the case in the previous example; the direct 
functional relationship between (Ko a, 8) and (t O t) is not easily specified. 
First, the transient solution for given Ko» a, 8 must be obtained by performing an 
inverse Laplace transformation. From the particular transient response resulting 
the values of ae oe: ue can then be determined. No direct relation between the 
time-domain performance-indices th oe t and the frequency-domain compensator 
parameters Se a,8 can be written. The end result is that experimental 


(1) 


optimization procedures become most attractive. Wilde describes several 


search techniques that could be used alternatively, e.g., Method of Steepest 
Ascent, Contour Tangents, Parallel Tangents or Method of Rotating Coordinates. 


ale 
PAS 


“Such results could not be located by the author. 


(wilde, D.J., "Optimum Seeking Methods, '' Prentice-Hall, Englewood, 


Cliffs, New Jersey, 1964. 


. 
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1.7? Summary 


Chapter One provides an overview of the main points considered in this 
study. The goals, accomplishments and hmitations of the investigation are 
sketched briefly and a concise description of several of the major techniques 
used is presented. Further detailed discussion may be found in the main body 


of the thesis. 
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CHAPTER TWO 
THE HUMAN OPERATOR IN CONTROL Sys tems 


2.1 Introduction 


The possibility of describing human operator 'behavior'' in terms of 
mathematical relationships began in earnest sometime during World War II. 
Confronted with the problem of designing complex systems such as gun turrets 
capable of tracking fast and maneuverable aircraft, engineers found themselves 
needing considerably more quantitative information about the human operator's 
transfer characteristics than was previously required in weapons design. Whole 
teams of professionals in mathematics, engineering, psychology and physiology 
set to work to determine more adequately the bounds of human performance and 


(1) 


the capabilities that could be expected. In addition to Tustin's work, efforts 


in the United States were being conducted along similar lines by Phillips; some 


mention of this is made in the now classic text of which he is a co-author. 


(3) 


The interested reader may find a good historical sketch in Sheridans's paper 
and the background interplay between the engineering and psychological philoso- 


phies given by McRuer. a The resurgence of interest in the human controller 


(5) 


probably may be marked by Russell's paper and the definitive work by 


|) rs 
, oT A 


ustin, A., An Investigation of the Operator's Response in Manual Control 
of a Power Driven Gun; C.S. Memorandum No. 169, Metropolitan- Vickers 
Electrical Co. Ltd., Attercliffe Common Works, Sheffield, England, Aug., 
1944, 


(2) James, H.M., N.B. Nichols and R.S. Phillips,"’;- Theory of Servomechanisms, 
McGraw-Hill, 1947. 
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Elkind. (1)(2) Military high performance aircraft, space capsules, deep sub- 
mergence vehicles and other systems whose dynamics are under the control or 
partial control of a human operator have served to intensify the interest in the 
human as a controlling element; additionally, the need for greater understanding 
of the man-machine interface has spurred research in this area. oe 

Almost all the human operator studies conducted to date, and this thesis 
as well, address themselves to the problem of a human controller in a tracking 
situation; in particular, these investigations have been done in the ''compensatory" 
mode, i.e., with the controller operating off the system error, Fig. 2-1. This 
topology lends itself most readily to analysis techniques; in the face of consider- 
able ignorance about the human mechanism and the fantastic description-defying 
complexity of that mechanism some simplifying assumptions are justifiable. 
Particularly is this the case when the simplifying assumptions still permit useful 
conclusions to be derived which succeed in predicting empirical results. 

Several characteristics peculiar to human-operator control are not treated 
in this study. Psychological aspects of human behavior (that is, stress, com- 
petitive situations, rewards, other similar task influences and their effects upon 
tracking) are not considered as such. Physiological loading influences (such as 
fatigue or alcohol) also have not been incorporated. This is not to discount the 
importance of these factors upon perceptual-motor behavior and their resultant 
effect on any mathematical model of the tracking situation; again, as many 


superfluous variables as possible must be eliminated if useful analytical models 


Tracking Response Characteristics of the Human Operator:, Human Factors 
Operations Research Lab. Rpt. No 40, PRP-2, Sept., 1953. 


i) Elkind, J.I., Characteristics of Simple Manual Control Systems, (Ph. D. 
Thesis, M.I.T., Apr., 1956) also issued as M.I. T. Lincoln Lab. Report, 
No; 111. 


(3) Hall, I.A.M., Effects of Controlled Element on the Human Pilot, Princeton 
Univ., WADC Tech. Dept. No. 57-509, Aug., 1958. 


(4) Bekey, G.A., An Investigation of Sampled Data Models of the Human Operator 
in a Control System, (Ph. D. Thesis UCLA, Jan., 1963) also issued as 
ASD-TDR-62-36 WADC, May, 1962. 
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of the human controller are to be formed. The completely successful explanation 
of human response in this situation of reduced complexity (a milestone still not 
yet attained) places the future investigator closer to the possibility of including 
these additional parameters in a more accurate description of gene1al human 
behavior. 


There have been recent studies attempting a time-variant description of 


(1)(2)(3) (4) 


Almost no work has been published on man as a monitor although the problem is 


humans and some experimental work relating to the adaptation process. 


apparently under active study. Much of the published material on the monitoring 


(5) 


situation has been essaylike' ', although there has been some hardware generated 


by the Federal Aviation Agency as part of the supersonic transport program. This 


system, known as ''Force Wheel Steering", allows the pilot to override the auto- 


(6) 


matic elements; flight tests have shown considerable promise. 


(7) 


A recent paper 


on monitoring reports a tentative theory to explain some experimental results 


but the theory is quite preliminary and, although some model fitting 1s attempted, 
nothing to compare with the demonstrated validity of the quasi-linear transfer 


function is produced. 


Operating as a monitor, the human controller observes the character of the error 
signal, but instead of remaining in the loop and processing the error signal con- 
tinuously, he enters only as frequently as he believes he is required to augment 
normally ''automatic’’ control. 


(1) Sheridan, T.B., Studies of Adaptive Characteristics of the Human Controller, 
Air Force Systems Command ESD-TDR 62-351, Dec., 1962. 


(2) Elkind, J.I., and D.M. Green, Measurement of Time-Varying and Non-Linear 
Dynamic Characteristics of Human Pilots, Aeronautical Systems Division, 
een. Rpt. No. 61-225, Dec., 1961. 


Elkind, J.I., et al, Evaluation of a Technique for Determining Time-Invariant 
and Time-Variant Dynamic Characteristics of Human Pilots, NASA TN D-1897, 
May, 1963. 


Young, Green, Elkind and Kelly, Adaptive Dynamic Response Characteristics 
of the Human Operator in Simple Manual Control, IKEE Transactions on 
Human Factors in Electronics, Sept., 1964. 
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(3) 


(4) 


Hickey, A.E. and W.C. Blair, Man as a Monitor, Human Factors, Sept., 1958. 


Control-Display Pilot Factors Program, Instrument Evaluation Project Report 
63-1, USAF Instrument Pilot School, Randolph AFB, Texas, Dec., 1963. 


(7) Senders, J. W., The Human Operator as a Monitor and Controller of Multi- 
degree of Freedom Systems, IEEE Transactions on Human Factors in Elec- 
tronics, Sept., 1964. 
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The following definition of quasi-linearization appears in Elkind. oe 
“Quasi-linearization is the process of representing the characteristics of a non- 
linear device by linear transfer functions whose parameters depend upon the 
environment of the device. These transfer functions are called quasi-linear 
transfer functions because they behave linearly when the set environmental para- 
meters remain fixed and nonlinearly when these parameters are allowed to change.'"' 

"Quasi-linearization is particularly suited to the kinds of nonlinearities 
observed in the human operator. His most significant departures from linearity 
occur when the input statistics or system dynamics change. Since its parameters 
depend upon the input and the system, a quasi-linear transfer function can handle 
this kind of nonlinearity. "' 

It is interesting to quote the following comment from a recent NASA report. 
The potentials of quasi-linear models have not been fully investigated and further 
studies along this line are warranted. (2) 

In this paper a technique is developed to display graphically a range of 
human-operator capability which has been derived from a variety of such quasi- 
linearization studies; these plots are then used, in conjunction with standard 
Nyquist contours which represent the controlled elements, to predict human 
control capacity. It is shown that these predictions parallel previously published 
experimental results. Following the stability studies, there is a chapter devoted 
to performance considerations, Chapter Six. The classical measures of step 
input response characteristics, rise time, overshoot and settling time are used 
to demonstrate that a ''vector'' performance index can be formed and satisfied. 

It is shown that the required compensator design for the satisfaction of such 
specifications can be determined using computer search methods in combination 
with an “artificial'' function which, when minimized to equal zero, identically 


meets the stipulated response requirements. 


; Elkind, J.I., Characteristics of Simple Manual Control Systems, (Ph. D. 


Thesis, M.I.T., Apr., 1956) also issued as M.I, T. Lincoln Lab. Report, 
Ivo, lil, 


(2) Summers, L.G., K. Ziedman, A Study of Manual Control Methodology with 
Annotated Bibliography, NASA Contractor Report CR-125, Nov., 1964. 
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eee, Controllability and Adaptation 


The utility of the quasi-linear transfer function approach to engineering 
system design has been established many times,over recent years; even in the 
face of its many acknowledged limitations it represents a technique which provides 
useful results. This thesis uses as its starting point the quasi-linear "best-fit"! 
functions and the empirical results that have been published in the literature. A 
fresh look at the analytical tools extant and some novel manipulation of system 
transfer functions afford a new method for incorporating the entire adaptive range 
of the human operator and his pure time delay in a stability analysis. To date 
servo analysis techniques used in man-machine studies have utilized a nominal 
transfer function where the specific values for the time constants are selected to 
reflect the gross form of the process, i.e., human-operator transfer functions 
for particular dynamic elements (plants or processes) have been experimentally 
determined and presented in a ''table"'’ which is then used by the systems designer. 
The table may be regarded as offering ''adjustment rules'' indicating how to ''set"’ 
the parameters in the human-operator transfer function for varied control tasks 
and input power spectra. It should be noted that these studies normally employ 
simple transfer function approximations, usually of a bilinear form, for the pure 
time lag introduced by the human operator. 

Some of the ''inappropriate'’ and lightly justified aspects of the quasi-linear 
modeling techniques are subdued in this ''new" theoretical approach to human- 
operator studies since a broad range of human-operator adaptation and control 
attributes are inherently included in the analysis technique. This broad range 
tends to reduce the sensitivity of analyses to inaccuracies in the nominal values 
for the human-operator parameters, which previously had to be specifically 
designated. Since no new experimental work is presented there is no attempt to 
put forth a new model or model matching technique; rather, the investigation 
should be regarded as providing another means of viewing the human operator in 
engineering terms with the physical justification of the theory obtained from pre- 
viously published empirical results. 

The mathematical model required to adequately describe the human operator 


for different tasks necessarily varies in correspondence with the human's actual 
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adaptive capabilities. To date, this variation has been characterized by selecting 
specific values for the quasi-linear description as a function of the particular 
task. The method offered in this thesis is sufficiently encompassing of the broad 
range a mathematical model of the human operator must take, that it compensates 
for some of the inherent rigidity this selection procedure introduces. The method 
accomplishes this compensation by its complete display of human-operator capacity. 
Simultaneously, the method provides both a means of analysis and a channel for 
developing better intuitive understanding of human-parameter effects upon a con- 
trol system and system effects upon the human. 

For man operated systems the most important questions to be posed are: 
"Can the human 'handle' the system?'' and ''What can be done to ease his control 
burdens?" A comprehensive analysis technique for certain man-machine systems 
is constructed by capitalizing on studies already conducted and results previously 
published in the literature. Itis, perhaps, a first attempt at a theory for coming 
to grips not merely with the human controlled system, since this, indeed, has 
been tried before, but with that system under the added complication of the 
human-operator-adaptive range. Thus, to answer the question, ''Can a human 
operator stabilize this system?'' would be useful. The method to be offered will 
allow estimation of whether a human operator can stabilize any given system 
assuming he can adapt his ''parameters'" over certain ranges; in general, the 
investigator need not select specific values for these parameters. Insight into 
the relative human effort required to stabilize is also provided. No approximation 
to the pure delay introduced by all human operators need be made. The ultimate 
justification for the new approach is then provided by predicting whether the 
human operator can stabilize certain systems. The systems used were analytically 
equivalent to those reported on in the literature, along with their associated human 
operator evaluations; the published classifications were determined experimentally. 
The predictions parallel the empirical results remarkably well even for systems 
with open-loop instability. Qualitative information concerning handling qualities 
and engineering insight into how the systems may be modified to improve their 
stability and the human operator's satisfaction with their characteristics is also 


provided. The theory offers a convenient and illuminating display of the effects 





41 


upon the system of the various parameters which form the human-operator trans- 
fer function. Thus, it becomes easier to observe which ones (and their associated 
physiological characteristics, Section 2.5) are causing difficulty or limiting con- 


trollability. 


2.3 Thesis Organization 


The dissertation is effectively in two parts. The first (Chapters Two-Five) 
deals with system stability and controllability by the human operator and has been 
described qualitatively in the above paragraphs. The second part is concerned 
with system performance and introduces a type of ''vector'' performance measure 
which could have considerable usefulness in future human operator studies. This 
portion of the thesis departs from the stochastic input signal and studies the rami- 
fications of deterministic forcing functions. Much experimental work in the psy- 
chological literature pertains to step input tracking and the techniques described 
in Chapter Six may have direct application there. Too often, however, the psy- 
chologist's reports are qualitative and include varied methodologies thus making 
cross-experimental comparisons with the engineer's work almost impossible. 
With all its difficulties the engineer's approach at least provides a consistent 
language. Chapter Six is offered not as a solution to the performance problem but 
as a stimulant to further study of performance measures. In conjunction with this 
"vector'' type performance criterion, a ''direct search" method is used to obtain 
a specific compensation of a given class which provides "good overall" performance 
in the presence of a plant (process) having variable dynamics. The results are 
applicable to automatic control systems, in general, and systems with a human 
operator, in particular. 

Additionally, in this last chapter, the displays of human operator capacity, 
offered in the first part of the thesis are shown to have utility in the search for a 
solution to the performance problem; indeed, it is demonstrated in an example 
that these displays can be used as a tool in the design of compensators for 


"inanimate'' systems. 
2.4 The Quasi-Linear Model 


In general, human-operator response characteristics are strong functions 
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of the input signal form. For signals such as single sine waves the operator tends 
to ''synchronize"' his response with the input; such situations where the human may 
exercise his capability for prediction have been called ''precognitive'’. Precogni- 
tive response can be considered a form of open-loop control rather than tracking, 
although certain tracking tasks do permit this type of behavior. If both the opera- 
tors corrective responses and the input are displayed separately, the form of 
tracking is called "pursuit"; the classic example of such tracking is aiming a 

rifle at a moving target through open sights. In ''compensatory'' tracking the 
operator cannot distinguish between the effects of his response actions and changes 
in the input, except under zero input conditions, i.e., the operator views an error 


compensatory’ tracking is aiming a 


display only. The analogous example for 
rifle at a moving target through a telescopic sight with a narrow field of vision. 
As mentioned above, this study restricts itself to the compensatory case. 

A very brief derivation of how the quasi-linear human-operator transfer 


function is obtained, follows. From Fig. 2-1, 
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where N(j) is the transform of n(t), a term arbitrarily added at the operator's 
output to account for that portion of the output not described by the linear transfer 
characteristic; this term is called ''the remnant. " 
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where the asterisk signifies the conjugate and the subscript T indicates the 


Fourier Transform has been taken on the associated truncated time function, e.g., 
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atts ulti 


u(t) = 
Oren Gee 
and a 
Uj) = J upltle 3? Tar 
Thus, 
> din vl fee Cale + Gu ee 
Uy 7 Cease | eG: I+H G 
J 
a l [ i l l a 
= ~~ lim 5% [RR ]+——-—~— lim 5 [NR ] 
mae rN oF To © 20 Pt oan oe gue 
the quantity 
lim ae 
T« Zl) 
represents the input power spectrum cae and the quantity 
tim NRO 2g 
T7 @ 
under the assumption that the signal and noise are uncorrelated. Therefore, 
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The amplitude and phase characteristics of the human operator model are then 


found from b 
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l.e., the cross-spectrum between input and operator output divided by that 
between input and error yields a set of data in complex form which represents the 
human operator block of Fig. 2-1 (solid lines). 


These data can be fitted closely by a linear analytic relation of the form: 


Keo +jwT, ) 
H (jw) = . 2-1 
od (14+j0 TL 1+jWT, (cm) 
Jackson has shown that negligible increased accuracy results if (2-1) is 
made any more complicated. o Equation (2-1) provides a ''best-linear-fit" in a 


minimum mean-square error sense assuming random inputs of a Gaussian type, 


(3) 


i.e., normal amplitude distribution. The human operator adjusts his para- 


meters to some unknown internally-generated criterion; comparisions with the 


mean-square-error criterion seem to support its validity as a substitute analytic 


(4) 


index. The Gaussian assumption has proved to be reasonable and realistic in 


the light of the low-pass nature of most control systems, i.e., low-pass char- 


acteristics in feedback systems tend to make the amplitude distributions Gaussian 


)(5)(6) 


(central limit theorem Additionally, the Gaussian signal lends itself to 


This is not so restrictive as might appear since any criterion which is a 
monotonically increasing, even error function is also minimized by the resulting 


(2) 


linear fit. 
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Jackson, A.S., Synthesis of a Linear Quasi-Transfer Function for the Operator 
in Man-Machine Systems, Wescon Convention Record, Sept., 1958. 
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Sherman, S., Non-Mean-Square Error Criteria, IRE Trans. on Information 
piheory, Sept., 1958. 
Booton, R.C., The Analysis of Non-Linear Control Systems with Random 


Inputs, Proc. Symp. on Non-Linear Circuit Analysis, Polytechnic Institute of 
Prooklyn, N. Y., Apr., 1953. 


Roig, R. W., A Comparison Between Human Operator and Optimum Linear 
Controller RMS-Error Performance, IRE Transactions on Human Factors in 
Electronics, Mar., 1962. 


Gibson, J.E., ''Nonlinear Automatic Control'', McGraw-Hill, New York, 1963. 
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(6) Arthur, G.R., A Note on the Approach of Narrow~-Band Noise after a Non- 


Linear Device to a Normal Probability Density, J. Appl. Phys., Vol. 23, 1952. 
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experimental work; also, Booton has shown that ''the best mean-square 
error approximation to any nonlinear devices having a Gaussian input will 


be realizable and is determined by (2-la)". 
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Greatest confidence in the description (2-1) occurs when it represents 
the human operator "adequately", i.e., when n(t), the remnant term, is 
small. A figure of merit termed "linear correlation’ has been generated 
which relates the amount of linearly correlated operator output to total 
operator output. From (2-2) is found the fraction of the linearly correlated 


operator output power spectrum to total operator output power spectrum 
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1) @ikind, J.1., f Simple Manual Control Systems, (Ph. D. 


Characteristics o 
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If the square of the linear correlation is written 
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Of course, the closer p is to unity the greater the confidence in (2-l) asa 
description of human-operator response characteristics. Frequently, the 


y (1) 


linear correlation is good (p 70.8 thereby justifying use of the quasi- 
linear model for low-frequency conditions. Further discussion of remnant 
characteristics and possible sources for the remnant term can be found in 


McRuer. (2) 


2.5 Adaptive Ranges of the Human Operator 


From various published data in the literature 


(5 )(4} 


the values assumed 
by Ky 


of these values ranges for each parameter were established. These data 


we Ty» qT, and Ty were assembled and from the largest and smallest 


were associated with controlled elements (system dynamics) having forms of 
pure gain, lags, integrators and second-order systems. Even simululated 
approximations to human operator control of the short period longitudinal 
dynamics of an aircraft and the dynamics of an F-80A aircraft engaged ina 

tail chase have been studied. In general, the most complicated human-operator 
transfer functions that could be justified by these data have the form (2-1). 


Some double leads and lags were used in the linear fits by Pegi! but these 


4 A IL., D.T. McRuer, A Theory of Handling Qualities Derived 


shkenas, 
from Pilot- Vehicle System Considerations, Aerospace Engineering, Feb., 
iyo7, 


val McRuer, D. T., and E.S. Krendel, Dynamic Response of Human Operators 
Wembe- Tech. Rpt. 56-524, OCect., 1957, 


(3) hia 


(4) Hall, I. A. M., Effects of Controlled Element on the Human Pilot, Princeton 
Unive, WADC Tech. Rpt. £7-509, Aug. Toc: 
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were isolated cases and they are of no great concern here. Again, the 
human operator adjusts his parameters to conform to each control task; the 
tasks themselves are changed by modifying the controlled element or the in- 
put power spectrum. 

Since a considerable amount of work has been expended in obtaining 
these data and they are highly indicative of human-operator capabilities they 
should be sufficient to permit upper and lower bounds to be drawn upon the 
human adaptation process. It is not maintained that the ranges given below 
definitively bound human-operator capabilities, but they are suggestive of 
limits obtained for a variety of forcing functions and controlled elements. 

The question as to whether a human operator can freely adjust these para- 
meters independently to optimize is not answered although some studies seem 
to imply he does. : The "adaptive behavior'' of the human operator refers to 
responses to tasks AF TER the learning or training period has been accom- 
plished. Inherent in the final adaptive-range potential is the assumption that 
the operator has gained some familiarization with the system. Studies are 
currently being conducted on the learning or adaptive process; some results 


(2)(3) 


have already been reported in the literature. 
Thus, the human-operator transfer function with the range of variation 
indicated below should provide a very broad, useful set of results until fur- 
ther definition of bounds is empirically determined. The validity of the ap- 
proach is not conditioned upon the absolute accuracy of these ranges which 


always can be modified when desired. 


: Roig, R. W., A Comparison Between Human Operator and Optimum 


Linear Controller RMS-Error Performance, IRE Transactions on Human 
Factors in Electronics, Mar., 1962. 


(2) Moone. L,Kk., D.M. Green, J. 1. Elkind, and J. A. Kelly, Adaptive 
Dynamic Response Characteristics of the Human Operator in Simple 
Manual Control, IFEE Transactions on Human Factors in Electronics, 
Sept., 1964. 


) Gibson, J.E., et al, Philosophy and State of the Art of Learning Control 
Systems, CISL, Purdue University, Nov., 1963. 
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The “best fit'' (in the minimum mean-square sense) human-operator 
transfer function based upon Gaussian input spectral analysis takes the gen- 
eral form (2-1) , 

oe 
ee ne (1+ jw T,) 


2-1] 
O (T+ joT 0+ jwT) ( 
The ranges on the various parameters are 


K 27086 =—250 


7) - 2, 20e - 0.4 
Ty: 0 - 5,3 
T, © 07025)4..25 
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where the units for KR depend upon the units of the input to the human-opera- 
Permeolock and the units for T, Th» Tp and TN are in seconds. 
_To date, little work has been done on the monitoring problem, and few 


of the techniques of control system analysis have been applied to either multi- 


axis Or monitoring tasks involving the human operator. (2) The stability 


analyses presented below are also applicable to the monitoring situation, The 
analysis is conducted to establish stability conditions only for the situation 


where the human enters the system and remains as acontroller. Further 


Some work has been done on sampled data models of the human operator 
based upon the philosophy that human performance is of an intermittent 
nature rather than continuous. This concept is not used in this study; 


ab) 


details are available in the literature. ° Bekey's work also contains a 
delightful critique on the psychological literature. 


te Bekey, G.A., An Investigation of Sampled Data Models of the Human 


Operator in a Control System, (Ph. D. Thesis UCLA, Jan., 1963) also 
issued as ASD-TDR-62-36 WADC, May, 1962. 


Dander, V.A., Predicting Pilot Ratings of Multi- Axis Control Tasks from 
Single Axis Data, IEEE Trans. on Hum. Factors in Electronics, Sept., 
| Shey eh 
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work on periodic entrance and exit could be useful in application to time- 
sharing a human operator, i.e., the human is responsible for N systems 
which he may "'trim up" ona "multiplexing'' basis. As an amplification of 
the sense of the word ''monitoring" as used in this study, the following quote 


(1) 


is offered. "Force wheel steering operationally maintains the human 
pilot as an active control element with an automatic system for the purpose 
of having the human pilot monitor, supervise, and act as a redundant element 
for those contingincies which require pilot action... THE AUTOPILOT ACTS 
ON THE SAME GUIDANCE INFORMATION THAT IS PRESENTED TO THE 
PILOT; therefore, the pilot can easily monitor autopilot performance. For 
those approaches when the autopilot is operating satisfactorily, the pilot 
makes a minimum number of inputs with the Force Wheel. But the pilot can, 
upon assessment of the total situation, modify the autopilot performance using 
the same information as the autopilot by making normal control actions... 
pilot inputs can be made WITHOUT disconnecting the autopilot.'' A possible 
block diagram display of such a monitoring system is shown in Fig. 2-2. In 
order to incorporate the "hands off'' mode in monitoring, i.e., when the 
human controller effectively eliminates himself from the loop as during the 
"off duty'' period in true monitoring, the lower bound on KL has been reduced 
to zero. To facilitate certain manipulations, and ease the algebra in some of 
the mapping theorems, the lower bound on T, also is considered as zero; 
this in no way restricts the generality of the results. The broadened range 
does ensure that when the stability conditions are met, greater confidence in 
the physical realization of stability follows. On the other hand, the number of 
cases which are stable is reduced; the modification required on T, is not too 
significant however, and is justified by the simplification provided. 

There has been an attempt to associate some of the parameters in the 


human-operator transfer function with physiological characteristics; the con- 


sensus is generally as follows: 


1) USAF Instrument Pilot School, Randolph AFB, Texas, Control Display 


Pilot Factors Program, Instrument Evaluation Project Report 63-1, Dec., 


1963. 
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Fig. 2-2: The Human Controller as a Monitor 





KR : human-controller gain; normally adjusted for overall 
system stability and low-frequency performance 
T : pure delay (characterizes human reaction time) 
Th Ty : equalization; these provide the human operator's 
adjustable compensation similar to lead or lag networks 


in "inanimate" systems 


Ty : neuromuscular lag; appears to be partially adjustable 
for the task 


This study does not enter the debate on the validity of these assignments. 


Zz. 6 Summary 


Normally, the human-operator transfer function is specified as a 
function of the controlled-element form and type of input power spectrum. 
Thus, any single description for the human operator really is valid only for 
a set of limited conditions. This study demonstrates the human-operator 
transfer function approach may be generalized by using a variety of published 
experimental data to form ''capability bounds" upon the transfer function para- 
meters. With such ranges defined, the human-operator capacity for adapta- 
tion can be displayed and used in conjunction with controlled-element system 
dynamics to determine certain ''man-machine'' closed-loop stability char- 
acteristics. This technique, to some extent, depends upon the assumption 
that the human operator uses as a basic criterion the maintenance of system 
stability; this assumption has been verified as a reasonable one after many 
experiments. A man-monitored system configuration also appears to be sus- 
ceptible to such stability analysis. The difficulty experienced by a human 
operator in attempting to maintain control of the system is implied by the 
"margin of human capacity'’ remaining; this capacity appears graphically as 


a part of the analysis technique. 








————— 
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CHAPTER THREE 
HUMAN OPERATOR ADAPTIVE RANGES AND STABILITY THEORY 


3.1 Introduction 


In this chapter the system transfer function mathematics are presented 
and the relations are manipulated so the adaptive properties of the human operator 
are isolated. The concept of a nominal transfer function is introduced to facilitate 
the separation of the closed-loop-system function denominator into a fixed and a 
varying part. Finally, a brief discussion of stability interpretation for this 
modified Nyquist display is given. 

Many of the concepts in this study may be found in different form through- 
out the engineering literature. Several such published papers led to considering 
the possibility of a study on the effects of variations in the parameters of a trans- 
fer function the FORM of which was known. OTS) Also, techniques similar 
in certain respects to those presented here are offered in Horowitz's text. i) Of 
considerable aid in contemplating the theoretical ramifications of a pure time - 


(7) 


delay was Choksy's work The utility of the method described lies in its inclu- 


(1) Blecher, F.H., les for Multiple Loop Transistor Feedback 


Design Princip 
Amplifiers, Bell Labs Memo. 57-2172-6, Nov. 8, 1957 also appears in Pro- 


ceedings of the National Electronics Conference, Vol. XIII, Oct. 7-9, 1959. 
(The stability criterion in this paper is due originally to M. L. Curtis). 
— Myers, B.R., A Useful Extension of the Nyquist Criterion to the Stability 


Analysis of Multiloop Feedback Amplifiers, Proc. 4th Midwest Symposium on 
Circuit Theory, Marquette University, Dec., 1959 





(3) 


Jones, P., Stability of Feedback Systems Using Dual Nyquist Diagram, Trans. 


IRE Prof. Group on Circuit Theory, Mar., 1954. 
(4) Satche, M., Discussion, J. Appl. Mecho. p. 4190.) Pec.) 744. 


(5) Stewart, R. M., A Simple Graphical Method for Constructing Families of 
Nyquist Diagrams, J. Aero Sciences, p. 767, Nove., 1951. 


(6) Horowitz, I. M.,'"'Synthesis of Feedback Systems,’ Academic Press, N.Y., 1963. 


me) Choksy, N.H., Time Lag Systems,''Progress in Control Engineering", Vol. l, 
Academic Press, N. Y., 1962. 
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sion of a range on the parameters of interest so a specific expression of a trans- 
fer function need not be selected. The case of parameter variation in a transfer 
function of given form is precisely the state of affairs in certain man-machine 


systems where the human operator provides an adaptive transmission link. 


5.2 System Transfer Function Manipulation 


The closed-loop system function for Fig. 2-2 is found to be 

vo. HG +G_G 

R  TFH SG FG IG, (3-1) 
It is noteworthy that in solving the stability problem for this system several other 
configurations are handled simultaneously. For example, removing G. in Fig. 
2-2 and placing it in feedback around G_ leads to the same characteristic equa- 
tion as in (3-1). The stability of (3-1) depends upon whether any of the zeros of 
its denominator lie in the right-half plane. For stability, there must be no zeros 
of 

lich bs <= LE pee (3-2) 
in the right-half wane 
Equation (3-2) can be written 


sane, 
? Oo 
(6+G,.G)0 + pa : m 


eats, 


or 


1+L 


1 


ul oa + L,) 


where L_ is the open-loop transfer function corresponding to the compensator- 
Je 
plant and Ly is the open-loop transfer function corresponding to the human con- 


troller-plant path where the plant is somewhat more complicated than oF and 


ste 
7 


That is, no roots of 1 foe arias = 0 in the right-half plane. 





ae = — 


abe) 


is represented by G 


Saat) eee 
i eG 
se 


Of course, 1+L_ presents no serious problem in stability analysis; the tech- 
niques for stability checks and performance siudies in this situation are well 
known since the structure is FIXED PARAMETER. 

The term containing the VARYING PARAMETERS is given by, 1+L 


h 
where HG . 
an ~ 1+GG 7 one 
c p 
that is, G 


; 


G = SES ae 
p tte ieel 


The symbol, Ho , represents that part of the L, expression which has 


h 
variable parameters and G. that part with parameters nominally fixed. ; 


Since; 


+ 
1+L, 


this can be written as H 
1+G'H = 
p o 


n O 
n 


ie eee (3-3) 
Oo p 





where Hy represents some nominal value of Ho which is obtained by picking 
n 
any set of values of the varying parameters in the given range of their variation. 


Now, the zeros of 











H 
eG us 2 (3-4) 
jo ) omen) >! 
n O 
n 
are the zeros of 
le Jel 
men é ean L (3 5) 
= + = 
Ho ear A n 





The examples in this study are for Gg=9 inorder to simplify the presentation, 
permit an intuitive grasp of the method, and allow comparisons with experimental 
results reported in the literature, since few empirical data have been reported 
for the monitoring configuration. No loss in generality occurs, the condition of 
an lL} with right-half plane poles, a situation representing the worst effect ofa 
non-zero Go: will be discussed in an example. 
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where 


be = (Gear (3-6) 
n Be foe 


is the nominal loop-transmission with FIXED PARAMETERS, Fig. 3-1. The 


stability problem is now one of shaping the function with non-varying parameters 


O O 
Le such that in the face of the parameter variation in iT the zeros of a Ls 
are all in the left-half plane; the criteria ensuring stability will be discussed 





below. 


3.3. The Nominal Function 


Of course, A should be chosen to facilitate the analysis. Equation 
n 
(2-1), representing the human operator is 


-ST 
Ke (1+sT, ) 


H = s = ]W (3-7) 
, J 
O (1+sT,)(1+sT)) 


where the final ranges considered in this study are 


O< K 
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As an example, the nominal expression for A is taken as 


n 
_ _-.1s 
HO =e (3-8) 
n 
where it is clear that (3-8) is (3-7) with 
K, =1,7 = 0.1, T) = T, = Ty = 9 


Thus, 
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(a) 
oa nee “| a <n an eo eat em eS rg a 
| VARYING | : | 
PARAMETER | ) FIXED PARAMETER STRUCTURE 
| = STRUCTURE : 
| : 





(b) 


Fig. 3-1 : Separation of Fixed and Varying-Parameter Elements 


(a) Original System, (b) System after introduction of nominal function, HA 
n 
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H 

O -.ls 
fines e 

Ho >| 5 re 
re) Ke (T, st) 

(l+sT A 1+sT, 
or 
H 
OT (Ts + 1)(T,s +1) 





A, RATS pe 


It is noted that the human-operator gain K, can be lumped with the fixed elements 


h 
Le of (3-5). This results in a design for a single overall open-loop gain, Ka 
from which the human-operator gain needed, K. or the plant gain needed, as, is 


h 
operator study the plant gain would be fixed at some nominal value since the in- 


determined upon specification of one or the other, i.e., K=K — In a human- 


vestigation is focused upon the human operator aspects of the analysis. In general, 


this approach is used throughout, thus from (3-4) 





-TSs 
H K_e (l+sT, ) 
h L 1 
1+G H = = 14GH esas a = 0 
Tee a 
Pp OD Ho p oO. (l+sTy Its I Ho 
n n 
(l+sT,) 


I 


li KhekeG (i WtstoiteT)| 7 0 (3-10) 
h p p OF ItsTy Il+sT, 


The study regards only fixed values of 7 and covers the entire range of T 
variation by repetition over a sequence of values uniformly distributed, or 
otherwise, across 0.1 <1 < 0.4. Since the plots are in the Nyquist plane the 
additional burden is negligible and consists only of a phase shift in L the 


-j WT : : 
magnitude of e a being unity. 








where the plant gain, ~ , has been removed from G. and H has been set 
O 


STG , a 
equal to e . Equation (3-10), put into the form (3-5), is written 








H 
O7 (1+s Ty)(1+s T,) 
TS 
FH i aa ee +K To, e- = © (3-11) 
ae 
Throughout the remainder of this work the expression for va is written as 
Oo 


H 
On (1+s T,)(1+ s T)) 


oa = (Faye) (3-12) 


Oo 





and the human-operator gain is lumped into Lo as in (3-11). It should be noted 
that the Ls term of (3-11) has a FIXED PARAMETER structure, Fig. 3-1; each 
analysis considers only fixed values of 7 across the range of interest, since the 
plots are in the complex plane the effect upon ao consists only of an additional 


phase shift. 


3. 4, Stability Criterion 


The standard Nyquist criterion for fixed parameter systems gives informa- 
tion about the characteristic equation of the closed-loop system. Thus, if the sys- 


tem function is 
6 KP(s) 


oO . 
ge 1+KFP(s) 


then system stability is determined by whether the characteristic equation 
heh ec) nO 


has roots in the right-half plane. If the Nyquist plot of the open-loop transfer 
function encloses the -1 point, the number of such encirclements and their sense 
gives the number of zeros minus poles that 1+KFP(s) has in the right-half plane. 
If the system is open-loop stable, i.e., no poles of KFP(s) in the right-half 
plane, the number of encirclements gives directly the number of zeros in the right- 
half plane. Thus, in this instance, any clockwise encirclement is sufficient to 


guarantee instability. If the system is open-loop unstable there are poles in the 


Ss: 


- 
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right-half plane, and these must be accounted for in interpreting whether or not 
zeros are in the right-half plane. Clockwise encirclements are regarded as posi- 


tive and counterclockwise encirclements as negative. For example, if 


K 
ir = 
a=) s(-1 ei 
the system is open-loop unstable (there is a right-half-plane pole). However, the 


Nyquist plot yields one clockwise encirclement, i.e., 


number of zeros - number of poles =Z-P = +] 


but P is known to be unity thereby yielding Z=2. In this case, the system is 
closed-loop unstable, the characteristic equation has two roots in the right-half 
plane. 

Can the Nyquist criterion be applied to systems containing time lags? 
Such a condition is always present with a human operator in the loop! The Nyquist 
criterion indicates the number of zeros of a function in the right half of the s 
plane by the number of revolutions the mapping of the function makes around a 
critical point of another complex plane for all values of s bounding the original 
right-half plane. The primary condition ensuring the validity of the Nyquist 
criterion is the analyticity of the function; since the exponential term representing 
the time lag effect is analytic in the region of interest the criterion is directly 
applicable without modification . ‘ 


The stability criterion used throughout this study follows from (3-11) 


(l+sT.)(1+sT)) _ 


= + , ~ = - 





r 1)(2)(3 
More detailed justification can be found Bie ermen 2)(3) 
2) Satche, M., Discussion, J. Appl. Mech. ,p. 419, Dec., 1949. 


(ay Dzung, L.H., The Stability Criterion, ''Automatic and Manual Control”, 
Tustin, A., Cranfield Conference papers, Butterworths (publishers), 1952. 


(3) Choksy, N.H., Time Lag Systems, ''Progress in Control Engineering", Vol. 
1, Academic Press, N.Y., 1962. 
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The stability of the system depends upon the number and type of encirclements 
the plot of (3-11) makes of the origin of its complex plane. ” Variation in Th» Ty 
and Ty causes ({) to take on the form of an area, however, and rctation of the 
vector" representing (3-11) as the values of s traverse the bounds of the right- 
half plane must be interpreted with care. All points within the region of variation 
which result in a rotation of this ''vector"' such that instability is established 
represent combinations of values for the triple (T, . T) Ty) which yield an un- 
stable system. If no points exist which provide a stable result this reflects a 
Situation where the human operator cannot stabilize (or control) the given system. 
The concept will be examined in detail later by building from simple examples. 
Thus, the region (3-12) is plotted for all possible values of the triple(T, , Ty) Ty); 
O 
n 
i 
re) 
the 'hull'') are universely applicable for the selected range of variation on the 





the boundaries of the regions {) = (jw) (for various fixed W each is called 


time constants Th Ty Tw This range was obtained from empirical results pub- 
lished in the literature and was established by selecting the largest and smallest 
values that Th Ty and Ty ever assumed, as indicated from the many different 
experiments conducted to determine them. The universal plots may be used in 
conjunction with any system of interest. The first generalized version of this 
technique seems to have been given by Tones eer regions of variation were not 
considered. 

Rather than plot the sum (3-11) and observe any rotations about the origin, 
the following technique is used. The negative of Le is plotted, i.e., -L as well 
as the region of variation 0; Fig. 3-2(a) displays a hypothetical situation for 


explanatory purposes. If a vector is drawn from the -L (j) locus to the 02 


to 


The origin here plays the role of the ''critical point". 


(1) Jones, P., Stability of Feedback Systems Using Dual Nyquist Diagram, Trans. 
IRE Prof. Group on Circuit Theory, Mar., 1954. 
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"locus'', each end of the vector being at the same value of Ww, that ''difference 
vector'' is actually the function (+ L. and the rotations it experiences give sta- 
bility information for that function. That is, stability is determined from the 
nature of the origin encirclements, if any, that occur in the ''difference vector" 
plot, Fig. 3-2(b). The vectors shown in Fig. 3-2(a) are for specific values of 
the varying parameters as well as for particular frequencies. Thus, the locus 
of (3-11) may be visualized without carrying out the summation. The replace- 
ment of (1 by unity results in an ordinary Nyquist plot with the head end of the 
vector fixed at +1+j0 andthe other tracing out the -L G4) locus. Thus, asa 
dividend the plot of -L, allows immediate observation of the plant-stability 
characteristics in those cases when G = 9; the critical point from which stability 
is determined in such instances is at +1+j0. Examples will be offered in the 
next chapter. 

It is not necessary, of course, to introduce the concept of the nominal 


human-operator transfer function H , since (3-3) 
n 


1+H G =0 (3-3) 
ey ie 


could have been plotted in the form of either 


an 
Hee neath ee 
p 
or 
r _ 1 
-G) = ae (3-14) 


and the "difference vector'' between the two loci interpreted similarly to the 
explanation above. In (3-13) the design engineer using the techniques of this study 
would be required to plot the inverse frequency locus and in either (3-13) or (3-14) 
the general contours to be developed which define the region of variation would be 
complicated in their calculation by a linear phase shift due to the delay time. The 
use of the nominal expression for the human operator removes the delay from the 
varying portion of the plot and places the burden for its consideration on the FIXED 


PARAMETER portion, i.e., the user of the study always plots this locus since rt 
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represents his particular system. His plot, however, need be modified by the 
addition of a linear phase shift to a single contour only instead of to a bounded 


region. 


3.5 Summ ary 


Under the assumption that the human operator can be described by a 
transfer function composed of five variable parameters, the system mathematics 
are manipulated so that three of these parameters (time constants) are displayed 
together. The remaining two parameters, gain and pure delay time, are lumped 
with the rest of the system; their effects upon the closed-loop stability character- 
istics are easily determined, since their influence upon the Nyquist contour of 
the plant is readily observed. It is shown that stability analysis of certain man- 
monitored system configurations also is possible. Stability estimates are rea- 
sonably accurate as long as the human-operator description used continues to 
represent the human controller adequately; this condition seems to hold for 


angular frequencies below 18 to 20 radians per second. 
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CHAPTER FOUR 
THE STABILITY CRITERION AND THE VARIABLE PARAMETER MAPPING 


4.1 Introduction 


One intent of this study is to provide a technique for viewing human con- 
troller-system response characteristics which is more intuitively satisfying than 
previous approaches. As a result of some interesting system transfer function 
manipulations, certain man-machine systems can be studied for stability on the 
familiar Nyquist plane. The time delay inherent in the human operator is easily 
handled and the effects of many of the human operator's physiological character- 
istics are plainly observed. The goal is to introduce a new approach to handling 
the quasi-linear human-operator transfer function without excluding the enriching 
properties of human adaptive range and pure time delay. 

In this chapter a complete analysis of a very simple single-varying-para- 
meter system is conducted to provide familiarity with how the theory discussed in 
Chapter Three actually is used. Following this presentation, the multidimensional 
variable parameter mapping is discussed in detail. The mapping function is that 
portion of the human-operator transfer function which previously was isolated as a 
separate term (the three variable time constants). Certain properties of the map- 
ping which relate to necessary and sufficient conditions on system stability are 
considered. Finally, plots of the mapping are presented (also on transparent over- 
lays); it is these plots which are used in Chapter Five to predict human control 


capabilities. 


4.2 Analysis of a Simple System 


To facilitate an understanding of the modified Nyquist diagram used in this 
study consider the system depicted in Fig. 4-1. The root loci shown in Figs. 4-2 
and 4-3 are representative of this system, i.e., for the same open-loop gain K=K, 
it is possible to have values for the variable time constant x such that a stable 
closed loop occurs, Fig. 4-2; or an unstable closed loop prevails, Fig. 4-3. The 


encircled crosses depict closed-loop poles while the others represent open-loop 
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K : 


(xs +l )(astl)(bs +1) 





Fig. 4-1 : A System With One Variable Parameter 
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mae. 4-2 : Stable Closed Loop Fig.4-3 : Unstable Closed Loop 


K(STABILITY LIMIT) = 90 
K (STABILITY LIMIT) = 60 





meet -4 ; Effect of Variable Parameter 
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poles. Such a system, with only one varying parameter, can be studied by way of 
example to gain an understanding of the dual Nyquist plotting technique. Generali- 
zation to a three varying parameter mapping follows easily. 


A specific example of Fig. 4-1 where b=1 and a= = follows. The closed- 


2 
loop function has a denominator of the form 
K 2K /x 
D(s) = 1+ fh _____. — ] + I 
(sxt1)(> + 1)(s +1) (st—)(s+1)(s+t 2) 
If ae then 
a e 
8K 
= Vv = K _ 
25) Ss re PERE E AGS (ETS Inca © eo 
Therefore, 


Di= eae ee lee i 


it 
>) 


is the characteristic equation. From the Routh array, the stability limit for the 


root locus gain is found 


5° l 14 
ae 7 (8+ K) 
a. | GR SE 
7 
0 


fopbe K=90, .. K = ae = 11.25; the point at which the closed-loop poles cross 


the jW axis is 


Tie = -98 
ee = -14 
S = Pty s.74 


I 
This situation is shown in Fig. 4-4 by the solid line loci. Now suppose, x =3 
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instead of 2 then 
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UI 
© 


is the characteristic equation. Again, using the Routh array, the stability limit 


is found 


5° l ll 
aa b (6+K) 
a! BiG Kine 
elloae ath 
0 
S 


60 _ 


tobe K=60, .’. K are 10; the crossing of the jw axis by the closed-loop 


poles occurs at 


6s = -66 
ae = -ll 
SH =4 =) G32 


The loci for this system configuration are shown dotted in Fig. 4-4 as is the pole 
at - 2 -3. 
x 

lim 1s allowed to vary over the arbitrarily selected range 0< x= > the 
system will be characterized by both stable and unstable response properties 
depending upon the value x assumes, provided the gainis "properly" fixed. Thus, 
for K=90 the system is at its stability limit for = = 4, andis surely unstable 
when = = 3 for the same gain, since the gain driving the closed-loop poles onto 


the imaginary axis for == 3 is only K=60. Thus, at K=90, there is system 
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stability for => 4 and instability for = < 4, (S€e Pip 44-5. 
The plot of these situations using the dual Nyquist approach is now con- 


sidered. Again, letting 


eS 


D(s) = 1 + 
(sxPINet ais +1) 


and with the previous knowledge that 


Kee 2) bc eanaeac = 2 
Vv 4 
corresponds to the stability limit, while 
1 
K =11.25 and x== 
Vv 3 


corresponds to instability, the denominator of the closed loop function, D(s) is 
rewritten. 

D(s) = 1+ V(s)R(s) 
The V(s) portion will contain the varying parameters andthe R(s) portion the 


fixed ("'rigid'') parameters. 


Thus, 
i 
es eo) 
and 
R = A252 (4-2) 
(st 1)(5 +1) 
so 
Me 
D = 1+RV v (4-3) 


The nominal function, Me , is taken as unity, i.e., x=0, in the expression for V. 


Then the roots of 


are those of 


(4-4) 
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or 


(exe) eee, (4-5) 


(s+ IS + 1) 


In equation (4-5), the expression (4-2) plays the role of Le in (3-5). Thus, 


aes. 





Oks car cael (4-6) 
(s+ (5 + 1) 
and similarly, 
Vv 
= ay = (sx+1) (4-7) 
fe) 
plays the role of — = 9 
Oo 


In analogy to the method described in Chapter Three, the expressions for, 
-L representing the fixed parameter contour and, ata the region of criticalzn, 
(here a straight line) representing the varying parameter structure,are plotted. It 
is noted that this varying parameter structure is characterized by the mapping ta 
over the range of x and for specific frequencies W™. The remainder of the exam- 
ple clarifies this notion. 

The expression (4-7) plots as a straight line for any given wW with its bounds 


defined by the end points on the range of x. For example, 


let 
s=jW=j i.e. W=] 
then 
V 
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V 
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The results are plotted in Fig. 4-6 and are represented by the vertical line running 
from 1+j0 to 1t+j® as W+e, Actually shown on Fig. 4-6 are the situations 
reflecting Y™=0.5, W=1] and W=2, Thus, when W=1, the ae plot may assume 
any point along the vertical line segment from 1+j0 to 1+j5 depending upon the 
particular value of x. This segment is defined in Fig. 4-6 by the vertical distance 
dimensioned W=1, and the others similarly. 

A plot for the FIXED PARAMETER structure, -L is obtained from the 


negative of (4-6). It is easy to see that 


-L_(j0) = -11.25 = 11, 25 hi 80 Sl, 2 |-180 


and that the 
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lim L-L (jo) = 0 lo 


WwW © 


thus providing the end points of the contour. Other points may be found in the 
normal, tedious fashion or by transference from a Bode plot. This contour is 
also shown in Fig. 4-6. 

It is difficult to see at this point that the ''difference vector'' drawn from, 
ae to the > contour (head end at the latter), and connecting identical values of 
W oneach contour, does itself undergo a clockwise rotation if sie 5 . Without 
some experience viewing the plot, initial interpretation may be formidable; for 
this reason, consider x assuming only two specific values in its range. When 

3 


. 1 ; 
x == , the system is stable, and when x=-, the system is unstable; these two 


conditions are now studied. Thus, for 


11.25 
EF .(s) = —-——— (4-8) 
l S Ss 
and 
F(s) = 2 (4-9) 


there is stability with (4-8) and instability with (4-9). This can be checked quickly by 
(1) 


counting sign changes in a Routh test. 


It should be recalled that 


EF (s) = tees (4-10) 


l Ss S 
<a me +1)\(2+ 
represents a system precisely at its stability limit. 


Separation of the varying (in this instance, "fixed", but arbitrary) and 


fixed parts of the characteristic equation results in (4-5), repeated below. 


7) Gardner, M.F. and J.L. Barnes, ''Transients in Linear Systems", Vol. I, 


Wiley, New York, 1942. 





fee, 


(exch) = ec (4-5) 
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where —— or =. 
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Thus, for x = { 
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The -L_ plot is the same as inthe Fig. 4-6 plot and is reproduced in Fig. 4-7; the 
(=+ 1) plot is again a straight line (this time for varying “ and fixed x) upon 
which are marked, to the right of the line, the corresponding frequency points. 
It is noted in Fig. 4-7, that the difference vector (shown for W= 1) never under- 
goes a complete encirclement (360° or more). At “ = 0, the difference vector 
points toward, 0” (due right), while at slightly less than W= 2 (first intersection 
of -L_ and =), the difference vector points toward -90° (straight down). Some- 
where in the neighborhood of W = 3, this vector experiences its largest clockwise 
excursion (approximately -120°), after which it begins a counterclockwise return 
rotation until at “ = 4 it has almost returned to 0°. As W-®, it undergoes a 
counter-clockwise rotation approaching 90° in the limit. Clearly, a complete 
encirclement does not occur. 

For x = : 
The -Lo plot is, of course, unchanged. The (5+ 1) plot superposes over the line 
describing (= +1) but with different frequency correspondences. These are shown 
by the markings which appear to the left of the vertical line in Fig. 4-7. It can be 
seen from Fig. 4-7 that as “ runs from “ =3 to W=4 the difference vector 
between the two contours undergoes an encirclement. It is clear that whether or 
not encirclement occurs depends upon the relative values of frequency on the -L_ 
and 5 + 1) contours. 

A magnification of the situation is contained in the enlarged sketches of 
Fig. 4-8. In Fig. 4-8(a) is shown the hypothetical situation where W | (at inter- 

v n 

section with ~~) =Io.7 and). ve (at intersection with -L) = 4.7. The effects of 


“Vv 
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encircling the right half of the s-plane upon the mapping of the difference vector 
between -L_ and (sx+1) is shown in Fig. 4-9(a) and (b). Ina qualitative plot 

of the difference-vector's angular rotation, Fig. 4-10, there is no encirclement 
of the origin and the system therefore has no right-half-plane roots in its char- 

acteristic equation; since the open-loop system is stable this conclusion is per- 
missible. Thus, the system is closed-loop stable. 


A different state of affairs is depicted in Fig.4-1]](b) where Ww L (inter- 
n 


section) = 3.7 and Ws (intersection) = 2.8. The effect upon the mapping of the 
n 


V 
difference vector is such that, Fig. 4-1l1(a),(b) and the qualitative plot of Fig. 4-12, 
two encirclements occur indicating a system with two closed-loop poles in the 
right-half plane and unstable characteristics. 

In this simple case, it can be inferred that if the value of the frequency on 
the “Lo contour is greater than that on the locus of the varying structure, at the 
eritical intersection point, then instability occurs. In short, if the “Lo locus 
"encloses'' the varying locus there is an unstable condition; "encloses" is inter- 
preted to mean: if the frequency at the intersection of the two loci is W, on -L 
and if in approaching the point -LG 2) as ive -Lo contour is traversed in the 
direction of increasing frequency, the point —y~(j¥,) is passed to the right of the 
-L_ contour (called enclosure), then the system is unstable. Again, there is 
inherent here the asumption of open-loop stability. 

Before this simple example is left, the precise values of the intrinsic para- 
meter W (at the points of intersection) are calculated for each contour. This is 
done for completeness of the example. It is not required for the stability test 
under ordinary circumstances since, with practice, encirclement can be envisioned 
from the relative locations of discrete frequency points on the dual Nyquist plot. 


For the stable configuration: 
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Fig.4-9 : Hypothetical Crossing Situations (in the large) 
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Fig. 4-11: Hypothetical Crossing Situation (in the large) 
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IMAG. 


Fig. 4-12: Excursion of the Difference Vector 
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Vv 
Re (jw) = Re [-L (5%) =] 


where ''Re"' stands for ''the real part of". 


Thus, from 





it is found that Ww) = 1,87 and W = 3.75, where the subscripts indicate the fre- 


quencies of first and second intersection respectively. 

It is now clear (since in this simple case the locus of the varying portion is 
the same straight line but with ee frequency markings) that, when x == 
-L (33. 75) meoes not ''enclose" Se 75); when x = 5 ; -L (j3. 75) does 


"enclose" + (3.75). Of course, this statement can be checked analytically. Thus, 


(4-11) 





where "Im'' signifies "the imaginary part of". For w= 3.75, (4-11) evaluates to 


0.936. But 


Im [-L(jw)} = Im [ve (ju) 

at the point of intersection. Therefore, for x a 
Ve r 
Im (J) a se 0. 936 
or V 
w = 4,68, on ay (j¥) at intersection 
V 
and, clearly, 7 (33. 75) is to the left of the contour -L Ww), indicating a stable 
system. 
For the case when x = 5 , similarly, 


V 
ny; = wv = 
Im 7 (5) = 3 0.936 
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or V 
oan * A 
i = 2258 1>.0n a7 (3) at intersection 


V 
and, clearly, =~ (33. 75) falls to the right of the contour -L (jw), indicating an 


unstable system. 

A thumb rule for interpreting the stability criterion in this simple situation 
has now been established. Inthe original example, Fig. 4-7, an attempt at inter- 
pretation can be made using the additional insight just gained. It is known that the 
intersection of interest occurs at the point -L (33. 75), and instability obtains for 
any value of x such that 


Ww (intersection) < vo. (intersection). 


Vv 
—2 n 
Vv 


In the particular case where the locus of the varying portion is a straight line 


there is the condition that 


faecal) = ek x = 0.936 


But for stability, 


therefore, since 





for stability 





Thus, for stability, x<.25 which reflects precisely the conditions seen in the 


root locus plots of Fig. 4-4 and 4-5. 


4,3 Multidimensional Variable Parameter Mapping 


The system examined in Section 4. 2 had but a single varying parameter. 
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It is for this reason that the locus of the critical point governing stability moved 
along a single curve, in that example, a straight line. When more than one para- 
meter in a system varies (multidimensional parameter variation), the locus of 
critical points assumes the form of an area. In a well-behaved situation, for 
example, a mapping of the boundaries of varying @ and B, Fig. 4-13(a), appears 
as in Fig. 4-13(b). In this case, every element XY which scans the constrained 
“range of @8 space is inside the mapping's periphery 0'1'2'3', called the "hull". 
In such a case the mapping is ''complete" and with a suitable generalization of 
Section 4. 2 the stability criterion is necessary and sufficient. 

In the general case, however, end point mapping of the region of variation 
produces results as appear in Figs. 4-14(a) and (b). Here, for a given distance 
along the Q-axis the element a Sh Ga maps into a line contained inside the hull; 
but, further along it is clear that AYL has points which are not contained in the 
hull, e.g., XY, Fig. 4-14(b). Thus, although the end points of the —qp line 
segments are the images of the end points of the rectangular regionin 8 space 
the mapping may not be ''complete''. It is clear, however, that all elements xe aa 
are inside the triangle 0'1'3‘; this "polygon" is called the "convex hull". 
Obviously, the convex hull contains points which are not a result of mapping the 
a8 variation. It is for this reason that exclusive use of the convex hull provides 
only a sufficient condition for stability, i.e., for no "encirclement" a system is 
stable but, since an ''encirclement''might be of a point or points not part of the 
mapping, an "encirclement" does not necessarily indicate ''instability'’. This 
becomes clear in the examples. It should be noted, however, that a complete 
mapping can always be obtained. One need only plot a few points such as x nae 
of Fig. 4-14(b) and the true hull can be sketched quickly. It is shown later that a 
closed form expression (a parabola) can be obtained for the envelope of the lines 
X'Y’ in Fig. 4-14(b); this will have direct application in the human operator 
Situation. 

Frequently, the requirement for mapping additional points is precluded by 
the relative location of the fixed and varying contours. For example, assuming the 
situation shown in Fig. 4-15 holds in a similar way for all w of pertinence in 


addition to » =), it is obvious that whether 0'1'2°3" or 0 1 22 3 Misused 
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IMAG. 


MAPPING: w=w 





Fig. +-13.: Boundary Mapping Resulting in Complete Hull 


ee. MAPPING: w=w, 





(a) (b) 


Fig. 4+-l+ : Boundary Mapping with Hull which is Not Complete 


(also shows Convex Hull) 
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IMAG. 
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REAL 


Fig. +-15 : Hypothetical Variable Parameter Situation 
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makes little difference. The shaping of ale is carried out with respect toa 
region common to both. This property does not hold dogmatically when the human 
operator is considered, since being an adaptive mechanism, he will strive to 
assume characteristics in that part of his ''varying region" which afford stability. 
It is through this understanding, however, that insight may be gained concerning 
the amount of effort the human operator has to expend to maintain control. At 
times it is possible to demonstrate that a mapping of the end points in 4B space 
yields almost the true hull and is negligibly different from the complete mapping, 
so that practical use can be made of the end point mapping. This is the situation 


with the human operator and the mapping is called ''almost-complete". 


4.4. A Map of the Human-Operator Adaptive Range 


With the human-operator's gain-varying capability and delay time lumped 
into the fixed parameter structure, the varying portion of the characteristic equa- 
tion is given by (3-12), repeated here 
LH Clits. -)( 1 36} 
ee ae 2 : (3-12) 

H (1+ Sy 


O 





In order to study the human operator pure delay effect, 7, over its entire range of 
variation, several plots for different values of T must be made. Usually, due to 
the limited range of 1 and the prospect of interpolating middle range effects from 
the end bounds on T, one or two plots will suffice. The plots (of -L_. (3-11)) are 
all identical except for the phase lag introduced by the particular value of T used. 

The region of variation for Ty Ty and Thy is shown in Fig. 4-16, i.e., 
equation (3-12) is used to map this three-dimenstional volume into the complex 
plane plot of (jw). The region covered by ‘) as the triple (T) . Ty Thy) assumes 
all values on the bounds of the rectangular parallelipiped is shown in Fig. 4-17. 
(Later the completeness of this ''end point mapping'' is considered and it is shown 
that the end point hull is almost complete as described in Section 4. 3. ) For motion 
along the boundary AB of Fig. 4-16 there is obtained 
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(T,s + 1)(0. 75+ 1) 


Se rn 55-51, (4-12) 


but for a fixed w this is the equation of a straight line in T |. Fox example, let 


w=1, then 
(Tj + 1)(0. 73+ 1) 
(5. 33+ I) 
and after some manipulation 


T (4. 6 + j4. 71) + (4, 71-34. 6) 





which is recognized as the general vector equation of a straight line 
isp Ss Fass 


Since T; has a range of interest, its end points, T)=0 and T= 25, define a line 
segment in the space of the mapping 9; this line segment has end points given by 
the images of T=0 and TyF 25, Were 


at T;=0 


a= Ati ° = 0. 162-j0.159 


while at T= 25 


_ 25(4,6 +54. 71) + (4. 71-54. 6) 
Y= 29 


Q = 4,13+ 53.91 


This line is plotted as A’B’ in Fig. 4-17. It is obvious that lines GF,OE and 
DC map ina similar manner into straight lines in the image (or mapping) space. 
Also, variations in Th for fixed TD and Ty map into straight lines; thus, 
OD,GA, FB and EC can be mapped easily. The results of this mapping yield the 
end points in the image space for w=1 as shown in Table 4-1. | 


These lines are mapped accordingly, Fig. 4-17; their images are shown by 





ee rr— Oo <3, a ee 


72 








TABLE 4-1 
Line Mapped TL T Ty End Points of Image Line 
pe 0 Onn 0. 162-j0. 159 
AB 
53 25 On 4.13+j3.91 
Sas 0 0 0. 0345-j0. 183 
GF 
See 25 0 4.61440. 679 
0 0 0 l 
OE 
0 2 0 1+ 525 
0 0 0.7 Oe 
DC 
0 25 0a Woes 258 0 
0 0 0 l 
OD 
0 0 0.7 Ler Ge 7 
558 0 0 0. 0345-j0. 183 
GA 
Spas) 0 OG 0. 1625-j0. 1586 
oe 25 0 4,61+ 0. 679 
FB 
5.3 25a. 7 4,12+ 53.91 
0 25 0 1+j25 
EC 


0 25 0.7 ~16.5+525.7 
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the same letter designations but primed. Thus, faces ODCEO and GABFG of 
Fig. 4-16 map into the triangles and line segment appendages of O'D' C Eo: 
and G’A’‘B’'F’'G’ of Fig. 4-17, respectively. 

The mapping of lines AD, GO, FE and BC into 9 space can now be con- 
fs 


sidered. The most general case, i.e., T #0, is represented by line BC. 


| eS) 


Again, for w=l, 
Ce (jZ5 +1) 0.87 1) 
(T, 5+ 1) 
= ~16.5 4525.7 (4-13) 
1+ jt, 


Equation (4-13) is recognized as the equation of a circle; a proof is given in 
Appendix A. i Since three points determine acircle, the arc corresponding to 


(4-13) can be found from the two end points on the range of T, and the property 


L 
that as tT © the circle passes through the origin, i.e., 


Q(jw)?0 as Tg 


Thus, on BC the three points are: 











H 
pe AW (jl) = 0 
Oo 
H 
O 
T, =0 (jl) = -16.5+525.7 
O 
H 
On 
1 See (jl) = 4.134 53.9 
O 


1) Churchill, R. V.,'"Complex Variable and Applications", 2nd Ed., McGraw- 


Hill, New York, 1960. 





Alternatively, the equation of the circle is given by (A-8) where in this case, 
m--10.5, b = 25.7 yielding 


ft 6o 25)° Piya es cue (4-14) 


Of course, only the arc of the circle corresponding to the range 0< T, <5.3 is 


L 
of interest. In a similar manner, the mappings of lines GO, AD and FE are 


obtained; the information, for w=1, is assembledin Table 4-2. 








TABLE 4-2 
. End Points of Equation of 
Line Mapped T, T, Ty Image Arc Corresponding Circle 
GeeeeS Oey 116. 5 aang AG alewee 2, asc 
BC 
=(15, 28)° 
5.3 25 0.7 4,.13453.9 
Bese 14 525 (u-5)°+ (v- 12, 5)° 
FE 
=(12. 53)° 
5.3 25 O 4.61+50.679 
0 0 0 ] (u -s)* tv" 
GO TAZ 
=(5) 
ae. 0 0 0. 0345-j0. 183 
pme0. On7 144087 (u-5)° +(v-0. Pay 
AD may eiles) 


5.3 0 0.7 0. 1623-j0. 1586 





It should be observed that, interestingly enough, the equations of the 
circles which represent variation in Ty are not themselves functions of TL 


although the end points of the arc segments reflecting the range of Ty surely 








75 


are. Hence, as long as the ranges on Ty and Ty remain unchanged so do the 
equations of the circles. This property proves to be extremely useful in this 
study later on, and will provide an opportunity for obtaining insight and intuitive 
understanding of the effect of human operator adaptation upon any given system. 
The end result of Tables 4-1 and 4-2 is shown plotted in Fig. 4-18 and 
represents the mapping (3-12) plotted for the single frequency w=1. A sufficient 
number of such plots must be made for selected frequencies to permit determina- 
tion of whether the locus of (3-11), the difference vector between -L_ and 
Ea a 
= 0), undergoes encirclements of its origin or not. From this, stability 


H 
fe) 


information may be inferred. Since the human-operator transfer function is valid 
only over a limited frequency range, the plots are not infinite in number; addition- 
ally, their tendencies are obvious and a few judiciously chosen discrete frequencies 
are seen to suffice for the purposes of this study. Several points are worth noting 
at this time: 

1. The human operator's pole variations map into straight line segments 
while his zero variations map into arc segments. 

2. The formation of a largest polygon (with some arc boundaries) around 
the bounds of the region in Fig. 4-18, i.e., formation of the convex hull, defines 
an area of possible criticality re stability. 

3. It must be determined if the end-point mapping results in a complete 
mapping, i.e., every point of the mapping {2 must be a point of the end-point 
mapping if the mapping is to be complete 


"1 


4, If the hull is not complete, rather than use a convex hull, can the “end 


point" hull be shown to be almost complete? 

5. The situation has changed progressively through the usual Nyquist con- 
dition where the critical point is -1+j0, to a one parameter variation where the 
locus of critical points is a curve, to a final condition where the critical points 
are defined by an area °). 

It is shown that Point 4. pertains, but that the hull of the complete mapping 
can always be found and may be described analytically; however, the end point 


mapping of Fig. 4-18 differs so little from the complete hull (indeed, very little 
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from the convex hull also) that for all practical purposes it may be used. 

The interpretation of Point 5. is: if there is some area of 2 unenclosed 
by the contour -L_ then presumably the human operator can assume such para- 
meter values and stabilize the system. The parameter values necessary to ac- 
complish this give insight into the amount of effort the human operator must 
expend to retain control. If the entire ( region is enclosed, the human operator 
cannot control under any circumstances (for the capability ranges assumed and 
plotted). If no portion of the region is enclosed, there is stability no matter 
what the actions of the human operator (within the sense of the transfer function 
description). Later, examples and discussion clarify these notions. 

For the study to proceed, the regions of varying parameters, (, must be 
plotted for selected frequencies. Generally speaking, the human-operator trans- 
fer function has validity only over a range of frequencies, w<2Q0. It can be seen 
that beyond an w = 3, the plots offer little additional information and their trend 
is obvious. The region {© is plotted therefore, for w=0.5, 1.0, 2.0, 3.0 and 
also for an w = 10; in order to do this it is necessary to derive the mapping 
relations for the various bounds of Fig. 4-16 as functions of w. Two such 
derivations are given below; the first for a human-operator pole and the next for 
a human-operator zero, the remainder are obtained in a similar way. 

The human-operator pole variation, Ty , along AB of Fig. 4-16, is found 
from 

(Tp jw +1)(0. 7j w+ 1) 


After some manipulation, 


(T jut 1143. 710°- 4, 65) 





Qos 5 (4-15) 
1+28u 
Thus, at T= 9, 
143, 71w* - 4, 6j: 
QO = 1+3, fiw - 4.05 (4-16) 


ees io 
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mad, at I,=25, 


I 


5 = L#118. 71+ jw(92. 802+ 20, 4) 
~ 2 (4-17) 
1+ 28w 


The mapping of AB is a straight line between the points obtained from (4-16) and 


(4-17). For variation in the human-operator zero, T,, along BC in Fig. 4-16, 


| re 
q = (2djwt 1)(0. 7jw+ 1) 
(T, jwt J) 


_ (-17. 50%+ 25 jwt 0. 7ju +1) 
1+T) jw 


If the complex number Q(jw), is represented by utjv, then 


Zz ; : 
1+T, jw e 1-17. 5Sw (aele Ue ey 


We jv u-jv 


a tallte SH es Cyne ete Il tie Bh ove 25. (UN, 


u tv 
The real parts may be set equal to each other (the imaginary parts permit recovery 


of T, for given w, Appendix A): 


ie u-17. Sane on 7 wv 


Woty 


If the square is completed, one obtains 


Z 


ara 2 Re 
eae ) |e -25 2 | (ee (25%) aie 


Relation (4-18) is the equation of a circle, an arc segment of which, BC maps into, 








as a function of frequency. 
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At To : 
= 1-17.51 +25. 7j (4-19) 
and at T=. Ye ; ; 
ns cy ions sea. -8i + 20. 
j ee 20. 4) (4-20) 
1+ 28: 


In a similar manner, corresponding relations are obtained for the remaining 
bounds of Fig. 4-16; these are displayed in Table 4-3. 

It is now possible to plot the mappings of the parameter variations for 
Seeeiewlar frequencies; these plots are shown for « =0.5, 1.0, 2.0, 3.0 and 10.0 
in Figs. +-19, 4-20. 4-21, 4-22 and 4-23, respectively. The question of the com- 
pleteness of these end point mappings may now be considered. As is clear from 
Fig. 4-24 the end point mapping is not significantly different from the convex hull. 
Mhis is particularly true in the light of system plots of “Le contours for typical 
physical plants, as can be seen in later examples. The difference between the end 
point hull and the true hull (complete mapping), however, is even less than between 
the end point hull and the convex hull. In the complete mapping, a small parabolic 
arc extends between segment O’D' and D’'C’, as well as between G‘A’ and 
6B . For the purposes of this study, this situation may be ignored. Details and 
feereoot are contained in Appendix B. 

One further point of some significance relating to the mapping should be 
noted. As stated above, face ODCEO of Fig. +-16 maps into the triangle plus 
line segment © DC E’O , shown in Fig. 4-17; similarly, face GABP Gun 


into G’A’B’F’G’. The former is associated with a T, =0, while the latter 


L 
corresponds toa i ee 5.3; this value represents the upper bound on the range for 


TL used inthis study. It is interesting to observe, and useful later in the physical 


/ / / 
interpretations, that as ee varies from 0 to 5.3, face O'D'C'E'O’ "rotates" 


. 2 ¢ 2 / / / u / ° . 
maporcoincidence with face G AB F G3; the “rotation” occurs by virtue of con- 


strained motion of the significant apex points of the "triangle” along the circular 


mee sepinents of Fig. 41-18. 


It may be preferable to reflect upon the effects of a two parameter mapping 
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TABLE 4-3 
End Point Value General Equations 
Line Mapped of Varying Parameter of the curve 
0 1+3. 71w*-4, 6ju 
1 ee 
AB a 
! Z 2 
25 1+118. 71lw +jw(92. 8w “+20. 4) 
1 + 28w° 
0 1-5. an 
1 + 28w 
GF T, = 2 
25 1+132.5w +19. 7 jw 
1+ 28 ihe 
0 l 
OE Ty = 
25 1+25jw 
0 Maes (0 
DC Ty = 
2 
25 1-17. 5w +25. 7iw 
5 1 +132, 5u“+ 19. 7a 
ae 28° 
FB . Ty = 
nae 1+118. 71u“+ jw(92. 80+ 20. 4) 
1 + BAGG 
a dn i 
OD Ty = 


0.7 aes (arts 
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TABLE 4-3 (cont'd. ) 





n Oint Value eneral Equations 
Line Mapped of Varying Parameter of the Curve 
0 1-5. 3jw 
GA 7 7 Leen 
Ne 2 
et So lite = 4 
0.7 SEEN Aa 
1+ 28w 
0 1725 jit 
EG Ty = : 
0.7 1-17. 5w te Ot a) 
0 ei me ob, Gaal 
oe 2 
Cu-(. 5-8. 75w°)] +€v-12. 85w 
BC Ty = 


2 
Sige gc eva, isi 


a 1 +118, 7lu“+ jw(92. 8w~+ 20. 4) 
1+ 28w 


cot e-ap + = BD 


5 3 —— Ne 
1+28w 
0 ears (tl 
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| Mis ok 21w*-4, 6jw 
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Fig. 4-19: Variable Parameter Mapping 
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Fig. 4-21 : Variable Parameter Mapping 
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before continuing on to the practical results, derived from the three parameter 
mapping, which are described in the next chapter. To this end, and also to offer 
further stability interpretation practice, consideration of Appendix © may be 


appropriate at this time. 


4,5 Summ ary 


After the analysis of a simple system containing only one variable para- 
meter, the situation of multi-parameter variation is considered. With the famili- 
arity gained from the single parameter example the generalization to a multiple 
parameter problem follows more easily. In particular, the three-varying para- 
meter problem (one zero and two poles) is considered in detail, since it relates 
directly to the three time constants in the human-operator transfer function that 
are within the adaptive capacity of the humancontroller. It is shown that these 
parameters, forming a type of variable structure, can be used as a mapping func- 
tion in the complex plane to display the boundaries of human-adaptive capacity. 
These boundaries contain a collection of ''critical points", akin to the -1+j0 
point on the usual Nyquist diagram; when the plant Nyquist contour is plotted in 
conjunction with these bounded areas (at particular frequencies) interpretation of 
system stability characteristics is possible (Chapter Five). Several plots of such 


boundaries are presented for important frequencies in the human range. 
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CHAPTE Kemi» 
APPLICATIONS OF THE, Tower vy 


5.1 Introduction 


The groundwork which is laid in Chapters Two, Three, and Four now 
makes it possible to consider actual physical systems which are operated by 
human controllers. If the empirical results can be predicted (or implied) by the 
approach outlined in the previous chapters, then this technique for viewing the 
response characteristics of systems with human operators in the loop is justified, 
and difficulties with human controllability of such systems can be anticipated via 
ime theory. 

Of the several experimental studies performed with human controllers, 
two are particularly well suited to the purposes of this investigation. The Hall 


(1) 


study is notable for its comprehensive nature and for the use of a non-trivial 
controlled element. The predictions obtained from the theory presented in the 
previous chapters can be checked for their physical validity using the results of 
the Hall experiment. The Smith ee is derived from human attempts to con- 
trol an unstable system; with these empirical results, the effect of open-loop 
instability upon the theory also may be compared with actual physical accomplish- 


ment. 


5.2 The Hall and Smith Studies 


In 1956 and 1957 a simulator study of the human pilot was conducted by 
Ian Hall at Princeton University. (1) The intent was to investigate the effects upon 
the pilot of various controlled-element-longitudinal dynamics. A very excellent 


analysis by McRuer etal, uses the data of that report as a basis for developing 


(7) Hall, I. A. M., Effects of Controlled Element on the Human Pilot, Princeton 
Univ., WADC Tech. Rept. No. 57-509, Aug., 1958. 


‘2. Smith, R.H., An Experimental Determination of the Limits of Pilot Con- 
trollability for Unstable, Oscillatory Second Order Systems, Thesis for degree 
of Aeronautical Engineer, Univ. of Cincinnati, 1961 (also appeared as a paper, 
On the Limits of Manual Control, IEEE Trans. on Hum. Factors in Electronics, 


Sept., 1963). 





a ''flying qualities" theory. 2 


Briefly, the McRuer approach follows the "adjust- 
ment rules" procedure outlined in Section 2.1 of this study. Another experimental 
study conducted at the University of Cincinnati by Ralph Smith probed for the 


(2) The results of 


human limits in controlling certain open-loop-unstable systems. 

both these studies are used to significant advantage here; indeed, the empirical 
2 

evidence of Smith! and seal) furnished the experimental justification for the 


results of this study. 


5.3 The Hall Study 


The experiments were conducted in a mockup cockpit simulating a Navion 
aircraft. The subjects were two Navy pilots, one with considerable jet time 
(Pilot I) and the other with a preponderance of propeller aircraft experience 
(Pilot i. The seemingly subjective pilot rating scale of the Hall study used the 
following definitions (quoted): 

"GOOD (G): aircraft response very suitable for the given tracking task, 
ACCEPTABLE (A): possessed of satisfactory flying qualities, but less 
easy to control than a ''good" configuration; aircraft response either 

faster or slower than that which the pilot visualizes as most desirable. 
POOR (P): possessing characteristics which require considerable con- 
centration to accomplish the given tracking task; configurations 

are either quite touchy, unresponsive, or poorly damped. 


UNACCEPTABLE (U): quite unsuitable. 


ate 


The author discussed the study with Pilot Il in November, 1964. 


i) McRuer, D.T., I. L. Ashkenas, C. L. Guerre, ASystems Analysis View of 
Longitudinal Flying Qualities, WADD Tech. Rept. 60-43, Wright Air Dvlpmt. 
Div., Jan., 1960. 

(2) Smith, R.H., An Experimental Determination of the Limits of Pilot Con- 
trollability for Unstable, Oscillatory Second Order Systems, Thesis for degree 
of Aeronautical Engineer, Univ. of Cincinnati, 1961 (also appeared as a paper, 
On the Limits of Manual Control, IEEE Trans. on Hum. Factors in Electronics, 
mept., 1963.) 

(3) Hall, I.A.M., Effects of Controlled Element on the Human Pilot, Princeton 
Univ., WADC Tech. Rept. No. 57-509, Aug., 1958. 





The subjects differed only between P and U. Pilot II tended to classify aircraft 
which were unsuitable for the given tracking task, but could conceivably be flown, 
as P, reserving U for configurations which he considered as unsafe to fly under 
any circumstances, while Pilot I tended to classify all configurations with which 
tracking was very difficult as U.'' Due to the use of plus, ‘'somewhat more than", 
and minus, "somewhat less than'' modifiers on the G, A and P ratings the resul- 
tant number of classifications was ten, i.e., A+ 7-4, G72, Gt+-1l1, etc. Those 
familiar with the military aircraft rating scheme are aware that the standard 
Cooper scale also has ten gradations; they are not exactly the same. vo It may 
appear that the entire rating system is seriously suspect by virtue of its apparent 
subjectivism. The consistency of pilot opinion, however, is remarkably good and 
the variation between pilots astonishingly small. (2) The inference, of course, is 
that objective factors are coming to play and it is there that the prospect of a pos- 
sible quantitative theory for handling qualities arises. 

The experiment was conducted in the compensatory arrangement as shown 
in Fig. 2-1 with the plant (aircraft) dynamics having the generic form (short-period 
approximation): 

5(1+4+0. Bains 
G = —— (5-1) 
P s(s*+26w_ s+w*) 
n n 

Hall used twenty combinations of € and w, to test pilot response. Of 
interest here is a representative sample covering the spectrum of pilot opinion. 
The configurations designated by Hall as C30, €25, C32, C31, and C21 are to be 
studied, Table 5-1. 

The following pertinent performance comment is a quote from the Hall 


report. ) ''As a measure of tracking efficiency... an average absolute value 


1) Cooper, G.E., Understanding and Interpreting Pilot Opinion, Aero. Eng. Rev. 
won lo, No, 3., Mar., 1957. 


(2) Ashkenas, I. L. and D. T. McRuer, A Theory of Handling Qualities Derived 
from Pilot- Vehicle System Considerations, Aerospace Engineering, Feb., 1962. 
(3) 


Hall, 1.A.M., Effects of Controlled Element on the Human Pilot, Princeton 
Univ., WADC Tech. Rept. No. 57-509, Aug., 1958. 
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TABLE 5-1 


Opinion Rating 


” C aeoee Pilot | Average ce a Average 
umber n IT 
apet > pst ; 
= ce Ce 
7 a 


criterion was preferred to a root mean square error, since it was believed that 












Configuration Mean Absolute Tracking 





Error |€| (degrees) 


W 






coe 
= 





the occasional large errors, which essentially decide the value of RMS error, 
more often were due to lapses in concentration rather than conditions associated 
with longitudinal dynamic characteristics... In discussing average error it is 
interesting to keep in mind that without any pilot output this average would be l. 36 
degrees... due to longitudinal... excitation signals... Average longitudinal errors 
of less than 1. 36 degrees indicate that the pilot is doing some good—something 
which must not be taken for granted. "' 

Additionally, the system excitation was random appearing with Gaussian 
amplitude probability distribution which resulted ina realistic target course; 
specifically, the input was white noise cut off at 18db/oct after 1 rad/sec. Itis 
noteworthy that "attempts to increase the cutoff frequency to 2 rad/sec. and 4 
rad/sec. had resulted in loss of the pilot's faith in his ability to ever catch the 


(1) 


In the case of C21, a Routh test shows this configuration to be closed-loop 


target. "' 


unstable (without a human controller in the loop). 


1) Hall, I.A.M., Effects of Controlled Element on the Human Pilot, Princeton 


Univ., WADC Tech. Rept. No. 57-509, Aug., 1958. 
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the two changes of sign in the first column signify two positive roots and therefore 


closed-loop-system instability. The human operator's task is clear if any attempt 


is made to control this system. 


Next, the closed-loop-stability characteristics of C32 in the Hall report 


are determined. 


1+ F(s) ie 
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The Routh array is given by 


Ss 2.64 


1] -(71-150) _ 
TS 


59 | -(-89)(71) 7 + 


Sein. 


at 


since there are no sign changes, there are no positive roots, and the system is 








closed-loop stable (without the human operator). It should be expected that less 
human effort is needed to control this system. 

In plotting these systems in conjunction with the variable pavameter map- 
ping representing the human adaptive range, one hopes to see the -L AG w) contour 
for C21 cross the ) region in a less acceptable manner than the plot of -L(j w} 
for C32 crosses it. The term "a less acceptable manner" may not be adequately 
defined at this point, but the two extremes of stability are clear, i.e., the case 
of stability regardless of the pilot's presence, or the other extreme of instability 
irrespective of the pilots attempts. 

The configuration C21 and C32 are plotted in Figs. 5-1 and 5-2, respec- 
tively. The solid line through the specified frequency points is merely a suggested 
locus and is offered mostly for visual ease. As indicated in Chapter Four, the 
important properties of each system are interpreted from the characteristics of 
the difference vector drawn between all frequencies of the system locus, -L: 
and corresponding frequencies in the variable parameter mapping, 2. The dif- 
ficulty, of course, is just what sub-region, set of points, or point is used in the 
large expanse of (1) ? It is recognized that if a point in the space of the ] map- 
ping is selected this is equivalent to selecting a set of values for the triple 
(T, Tp Ty): It might seem that this amounts to using a specific set of parameter 
values and working with a nominal transfer function, so the human adaptive range 
is not really used inthe study. This is rather a pedantic point of view. It is seen 
that all the possible capabilities available to the human (by way of parameter 
variation) are clearly displayed in the plots and the net effect of choosing any set 
of them (a triple) is equally visible. The usefulness of such information is un- 
deniable, particularly when experience and experiment allow correlation of cer- 
tain of the sub-regions of the ] space with human effort. Such implication is 
suggested even by the few examples viewed here. In any case, selecting a set of 
parameter values corresponding to a point in ] space, or even a cluster of para- 
meter values (T,. Tp Ty): each over a tiny range whose image yields a sub- 
region in { space, is really no more than the human operator is forced to do. 
The human inability, or reluctance, to change parameter character (break fre- 


quencies) discontinuously, corresponds to a constraint upon jumping about in ? 
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Fig.5-1: Plot of Configuration C2) 
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space once some "neighborhood" or sub-region of that space has been chosen. As 
to how such a neighborhood or sub-region is chosen, some insight regarding this 
devolves from the examples along with the inherent assumption that the human 
operator will always try to dohis "best''. It is not unreasonable to associate 
"doing one's best" with the maintenance of stability, which is the prime require- 
ment of any system. Once in this sub-region it seems entirely proper that signifi- 
cant departures from it do not occur and interpretations as to whether any portions 
of this sub-region are ''enclosed" or not by -Lo are possible. 


From (5-1) is obtained the FIXED portion of the system, the negative of 


5(1+.6s)w0~ 
ee (5-2) 
= s(aee 2 Ow Sa 
n n 


(5-2) is plotted to obtain -L_ (180° phase angle is added). Since the delay time, 


-TS 


e , and human gain, K,, are lumped with (5-2) it should be borne in mind that 


h’ 
Fig. 5-1 actually represents a situation where K,e | °=1, i.e., 7 =0, KF LV; 
The significance of these conditions comes to focus shortly. Equation (5-2) may be 
written | 2 
-5(1 +0. 6jw)w 

i ee eer oe (5-3) 
‘ sale wet 2 Cu jw) 

J n n 


Equation (5-3) is rationalized to give 


10 Cw w (1 tO e G4 ato af wlut - al +0. 6jw) 


-L a_i > A ae een 
n 4Cewe wet Pie WP) 


(5-4) 


n 
The imaginary part of (5-4) is set equal to zero 


2 
»(6Cw -5)4 Beto 


n 


which for C21 (€=0. 2, w= 1.57) gives a real axis crossing at w=+1,99. For C32, 
the solutions are imaginary, reflecting no real-axis intercepts. The conditions 


are portrayed in Fig. 5-1 and 5-2, respectively. 


Transparent overlays of Figs. 4-19 through 4-23 are enclosed in an envé- 





lope fastened to the inside back cover of this study. These overlays, which are 
universally applicable, for the given scale, facilitate considerably the determina- 
tion of system stability characteristics. Detailed consideration of configuration 
C21 follows. 

A glance at Fig. 5-1 and the general trend of the Q mappings for the fre- 
quencies given (w=0.5, 1.0, 2.0, 3.0, 10.0), indicates that system stability is 
determined by the relation between the “Ls plot and the overlays for w=2.0 and 
3.0. This is obvious, since the -L_ plot leaves the © region of criticality for 
wd, thus implying that a condition of instability may pertain. It is not yet clear 
how to take the difference vector between -L (52) and ()(j2), i.e., it could con- 
ceivably be connected from -L _(j2) to any point within or on ((j2), and for in- 
stability it must be demonstrated that a clockwise encirclement must occur as 
w-®° , irrespective of this vector's selected orientation (all the Hall configurations 
were open-loop stable). Due to the time invariance assumptions, however, once 
a settled condition of operation has been adopted by the human controller, any 
selected orientation" of the difference vector defines the sub-region in space to 
be used for the remainder of the stability analysis. One notes, for example, that 
if the human operator assumes a set of parameters near point B’ of Fig. 4-21 


(corresponding to a neighborhood of max T,» max T., max Ty): there results, 


from Fig. 5-1, an encirclement by the difference vec and instability. If the 
point F‘° is assumed, the human can “hold out"! longer before enclosure, but 
eventually an encirclement again occurs. Unfortunately, the prospects of being 

at F’ are not encouraging since, although the point corresponds to max Ty (lead) 
which is obviously needed, it also implies T= 9 which, giving some credence to 
the physiological correlations, corresponds to no neuromuscular lag. It is soon 
discovered that if there is any hope for system stabilization or controllability, the 
human operator must assume parameter values somewhere between arc Gi @, 

and arc A‘D’ and inthe vicinity of the origin of Fig. 4-21, also Fig. 4-18 (detail). 
Particularly is this the case whenthe human time-delay is imposed upon the sys- 
tem. Thus, the characteristic of C21 projecting beyond the © region is further 
aggravated by the introduction of the time lag, i.e., the point -L (32), at approxi- 


mately 5.0 on the real axis, rotates with that distance as a radius arm, a number 
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of degrees proportional to the time delay. For example, if T=.15 is assumed, 
the angle of end? (its magnitude, of course, is unity) for w=2 is (2)(, 15) = 
Me@radian = 17.2-, taken in a clockwise direction (lag), and similarly for other 
pertinent frequencies; clearly, conditions are not made better. It might be added 
at this point that the effect of the ew a term should not be considered as adding 
ever increasing phase, thus further complicating the plot; since the human- 
operator transfer function is not valid beyond a maximum frequency of w * 20, 
the corresponding phase angle contributed by this frequency can be regarded as an 
upper bound. 

Since C21 is more than a hypothetical situation, it is possible to check the 
implications of the above analysis. A "best fit'' human-operator description was 


calculated by Hall from the empirical data and a plot of this function can be made. 


The Hall fitted transfer function was (Pilot II): 


s,2 
0. 1(1+5) 15s 


H, = ——45 e™ (5-5) 


Several properties of (5-5) should be evident. First, this is one of the cases where 
Hall required a double lead; this characteristic is not accommodated by the gen- 
eralized human-operator transfer function assumed for this investigation, Equation 
(2-1). Nevertheless, the description is analyzed in order that the curve repre- 
senting the Hall fitted transfer function can be observed relative to the regions of 
criticality which have been used in this study. To be noted are the value of human 
gain assumed, i.e., K, = 0. 1, and time delay tolerated, T=.15; the system 


dynamics, -L are plotted for K,=1 although here it is KNOWN that KL actually 


was 0.1. The effect of this gain Scie must be imposed upon the system 
dynamics as indicated previously; in an actual study, where the approach is from 
the direction of synthesis, the plots imply that reduced human gain is a necessity 
for this man-machine system if there is to be any chance of controllability. 
Equation (5-5) represents those characteristics assumed by a real life 


human operator and when mapped properly, indicates those points in () space 
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actually assumed under the task conditions imposed by C21. Table 5-2 displays 


the pertinent information for 


Se? 
Cl4 0). 5 )a) 


TABLE 5-2 (C21) 





The last row in Table 5-2 displays the values of Hall's ''best fit" transfer 
function for several pertinent frequencies. A plot is not required to see that the 
human controller did, indeed, assume those characteristics previously implied 
by the analysis; furthermore, it is obvious that his capacity is being strained 
severely just to maintain control of the C2l system. That is to say, precious few 
points inthe variable mapping space {) are not enclosed by the plot of Fig. 5-1, 
particularly after introduction of the time-delay effects. It is also interesting to 
observe that the double lead does not impose any gross effect upon the experi- 
mentally determined © that is not already incorporated in the generalized map- 
ping (3-12). Additionally, it is noted that redveed gain, Fig. 5-1, would ease the 
instability tendencies considerably. 

Since the dynamics of C21 were rated 'unsatisfactory'' by the pilots in the 
Hall study, it is interesting to look in detail at an analysis similar to that above 
for C32, to which a rating of ''acceptable'' was assigned. In Fig. 5-2, the plot 
for C32 of the Hall study, it is clear that controllability should not be as difficult 
atask as presented previously by C21. It should be obvious what the trend is for 
the plots of Figs. 4-19 through 4-22 (the transparent overlays) as the frequency 
continues to be increased beyond w=3, Fig. 4-23. The plot for w=3, Fig. 4-22, 
is already so large as to preclude displaying the entire hull upon a page. For this 
reason, only the region of greatest significance is displayed in the plot. Thus, 


observation of the characteristics of C32, Fig. 5-2, should indicate that the human 








bel 


controller need not reduce his Ty to almost zero while maximizing his lead, thus 
operating in the neighborhood of the origin of Fig. 4-22, as was required in C21; 
but, that there is room for assuming parameter values in Q that yield a critical 
sub-region in the neighborhood of 3.5+jl, Figs. 5-2 and 4-22, If the mathe- 
matical modeling here represents in some small way the actual physical situation 
that obtains, it is to be expected that the human controller would feel less vitriolic 
about this system. 

For the system dynamics of the Hall report, C32, the ''best-fitted'' trans- 
fer function for the pilot was (Pilot II): 
H = 0,5 abet epee 


Oo Ss 
(1+ 9q) 


Ee) 





Again the value of gain actually assumed by the pilot was less than that used in the 
plot of Fig. 5-2; the influence of delay time T =0.15 seconds also is not included; 
both facts must be incorporated in the final interpretations. The associated map- 


ping function for (© is 
=e, 25) 
= Seon aa (5-8) 


The form of equation (5-7) may seem to imply a T y= 0 which would te physio- 

logical nonsense. Frequently, experiments connected with tracking studies and 
the servomechanism approach have incorporated some of the minor lags into the 
(1) 


delay time wherever such a technique did not prejudice the data. The pertinent 


data for (5-8) is given in Table 5-3. 


TAB ME 5-3 (C3A} 
pe hoe fo fe fe fot foe 
(jw) fa. 298 3.25 1.844 51.09 |2.594+j1 | 2.9 +). 82/3. 18+j. 40 43. 22+ j. 29 
i] 


, McRuer, D.T. and E.S. Krendel, The Human Operator as a Servo System 
Element, Journ. of Franklin Instit. ,Vol. 267, Nos. 5 and 6, May and June, 1959. 
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The results of actual human control shown in Table 5-3, viewed in con- 
junction with Figs. 5-2 and 4-22, again bear out the previous implications of the 
analysis. It is seen that the human controller did indeed assume those values 
for his parameter characteristics which the two plots of Fig. 5-2 and 4-22 appeared 
to indicate. Since not as large a value of TL was required to stabilize this sys- 
tem, the human operator did not assume any larger value than needed when the 
time lag is considered, and since experiment has shown that generating lead 
(anticipation or derivative error control) is regarded by the operator as a strain, 
the system C32 was regarded as superior to C21 and rated ''acceptable". Again, 
it is obvious from Fig. 5-2, that reduced gain will tend to stabilize the syste m,.this 
was recognized empirically by Pilot II in the experiment as demonstrated by the 
gain actually assumed, KL = O75. 

The plots for system configurations C30, C25, and C31 of the Hall study 
are shown in Figs. 5-3, 5-4, and 5-5 respectively. It is noted that the shape 
of the curve for the particular system being considered, as well as the relative 
location of its frequency markings, lead to estimates of pilot-derived handling 
evaluations with orderings equivalent to those appearing in Table 5-1 under the 
column Opinion Rating (average). A quick check of Figs. 5-3 through 5-5 and 
their relation to the mappings on the transparent overlays corroborates this point. 
Thus,C30 is regardedasa quite good "aircraft"', with C25 and C31 rated lower, 
in that order. The significance of these facts and their ramifications as discussed 
in detail for C21 and C32 follow directly. 

Any person conducting research in the field of control systems and the 
human operator soon encounters a preponderance of studies, published and un- 
published, using a nominal value for TL of approximately 2.0; indeed, many 
experimenters feel that human-generated lead cannot exceed this figure. It is 
informative if the effects upon the stability hull of assigning such a bound to the 
range of T, are considered in some detail. As in Chapter Four changing the 
upper bound on the range of TL has no effect upon the associated circles them- 
selves (in the complex plane), but only upon the extent of the arc segments. Thus, 
fora ThL= 2, one expects to see smaller arc lengths than for a T= 5.3, which 


has been used up to now. It is possible to infer already the effect upon the sta- 
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bilizing ability of the human operator without actually plotting the new specifications. 
One notes that points 5 '6'8°7' corresponding to face 5687 in three space, Figs. 
4-18 and 4-16 respectively, do not move along their circular loci to the extent 
previously, when es 5.3. Thus, line 7’- 8’ makes a steeper angle with the 
real axis and there are many more systems with amplitude-frequency character- 
istics which leave the hull and tend toward instability. It is interesting to note 
that line 7’- 8° relates toa Ty 0 which, as mentioned above, is usually not 
possible (unless the effects of Ty can be incorporated into an expanded T); 
hence, for Na 0.7, the stability criterion becomes associated with line 5 ae 
and further difficulties with stabilization are implied. It is in this way that the 
designer is given considerable insight into the effects of human-operator para- 
meter variations upon overall system stability and even, perhaps, handling qual- 
ities. For example, with maximum lead in, attention may be confined to face 
5687 corresponding to the augmented triangular mapping 5 ‘6 "8° 7", Clearly, the 


' 


outer bound of this ''sub-hull'' most conducive to stabilizing a system is line 


eS i.e. , relating to max Ty and min Tw It is plain to see now why all 
previous examples drawn from empirical results had stability determination 
dependent upon points clustered in a very reduced area of the large region °) , 
i.e., this reduced area corresponded to those extreme values of possible para- 
meter variation most conducive to stability. The consequences of various modi- 
fications in human-operator capability are all quickly discerned with this tech- 
nique, i.e., the region of criticality for stability, if Th 0, is the mapping of 

face 1375 which is quickly found to have a hull 37 5e (with two arc segments), 
and the ramifications of this state of affairs are rapidly seen. * Of course, if the 
ranges of T, and Tw must be altered, the hull in the complex plane is changed, | 
but the new points may be found expeditiously from relations such as (4-12) in 
suitably modified form. The effects of varying T, and Ty Appendix B.2, and 


the characteristic that the Ty circles always pass through the origin, can be 


Appendix C contains an example of the use of the degenerate hull in a two para- 
meter mapping problem as well as an instructive demonstration in stability 


interpretation. 





Way, 


used to good advantage in plotting any modified range on the parameters 
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ee [he Smith Study 
In September, 1963 Ralph Smith reported results from an experimental 
investigation attempting to define the limits of manual control for unstable- 
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second order controlled elements. 
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There are but two results given: 
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Both of which are for the controlled element having the form 


es) Cae eee 
G (5-11) 
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To plot Set (5-9), from (5-11) one obtains 


300(6. 25) : 1875 
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For s = j¥, and after some manipulation, it is found that 


2 
a ° 1k \ -]j ° 

a: 11718 1. Eee 3 5u (5-12) 
(6. 25-1)" + 12, 25 


and, in a similar manner for (5-10), but in a non-rationalized form 


3 
Be 2 = Oe. x10 (5-13) 


n (51.9 - w*) - 55. 04 


Qualitative plots of (5-12) and (5-13) appear in Figs. 5-6 and 5-7, re- 


spectively. The actual coordinate information for (5-12) and (5-13) appears, 


(Mc ith, R.H., An Experimental Determination of the Limits of Pilot Con- 
trollability for Unstable, Oscillatory Second Order Systems, Thesis for degree 
of Aeronautical Engineer, Univ. of Cincinnati, 1961 (also appeared as a paper, 
On the Limits of Manual Control, IBEE Trans. on Hum. Factors in 
Electronics, Sept. 1963.) 
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for selected frequencies, in Tables 5-4 and 5-5. A glance at the row for 
-L_ in these tables explains why Figs. 5-6 and 5-7 are presented qualitatively, 
i,e., the scale is such that the universal overlays, given with this thesis for 


the () regions, cannot be used. 


TABLE 5-4 
~246-j164 }-78-j243 | 43. 8-j167| 62. 2-j89. 2154. 8-j51. 2 








TABLE 5-5 





In the face of an open-loop-unstable plant the plots must be interpreted 
with even greater care than before. The successful application of the variable 
parameter mapping technique to systems having open-loop instabilities may be 
of direct use in extending this study to the monitoring situation where the 
human operator is called upon to control a destabilized plant (or plant plus 
compensator). In the normal Nyquist plot, if a system is open-loop unstable, 


e.g., has two right-half-plane poles, then for the system to be closed-loop 





> aa <a. ie. +h 





We 


stable the net number of encirclements of the critical point must be exactly 
equal to the number of those right-half-plane poles and the encirclements 
must be ina CCW sense. This results directly from the number of encircle- 
ments yielding the difference between the number of zeros and poles of the 
closed-loop-system-function denominator (characteristic equation when set 
equal to zero). Here, zeros (CW rotations) are regarded + and poles (CCW 
rotations) are regarded -. For no zeros to be present in this example, the 
net number of encirclements N = Z- P must be N=0-22=-2. That is, for 
stability, the plot of the so-called ''difference vector'' as s traverses the 
entire right-half plane must undergo two complete counterclockwise revolutions 
about its origin (or, equally well, one CCW revolution for the positive fre- 
quency branch, since there is symmetry about the real axis). If N =0, there 
are also two zeros in the right-half plane and the system is closed-loop 
unstable. What does a CCW encirclement appear as in the dual Nyquist plots? 
The single parameter variation as in Chapter Four and Figs. 4-8 through 4-12 
may be viewed to gain some insight into this question. Specifically, suppose 
the situation were as depicted in Fig. 5-8, where there is shown a hypothetical 
plot of some -L. Clearly, if ae ae , there is a CCW encirclement, while 

if ae < “ane there is no encirclement. If attention is now turned to Smith's 
results, it can be seen from the data points of both unstable plants that there 
are no encirclements, and since Smith used an unstable second order plant 
(two right-half-plane poles), it appears that this is an indication of two zeros 
in the right-half-plane and closed-loop instability. In helping to interpret the 
somewhat enormous plot of Fig. 5-7 (¢ = -0.7), the following quote from (1) 
appears pertinent. ''Controlled element gain had no apparent effect on the 

ope rator's ability to stabilize the system. Changes in K Ks in this thesis) 
from 300 to 900 were barely noticeable by the operator atW = 2.5 rad/sec."' 
For the ¢ = -0. 35 case, the following comment appears. ''The controlled 
element gain, K, (K_ inthis thesis) was originally set at 200. The final 

value was 1200. Changes in gain produced no apparent change in the operator's 
technique. The only noticeable effect produced by gain changes was to decrease 
the maximum control deflection.'' Since the Smith study addressed itself to 


the control problem presenting extreme difficulty to the human operator, 1.e., 


Ves mith, R.H., An Experimental Determination of the Limits of Pilot Controlla- 


bility for Unstable, Oscillatory Second Order Systems,Thesis for degree of Aero- 
nautical Engineer, Univ. of Cincinnati, 1961 (also appeared as a paper, On the 
Limits of Manual Control, IEEE Trans. on Hum. Factors in Electronics 

Sept. 1963.) 
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controlling an unstable-open-loop system, one further quote from (1) is appro- 
priate. ''Observations and operator comments during the conduction of the 
experiments for © = -0.35 indicate the appearance of a very important pheno- 
menon. In watching the operator at his control task, it seemed as though his 
reactions were intuitive. The changes in error were so fast and the error rates 
were so great that satisfactory control looked impossible; yet, somehow, suffi- 
cient stability was generated by the operator. After becoming accustomed to 
the control task the operator often remarked that he was simply outguessing 
the system. He made the statement that he was attempting to get his stick 
motion in phase with the error and trying to keep it there, i.e., system 
stability was his foremost requirement -- keeping error minimum was not 


consistent with maintaining stability. "' 


The task of interpreting Figs. 4-6 and 5-7 (actually the data of Tables 
5-4 and 5-5 in conjunction with Fig. 4-22), must now be considered; it is 
obvious that to stabilize the systems all the human gain must have been quite 
low. Unfortunately, the Smith study did not provide information relating to 
KL since Smith's stability criterion was taken to be positive phase mergin for a 
h the calculation technique is not 
clear. The upshot, however, as recognized also by Smith, is that the general- 


corresponding amplitude ratio in terms of K 


ized transfer function for the human operator (2-1) is a model incapable of 
stabilizing the systems under consideration. Yet the human operators in the 
experiment were able to stabilize the system (5-11) for the conditions (5-9) and 
(5-10), though just barely. That is, these results are supposed to reflect the 
"limits'' of manual control. The plant characteristics (5-9) and (5-10) define 
the boundary of marginality, i.e, that point at which the human operator can 
no longer control the system, for a plant having form (5-11). Greater confi- 
dence accrues to the predictive capability of this theoretical approach if it too 
displays such marginal conditions. In order to account for the experimental 
results, Smith hypothesized a second lead term for the human operator model 


so that stability results from the theoretical description, also. This tack is 


: Smith, R.H., An Experimental Determination of the Limits of Pilot Con- 
trollability for Unstable, Oscillatory Second Order Systems, Thesis for degree 
of Aeronautical Engineer, Univ. of Cincinnati, 1961 (also appeared as a paper, 
On the Limits of Manual Control, IEEE Trans. on Hum, Factors in Electronics, 


Sept. 1963.) 
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not taken here. It is sufficient to demonstrate that the techniques of this 


thesis imply a marginal stability situation. Thus, for“ =10 and K. =0.1 (not 


an infrequent value for human-operator gain) from Table 5-4, -L - 1, 8-j. 7, 
and it is conceivable that if KY were slightly smaller than 0.1, es -L. locus 
might intercept the small semicircle (which is not a function of frequency, 
Table 4-3) projecting below the real axis on all the 9 plots. Under such cir- 
cumstances, there is an enclosure and concomitant CCW encirclement indi- 
cating closed-loop-system stability. If it is felt that this discussion affords 
nothing more than an amateurish job of rationalization, it is also evident that 
the plots seem to indicate no clear cut conditions of one kind or the other. 

This substantiates that the mapping analysis adequately depicts the marginality 
of the situation. Application of the variable parameter mapping technique to 


the Smith data certainly demonstrates that the human operator is being pressed 


to the limit of his abilities. 


The possibility of applying the techniques of this thesis to the multi-axis 
Situation should not be discounted, i.e., the study provides a theoretical method 
of approaching the multi-axis situation where a particular ''nominal" transfer 
function may not apply with invariance over all channels. The plots display the 
entire adaptive range of the human-operator transfer function and, since the 
bandwidth dealt with is so limited (0< ¥< 20), the number of plots required to 
cover the necessary frequency range is pleasantly restricted. This may prove 
also, to be a more expeditious and fertile approach even for the single axis 
problem, since the automatic part of the system and the human operator are 
displayed separately. This isolation permits observation of the control charac- 
teristics of the automatic part of the system from the same plots used with the 
human controller analysis. Indeed, the positioning of the -L locus relative 
to the (+1+j0) point gives precisely the information normally obtained in a stan- 


dard Nyquist analysis with respect to the (-1+j0) point. 


No experimental monitoring studies as such have been available to the 
author and none exist to his knowledge. Information vital to multi-axis moni- 
toring by human operators may devolve from the Elkind studies of human 


response to abrupt plant changes, i.e., sudden changes in the gain and polarity 








of a controlled element that was simply a gain; however, considerable experi- 
ment with Tang coos situations is needed to generate confidence in any 
1) 4 


inferences. The theory presented herein works remarkably well for the 


single-axis situation and an equivalent simplified monitoring topology with 
open-loop instability, i.e., the theory predicts results consistent with the 
empirical evidence. Experimental results are needed to determine if the 
theory is totally applicable or, at the least, can continue to provide a valuable 
tool for design insight and ''rule of thumb" criteria in multi-channel task 


situations. 


The difficulty with the experimental studies reviewed here is that all 
parallel compensators, G.: have been equivalently zero, Fig. 2-2. As indi- 
cated above there are no experimental results for monitoring. The question 
is: can theory predict that the human operator will have a less difficult task 
in a monitoring situation and how can this be determined? It would be instruc- 
tive to observe an Ho controlling aG_ with G. = 0, then compare the human 
effort there with that required when an A is paralleled by a G. designed for 
good automatic control of the same on Such information is not available at 
present. The problem is compounded by the valid question as to what role the 
human would play when the system is performing properly, i.e., presumably 
the human enters only upon some catastrophic failure of one kind or other. 


Obviously, considerably more research is required in this area. 


55 Summary 


The variable parameter mapping of the three time constants (which 
man can adapt to a given task) in the human-operator description, is used in 


conjunction with a slightly modified Nyquist plot of the plant characteristics to 


: Young, Green, Elkind and Kelly, Adaptive Dynamic Response Characteris- 
tics of the Human Operator in Simple Manual Control, IFEE Trans. on Hum. 
Factors in Electronics, Sept. 1964. 


Bander V.A., Predicting Pilot Ratings of Multi- Axis Control Tasks from 
Single-Axis Data, IEEE Trans. on Hum. Factors in Electronics, Sept. 1963. 
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demonstrate that stability of the closed-loop system (with a man in the loop) 
can be predicted. Experimental verification of the theory is obtained by 
analyzing certain systems which have been studied empirically in the past and 
for which results have been published in the literature. To this end, trans- 
parent overlays are provided which facilitate stability determination. The 
theoretical predictions parallel the experimental results remarkably well; 
one of the sets of data relates to a system which is open-loop unstable. The 
correlation possible between boundary lines of the mapping and the physical 
capacity of a human operator is cited and, in the cases of stable control, the 
margin of human capacity remaining (displayed on the same diagram) is 


indicated. 
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CHAPTER SIX 
PERFORMANCE 


6.1 Introduction 


As in the discipline of inanimate systems engineering, several performance 
measures have been used for man-machine systems to provide a so-called figure 
of merit or index of quality and, as in the case of the all-hardware system, defi- 
ciencies have been found to exist with all these scalar indices. The difficulty in 
attempting to quantify a complex system's performance by virtue of a single num- 
ber is that satisfying the requirement on this unique scalar may still result ina 
system with other unsatisfactory characteristics. Of course, the designer of per- 
formance measures answers that all pertinent characteristics are duly weighted 
and included in the scalar index thus causing it to reflect those properties deemed 
important. Still, the final selection of those properties regarded as necessary for 
an acceptable system and assignment of relative weights are only two of many rea- 
sons for viewing single performance measures with some caution. Generally, the 
selection of performance index has often been predicated more on ease of calcula- 
tion and other conveniences than upon how accurately the truly efficacious char- 
acteristics of a system are reflected. In the psychological literature, time on 
target, an easily obtained number has been used extensively although studies have 
shown its deficiencies; it has even been demonstrated that there is considerable 
danger in inferring that system-measure changes imply system-component 
changes. me) Other indices used by psychologists studying manual control have 


been average error, average-absolute error and root-mean-square error; attempts 


1 Bahrick, H. P., P.M. Fitts, and G.E. Briggs, Learning Curves-Facts or 


fapeipacts, Psych. Bull. , 54, 256-268, 1957. 


(2) 
Taylor, F.V., and H. P. Birmingham, That Confounded System Performance 


Measure-A Demonstration, Psychological Review, 66, 178-182, 1959. 
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have been made to measure task difficulty but independent critics have dismissed 
them as measuring little about the specific requirements of a task. J The engi- 
neering approach to the performance measure quandary has been to employ the 
seemingly dependable (and quite convenient) mean-square error, as well as cer- 
tain other control theory measures. The period of time over which the measures 
themselves are calculated serves to influence the result and thereby confuse the 
problem still further. Historically, pilot opinion has been used as an acceptance 
technique in the case of aircraft flight control systems. In this regard, the fol- 


lowing quote is considered appropriate. ‘''It should be clear at this point that a 


single measure for a complex tracking situation cannot give information concernin 
(2) 


both system performance and operator performance'' (underlining theirs). 


In this chapter, the various performance measures are discussed anda 
brief survey is given of some techniques available for "optimizing" a system so as 
to best satisfy these indices. The formulation of a ''vector'' performance index is 
considered as a means of circumventing certain of the disadvantages of scalar 
indices. The concept of a "'best-overall'' compensator (fixed compensation), for a 
plant with variable dynamics, is introduced and a technique for specifying the para- 
meters of this compensator is offered. This technique depends upon the formation 
of an 'artificial'' function which, when minimized to zero, identically satisfies a 
set of performance requirements. A "direct search" procedure, used in conjunc- 
tion with this ''artificial'' function, which is sequentially modified, serves to isolate 
satisfactory values for the compensator time constants. In this manner, the space 
of acceptable compensator time constants is defined and correlated with the space 
of required performance. The entire procedure is especially amenable to the 
electronic computer. An example of the search technique is given at the end of 


the chapter. 


Taylor, F.V., Human Engineering and Psychology, in ''Psychology: A Study 
of a Science", Koch, S., (Ed.), Vol. 5, McGraw-Hill, 1963. 


Summers, L.G., and K. Ziedman, A Study of Manual Control Methodology 
with Annotated Bibliography, NASA Contractor Rpt. CR-125, Nov., 1964. 


(2) 
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6.2 Varying Parameter Plants and ''Best'!' Performance 


The question frequently arises how "best" to use the human operator in 
any given system; it 1s assumed that the stability criterion is satisfied. Given 
the requirement for a human operator one may ask, alternatively, how to design 
a satisfactory compensator for use in a system which includes this human con- 
troller? The design task may be complicated further by specifying a plant which 
itself contains some varying parameter. Generally, the controls engineer has 
tended to avoid computer design for the more familiar Bode diagram, frequency 
plot or root locus technique. An approach to the performance problem, is pursued 
in this thesis which, it is hoped, can be suggestive of a standard technique, or 
algorithm, for handling far more severely complex situations than can be treated 
here. The method is particularly applicable to systems containing human con- 
trollers since the human adapts his characteristics to some unknown internally- 
generated criterion (which is frequently taken to be some function, or functional, 
of the error). The situation is quite fluid in the face of systems with varying- 
parameter plants since the error itself then changes due to this parameter varia- 
tion thereby causing further adaptation on the part of the human. It is felt intuitively 
that some design compromise on an adequate compensator should be possible so as 
to provide the ''best'' performance over some range of parameter variation. The 
question immediately arises regarding a proper criterion (performance index); 
this is a subject around which much controversy still rages within the engineering 
profession. pa) Frequently, the physics of the situation suggests a performance 
index which leads to a mathematically intractable problem. The designer then 
must compromise on a choice of performance index which permits a meaningful 


solution, i.e., does not preclude analysis,yet still provides a useful quality mea- 


1) G.C., L.A. Gould, and J. F. Kaiser, "Analytical Design of Linear 


Newton, 
Feedback Controls", John Wiley & Sons, New York, 1957. 


(2) Gibson, J.E., ''Nonlinear Automatic Control", McGraw-Hill, 1963. 
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sure for the system. Usually, the indices yielding the greatest analytical ease in 
their calculation are the mean-square error and the integral-square error. For 
inputs having Gaussian characteristics and certain other properties which are 
frequently encountered (or at least frequently assumed) it can be shown that sev- 
eral error criteria offered as more appropriate lead to the same results as the 


(1) 


mean square error criterion. Of course, the advent of high speed computation 
and concomitant expeditious numerical integration offers great latitude in the 
choice of a performance index. 

Some work has been done indicating that minimization of the mean-square 
error is a fair approximation to what the human controller attempts, i.e., it 


(2) 


appears the ultimate effects are similar. Generally, the mean-square error 


is used as a criterion for stochastic inputs; indeed, this is the basis upon which 


(3) 


the human-operator transfer function has been generated. For transient inputs, 
the analyst frequently uses the integral square error criterion; of course, the 

ISE is used only in the design of systems having zero final error since continual 
finite error leads to the non-convergence of the ISE index. With either measure, 
however, stability checks must be made at discreet (also discrete) moments to 
insure that the values of the adjustable parameters which minimize the performance 
index correspond to a stable solution. This is elaborated upon below. Thus, one 
finds that the MSE is used as an index of, loosely speaking, steady state perfor- 
mance while the ISE is used to measure transient performance characteristics. 
Certainly, the selection of a performance index is at the discretion of the designer 
(or the customer) and modern digital computation offers the possibility of numer- 


ically computing heretofore mathematically unattractive but physically useful 


response criteria. 


: S., Non-Mean-Square Error Criteria, IRE Trans. on Information 


Sherman, S., 
Miacory, Sept., 1958 

(2) Roig, R. W., A Comparison Between Human Operator and Optimum Linear 
Controller RMS-Error Performance, IRE Transactions on Human Factors in 
Electronics, Mar., 1962. 


(3) Graham, D., and D. McRuer, "Analysis of Nonlinear Control Systems", John 
Wiley & Sons, 1961. 
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The remainder of this section is confined to the exposition of a possible 
approach to solving such a "'best'' performance problem. The attack is along 
conventional lines while that of the next section is considerably more original 
and exotic. The system studied here is somewhat simplified but it should be 
recognized that the simplification incurs no loss of generality. 

The system under investigation appears in Fig. 6-1. The configuration 
within the dotted lines represents a ''choice" situation, i.e., the human operator 
either chooses to operate upon the error with his dynamics or he does not. Again, 
such a situation where the human operator decides whether or not to enter the 
circuit may be considered as an alternate form of the monitoring situation. The 
representation here does not place the human operator in parallel with the original 
compensation but in cascade. It is recognized that the design of G. for good 
automatic control may be distinct from a good design with a human controller in 
the system; therefore, the human is cascaded always with some G such that 
bey = G. 2a provides good compensation for the human controller-plant com- 
bination. It will be observed that no loss of generality occurs if the situation is 
depicted as in Fig. 6-2 since the plant-human path is under study. The human- 


operator transfer function is considered in the simplified form 


5 


-.ls 


oie ara (o=1) 





where h can be varied adaptively over some range. The unknown compensation 
to be determined is cast in the form of a straight gain, 4, and the plant has the 
form 

G = + — (6-2) 

ees eee 

where 8 varies over some range. It should be noted as the study continues that 
any class of more complicated frequency-sensitive compensation, a plant with 
several varying parameters or a mathematical model for the human operator con- 
taining more adaptive variables, only makes the calculations more severe but does 


not change the general technique. Thus, Fig. 6-2 takes on the form of Fig. 6-3. 
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Fig. w-3 : Human Control of Varying-Parameter Plant 
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Since the study is restricted to deterministic inputs (steps, ramps, etc.) the 


performance index is chosen to be the ISE, 


r 2 
ISE = i= | e°(t) dt (6-3) 
O 


where e(t) is the error signal. 
A design is desired, i.e., a value for 4, such that the ISE is minimized. 
The standard analytic technique for solving this problem uses a rational transform 


e e -7T 
approximation for e ‘ 


— 
i 
= 
47) 


aa gee (6-4) 


Relee 
lo 


' This permits evaluation of the ISE as a function of the free parameters. One may 
then attempt to minimize the integral by taking derivatives of ISE with respect to 
each of the free parameters setting them equal to zero and solving for the critical 


value of each parameter. Thus 


ol. ol. ol. 
iy gcoemeoe: sec Sie 


The step above is often more easily described than accomplished, but the 
main difficulty is not merely that of explicit solution. In the case with a human in 
the system, one must assume there is an attempt to adapt operating procedure so 
as to minimize some function of the system error, irrespective of the value of a 


and 8. The 8 variation is not controlled and may represent many possible 


The evaluation is easily done by means of tables. >) 


(1) Newton, G.C., L.A. Gould and J.F. Kaiser, "Analytical Design of Linear 
Feedback Controls". John Wiley & Sons, New York, 1957. 
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phenomena; e.g., change of response characteristics with altitude by a high- 
performance aircraft, drift in plant parameters, etc. The problem may be ap- 
proached as one wherein Nature is in direct opposition to satisfactory perform- 
ance and a design is sought such that performance is at least as good or better 
than some value of the index of merit, in this case, ISE. For tutorial purposes, 
consider & and 8, each having such range that they are adequately described by 
a grid of five (5) values, each grid covering the respective ranges of @ and 8. 
Thus, for each pair (a@,8) there is obtained that value of h° which minimizes | 


the general performance index S; this index is a function of the three variables 


a,8,h and is taken in ISE form. The minimization of 
s(a, 8, h) fe eat (6-6) 
is then given notationally a 


= fat SEE Er, 8) (6-7) 


i.e., for each (4,8) over the chosen grid 


One finds an h° 
=o SEE eee 


(a, 6) (0, 1) 


01 
O 
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( j) —— hh ij 
This step represents the effect of assuming that the human operator always does 


the ''best'' he can, where ''best'' is the optimization of some error criterion (here 
the identical performance index, ISE, is retained throughout as the criterion but 
this is not mandatory, i.e., distinct performance indices can be minimized for 
h,@ and 8, respectively). 

In the next step, for each value © can take on there is found that value of 
B which maximizes the ISE. Here one observes the result of assuming the 
"cussedness" of the plant; this is, at times, called ''worst case design''’. The 
disadvantage of this approach lies obviously with the extreme conservatism of the 
final solution. Thus, 


2, max | e“ dt —————~ 8(a) (6-8) 
O 


B 
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1.e., for a= 1, run through 8=1,...,5 to get an optimal - 


oo ae T= 
a 


>, run through B=1,...,5 to get an optimal Be 


fo «) 
where the ''optimal'' 8B maximizes je? dt. 
O 


Finally, one looks at the five values of (6-6) for a=1,...,5 and picks the 


minimum of these, thus yielding the "'best''a. Hence, 





foe) 
3. min i ae dt ——$——_——_—. a” 1a° (6-9) 
O 
The above procedure insures that performance 
s(a°, 8. h°) < value i Ae ates Ss(a?, ahr ho (6-10) 
O 


where Be it will be recalled, refers to the ''worst'' situation. For those familiar 
with Game Theory, the use of a game theoretic approach in the technique is notice- 
able, i.e., a minimax solution is being obtained. Clearly, even in this extremely 
simplified example many calculations and evaluations are required and frequent 
comparisons between them must be made. The technique is especially suited for 
digital or other machine computation. 

Two additional points with regard to this approach should be mentioned. 
The first step (the adaptive action) yields an h°(a,8) from which one can proceed 
to evaluate S(a, 8) as indicated above, but care must be taken to insure that 
(a, B)€ Ro} that is, the pair (4,8) must correspond to a region of stability. 
Otherwise, S(a,8B) ~@ if (a,8) ¢ R The parameter space plot of Fig. 6-4 
should clarify this. If the region of stability has the form sketched in 4B space 
and the range on &@ and §8 is indicated by lines (1)-(2) and (3)-(4), respectively, 


then it is clear that during the step max S(a,8) for the value a = a the criterion 


is maximized by picking 8 = ee since this combination corresponds to an unstable 
system. Either a manual check must be conducted to ensure that the optimal 


values determined correspond to a stable condition or the computer program must 
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Fig. 6-4: Solutions Outside the Region of Stability 
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include a check for stability during the computational cycle. 

If the transform of e(t) is a ratio of polynomials then the evaluation of 
ISE is straightforward, by virtue of previously prepared tables,and only the 
minimization is carried out by computer. If e(t) transforms into an irrational 
form (as it does with a pure time delay) then evaluating the ISE is not so simple. 
The integral may be truncated, however, and evaluation carried out numerically. 
To follow such a procedure, however, reduces any feeling of commitment to ISE, 
or other mathematically attractive figures of merit, since their calculation would 
require as much work as that involved in other more physically appealing per- 
formance indices. Also, a coordinated computer program may be possible such 
that, while the performance index itself is being calculated, the information 
generated in that process is used to simultaneously minimize the index and pro- 


vide a stability check. 


6.3 The Concept of a ''Vector'' Performance Index 


A technique is now considered for determining the "best-overall'' com- 
pensator for any given system with a humancontroller. The approach in this 
section may be used for a plant with variable parameters; the only restriction is 
the usual aggravation imposed by the ''curse of dimensionality", i.e., the larger 
the number of varying parameters the more iterations required. (1) The solution 
to be presented depends intimately upon an initial stability study conducted in 
accordance with the procedure previously outlined in this thesis. : Thus, for the 
situation of Fig. 6-5, the result desired is values for K a@ and B which pro- 


vide good, overall performance, defined by specifications on the range of the 


l B R., ''Adaptive Control Processes-A Guided Tour", Princeton Univ. 


ellman, 
Press, Princeton, New Jersey, 1961. 


The author has been unable to find in the literature information relating to the 
step response of a human-operator-frequency-sensitive-plant combination, 
although there is considerable published on such responses for plants represented 
by pure gains. Thus, the study represents a theory which is yet to be checked 
empirically, although all previous experiments would seem to justify this 
theoretical extension. 
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triple (t_, 0. t), in the face of a plant which has a varying ¢. The triple 
(t_ Oo. t .) represents the values of rise time, overshoot and settling time, 
respectively, which follow in accordance with some arbitrary definition. Of 
course, the triple may be regarded as a vector quantity in three seven The 
stability analysis outlined in Chapters Four and Five is first accomplished; for 
the class of compensator being considered, one zero and one pole, it can be 
determined quickly what approximate ranges of values on & and 8 provide good 
stability. The value of K is not difficult to determine since a nominal value of 
KL may be used, thus indicating quickly what the overall loop gain should be, 
from this a first approximation to KO may be calculated. Refinement of this 
value is made along with the search for final values of 4 and 8. Gain adaptation 
is the easiest act for the human operator to perform; indeed, his task is regarded 
as simplest when he acts merely as an amplifier. (1) This aspect of the variation 
therefore can be discounted in the remaining analysis. The competent designer 
desires a final empirical test anyway, if only by simulator, overall operation can 
be double checked at that time. In Section 6.5 an example is given which includes 
consideration of K 

In addition to satisfying the general stability requirement the designer does, 
of course, attempt to use the freedom available in selecting 4 and 8 so that the 
plant's Nyquist curve passes through that part of the ] region corresponding to 
reduced task requirements on the human operator,(Chapter Five). The entire 
analysis so far has been conducted for a nominal value in the range of variable Ce 
Generally, there are some ranges on 4 and 8 which seem to meet the conditions, 


or it is apparent they are impossible to meet under the constraints imposed, ora 





Quality measures other than the three suggested are, of course, possible; the 
technique permits considerable latitude in employing multidimensional per- 
formance measures. 


) Birmingham, H.P. and F. V. Taylor, A Human Engineering Approach to the 
Design of Man-Operated Continuous Control Systems, NRL Rep. 5333, 7 Apr., 
1954 (also published as A Design Philosophy for Man-Machine Control Systems, 
Proc. IRE, Dec., 1954.) 








10 


different class of compensator may be required. With the ranges of a and B 
that provide acceptable stability available, nominal values for the human controller 
can be read from the plot using the time-invariance asssumption. These values 
must also retain the "high-frequency" end of the compensated Nyquist contour ina 
location satisfying the stability criterion, as reflected by the corresponding 1 
region. " That is, the triple (T, Ty) Ty) must have value such that it lies to the 
left of the plant's Nyquist curve and still be within the range of human adaptivity. 
Simultaneously, Th» T, and Tn should be as small as possible; particularly 
Ty» since this parameter is almost a direct indication of 'comfortable'' adaptation. 
All these conditions must also be consistent with physiological realities. 

With the analysis above conducted for a nominal ¢, it is then repeated for 
a gridof ¢ values covering the range of ¢. The intersection of the various 
ranges on © and 8 yields those ranges of 4 and 8 which satisfy stability over 
the expected range of variation in €. These final ranges of acceptable & and 8 
must now provide the specific values of &@ and 8 which yield the best performance. 
Again, a vector (rather than a scalar) definition of performance is being used. The 
requirement is to meet specifications on rise time, t overshoot, oe and settling 
time, o . There is no problem in determining what values of t oe and oe 
result from any given set (a,8); the question is how to go about selecting sub- 
sequent trials of 4 and 8 so that each iteration is at least as good as the last. 
What is sought is some "scientific method" for approaching the sub-space of 
acceptable performance within the space of acceptable stability. The following 
approach provides a routine technique, or algorithm which lends itself to com- 
puter solution. With the resolution of ranges on @ and B which provide satis- 


factory stability characteristics, there is conducted a search for the sub-range 





Inherent in this statement is the asumption that this is what the human does; 
i.e., it is hoped this adequately reflects his unknown performance criterion. 


That range over Q and 8B, respectively, whose end points are within all the 
other ranges of 4 and 8 determined for various ¢ is called the intersection 
(greatest common range of 4 and 8 over C). It should be noted that this 
step considerably reduces the experimental region remaining. 
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which yields good performance, Fig. 6-6(a) and (b). Quite possibly none may 
exist; this situation is discussed later. A grid search of 48 space could be 
performed to find acceptable (t_, Ou: t.), but that is a strictly trial and error 
method. Ifthe (@,8) range is small, straightforward grid search could be the 
most efficient technique; if there are more than two parameters available for 
manipulation, however, the algorithm is probably more efficient. 

The function 


_ Z Z 2 
f = (lee T) ous 0.) tas T.) (6-11) 
is formed where O<t <T,0<0 <O, O0<t <T ,i.e., T ,O, andT are 
r- rr S— § s=— os igre Ss 


the maximum acceptable bounds on t. ) and to respectively. It should be 
noted that the function f has the interesting property that when it is zero the 
specifications are met identically at their upper bounds, i.e., @ and 8 have 
values such that the surface boundary of the acceptable performance region in 
(t_, Ou t.) space has been attained, Fig. 6-6(b). The function is an artifice 
created specifically to facilitate finding some point on this surface and, in this 
instance, is still a function of two independent variables & and 8B, although its 
ultimate evaluation requires three numbers tO, t.- Thus, it is desired to 
minimize f over the space of 4B which yields stability, recognizing that due to 
the sum of squares form of f it cannot be negative, while to meet specifications 
just barely, min f must equal zero. The form of the function therefore precludes 
the occurrence of saddle points when the min of zero is attained. The search 
techniques for ''optimum' a and B must incorporate certain characteristics to 
keep the designer adequately informed. 


First, it is desired to check each of the terms t_-T_,o -O ,t -T 
r rs s’ s Ss 


algebraically; in the event that they are all negative the program should STOP 
and PRINT OUT such information since this is a reflection of an adequate design 
whose response falls inside the acceptable performance region. Second, when 
the search pattern finds min f a PRINT OUT should again be required, since it 
must then be determined by the designer whether the results are sufficiently close 
to zero to warrant an interior search of the acceptable performance region. If 


they are not, this is an indication that the requirements on tO te cannot be 
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met and a different class of compensator may have to be used or other compromise 
fashioned. The utility of such information should not be regarded lightly since it 
represents an analytical demonstration of inconsistent or impossible specifications. 
Too frequently there have been attempts at building systems to meet such impos- 
sible requirements and only successive futile results have led to their abandon- 
ment; any technique which uncovers such sets of impractical demands serves the 
important purpose of precluding the expenditure of funds and effort on an un- 
realizable quest. 

Having located the surface of the acceptable region of performance, one 
proceeds in the ''same direction" until the opposite side of the region is found in 
which case the entire line between (and the associated values of 4 and 8) must 
be contained inside. Many times the "opposite side" is known due to constraints 
on the parameters, e.g., @ and 3 must be 7 0. In this manner, the region can 
be crudely, or completely, defined by successive computation along lines per- 


pendicular to the one originally found. 


6.4 The Gradient Method 


Frequently, a steepest ascent (descent) technique is employed in such pro- 


(1) 


the description allows comparison with the procedure to be given in detail some- 


blems. This type of approach is described briefly, although it will not be used; 
what later. From (6-11) for an arbitrary (a,,8,) the function f (a, By) is deter- 
mined, i.e., for the values (a, By) the system response to the input function of 
interest (here a unit step) is calculated and to. t is found thereby permitting 
f to be enumerated. For some (4a,A8) the functions f (a, +4a, B,) and 


if (a,,8,+ AB) are determined similarly. The function below is formed 


f(a, +44, 3,) -f(a,,8,) 
Aa 


1) Wilde, D.J., ''Optimum Seeking Methods", Prentice-Hall, New Jersey, 1964. 
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and also, 
f(a, B +48) 7 £(a,,8,) 
Ap 


In this way, estimates are obtained on — and St respectively. 
But the gradient of f is 

Vf = ae i, tort 
where i) and my, are unit vectors in the positive 4 and 8 directions, respectively. 
The gradient of f, of course, points in the direction of greatest increase of the 
function f. Since it is desired to minimize f the program proceeds in the opposite 
direction. By surrounding the point (a, 8) by a pattern of, say, two points 
oriented as in Fig. 6-7, the gradient at each of these two points may then be 
determined and iteration continue in an appropriate direction. Thus, the direction 
of incremental change in @ and 6 yielding the greatest decrease in f (or, increase 
in minus f) would be followed to obtain the next trial point. The search procedure 
is then repeated as one follows one's nose to a minimum. When the optimal point 
(be it min or max) is apparently reached it is good practice to conduct a search in 
the local area surrounding this isolated point to ensure that it is not an apparent 
optimum. For this purpose, non-linear exploration of the adjacent region is in 


(1) 


order; Box and Wilson have discussed the problems inherent in such studies. 


6.5 Direct Search and Pattern Strategy. 


The upshot of the previous discussion is that any of several different type 
search techniques can be used to determine the "best'' @ and 8. The approach 


used in this section, and which seems to offer the greatest advantage for the sys- 


tems under consideration, follows that of Hooke and Jeeves. (2) This technique of 


(1) G.E.P. and K.B. Wilson, The Experimental Attainment of Optimal 


Box, 
Conditions, J. Roy. Stat. Soc., B13, 1(1951). 


(2) Hooke, R., and T. A. Jeeves, Direct Search Solution of Numerical and 


Statistical Problems, J. Assoc. Comp. Mach.; 8,2, Apr., 1961. 
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Fig. 6-7: Perturbation Pattern for Gradient Determination 
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direct search with pattern strategy is especially suitable to minimization of a sum 
of squares, but its overriding virtues are its straightforward simplicity and com- 
patibility with the electronic computer. 

The starting techniques in pattern search are quite similar to those of the 
gradient or steepest ascent method, however, it appears at first glance that pat- 
tern search neglects some information that is readily determined and used in the 
gradient method. The effects of this neglect are more apparent than real. By 
"direct search'' is meant a step-by-step examination of trial solutions and a 
method of comparing each to determine which is the "best'' in some sense. Addi- 
tionally, there is some method offered for choosing the next trial point. This last 
aspect forms the strategy part of the search procedure; a ''pattern" approach is 
used and from this is derived the strategy's name. 

The example to be used is shown in Fig. 6-8. It is assumed that the 
ranges of & and 8 have been discovered which yield good stability; specific 
values within that range are now desired which meet the step-response require- 
ments on rise time, overshoot and settling time. The compensator gain Ko is 
included as a third parameter whose value must be determined. The ranges Ko 
a and 8 form a space of points representing possible solutions. Again the func- 


tion, f, is formed 
erode eee eo a yo ili} 
Na TS S S Ss S 


For economy of notation, let the compensator gain, Ko be denoted, K, then the 
point (K,,a,,3,) is said to be better’ than Urea) if 

3 6-12 
The values K,,4,, 3, are determined by adding some increment 6 to each value 
Ky: a and 8 1? the increment size used for each parameter need not be the same. 
This portion of the search procedure is called the exploratory phase. It is under- 


stood that by f (K.,@.,8.) is meant the value of (6-11) found from the substitution 
fea 


a 


“ The inequality is given for a minimization problem. 
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of the t Oo. and un associated with the system step-response when the com- 
pensator has gain and time constants K.,a. and B respectively. The initial 
base point (K), Lore B,) is chosen arbitrarily or the investigator may exercise 

his "hunches". There may be none, one, or more points such that a min f of 
zero is found. The ramifications of this have been discussed above. Under the 
assumption (6-12) holds, there seems little reason in the face of such success, 

to alter the "direction of search"; therefore, exploration is continued by per- 
turbing the independent variables with increments whose proportions are based 

on the number of successes achieved along each dimension in the recent past. 
This strategy forms the so-called ''pattern'"' aspect of the search procedure; of 
course, great strides are made with pattern moves in relation to progress 
resulting from exploratory moves. An example is offered below. In those situa- 
tions where a min f=0 is obtained, the direct search method continues to offer a 
consistent philosophy. The designer unquestionably desires additional information 
with regard to the size of the acceptable performance region and does not wish to 
stop at the barely acceptable surface bound, e.g., other satisfactory triples 

(K., foe B .) may offer simpler design realizability. Direct search can forge across 
this region until it pierces the opposite surface; in this manner the performance 
region may be completely explored. 

Termination procedures are required when a pattern move failure is fol- 
lowed by consistent failure from all exploratory moves. A retrenching is occa- 
sioned by reducing the nominal step size of the increments and new explorations 
are conducted; this is continued until a new pattern is constructed or the step 
size falls below some pre-selected margin. The lower bound on this margin is 
usually imposed by the demands of the particular computation interval peculiar to 
the computer used. Termination does not imply an optimal solution so the wary 
investigator tests the solution generated in radial directions about the final base 
point distinct from the coordinate axes since pattern search restricts itself to 


motion along these axes. 


6.6 "Best-Overal]l"' Compensation of a Varying Parameter System - An Example 


The system for which a best-overall compensator is to be designed is 
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shown in Fig. 6-8. It should be noted that the human-operator transfer function 
has been omitted and a simplified plant assumed; the simplification is imposed 
for clarity only. In this manner, the techniques to be demonstrated are not 
obscured by technical complexities; no loss of generality occurs. The plant's 
damping ratio, €, is assumed to vary Over the range 0.01 < C = 1.0; the rise 
time, t is defined as the time required for the response to go from 10 to 90 
percent of its final value; the overshoot, Ou is the maximum departure from 
final value expressed as a percentage; the settling time, to is defined as that 
time after which the response is confined to a tolerance band of + 5 percent 
about its final value. In order to demonstrate the search technique for more than 
the two parameters @ and 8, the value of gain, K, is included as an unknown in 
the search. 

Sketches of the negative Nyquist contours for the plant, -L. for values of 
GC =0.01, 0.1 and 1.0 are shown in Fig. 6-9. The limiting ranges on @ and 8 
are taken tobe 0<a<5.3, 0<8 < 25. The variable-parameter mappings of 
Chapter Four can be used to aid in the design, since they are not required to 
represent the capabilities of a human operator, i.e., Ho has been assumed equal 
to unity in this example in order to illustrate the search method. Thus, the 
degenerate mapping shown in Fig. C-2, Appendix C, (for w = 1.0) can be used to 
represent the available region of variation for the compensator's time constants, 
a and 8. Similarly, the corresponding regions on the transparent overlays (in 
the pocket of the back cover) can be used to provide additional frequency effects 
other than for w = 1.0. 

It is evident from the overlays, used in conjunction with the more accurate 
plots in Fig. 6-10 (or root locus considerations), that lag network configurations 
maw< (8 ) severely aggravate stability when ¢ is very small, i.e., between 
fedleees- 0.1. Appendix C contains further stability interpretation elaboration. 
From the overlays it is apparent that after introduction of a compensator, stability 
is possible for small ¢ (€#0.01) only if the compensation is lead type. The over- 
lays also indicate that the furthest ''to the left’ the critical points governing sta- 
bility can be placed (for the curves representing 0.01 < ¢< l. 0) is along the 


circular arc segment 3'-7', since this arc is the image of the condition 
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0<a<5.3, 8B =0. Also, in passing it may be noted that the condition 
0<a< 5.3, 8 =0.7 is represented by arc 1'-5' (Fig. 4-18 indicates the definition 
of these points), i.e., B would play the role of Ty in the human-operator trans- 
fertunction and the tiny hull formed by line 3'l', arc 1'-5', line 5'7' and are 7'=3' 
would be of interest, (detail in Fig. 4-18). Since this region remains reasonably 
unchanged in the overlays for increasing frequency, it is possible to use the over- 
lay for s =1,2 or 3, respectively, for an overall estimate of what is needed for 
a,8 values. Furthermore, it is obvious that common values of @ and 8 which 
provide stability over all €, can be found from those values of G@ and 8B yielding 
stability for the condition € = 0.01, since it is this value which imposes the most 
severe stability requirement. Also, the plots suggest that increased gain might 
help the overall situation across the range of varying €, since the ¢= 1.0 con- 
dition most assuredly benefits from increased gain, while the effects of increased 
gain upon the € = 0.01 condition can hardly place that locus much closer to the 
real axis. Increased gain, however, may improve stability (for ¢= 0.01) by 
shifting the frequency points further out along the locus and therefore further away 
from the region of criticality, Fig. 4-18 detail. Thus, values of 475.3 only 
place thé critical points depressing ly close to the plant's Nyquist loci (all ¢), while 
for values of 8 > 0.7, a large portion of the circular segment for 0<a< 5,3 hes 
in a region of instability for much of the range of (7 0.1. The largest overall 
degree of stability across the varying range of § appears to be given by 
meee, 0-4 < 5. 3. 

It is reasonable, therefore, to begin the search for the @ and 8 which 
meet specifications with starting values of @=2, B=0.5. The specifications are 


taken to be 


t = 2 seconds = T 
r r 

o = 15 percent = O 
Ss S 

t = 4 seconds = [| 
Ss S 


Since the problem includes a search for the proper gain setting, Ko this parameter 


+s chosen to have a nominal value of unity, which represents the original conditions 
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for which this preliminary analysis has been conducted. It is noted that the plots 
suggest a final value of gain for a which is greater than unity, i.e., the Nyquist 
curves of Fig. 6-10 indicate that increased gain while not affecting the ¢ = 0.01 
and ¢ = 0.1 curves substantially, ought to improve the response for the C= 1.0 
condition considerably. Thus, the starting values for the search are taken to 
be a@=2, B =0.5, le 
Also, 
(uss 2) (t= Tce opt On) ane eeiaie (6-13) 
Cc a S s ss 
where 


lh 


r 


Zz O.=15, To-4 (6-14) 


The step sizes chosen for the intial exploratory move are as follows: 


For KO: take step 5 =) 2 (6-15a) 
For a4 : take step a Oo! (6-15b) 
For 3 : take step 6, = mos (6-15c) 


These exploratory increments have been chosen arbitrarily and, for explanatory 
purposes, they have been made distinct; obviously, they can also be set equal, 

if desired. The subscripts 1, 2,3 are assigned to represent step increments along 
the Ko QandB dimensional axes, respectively. Usually, the problem and the 
designer's experience suggest basic step sizes for each independent variable. 
During the termination procedure, the nominal step sizes of the increments are 
reduced until no further reduction is possible due either to solution indifference 
or equipment limitations. 

Due to the pressure of time the search was conducted on analog-computer 
equipment rather than digital; the most satisfactory equipment for this type of 
problem is a hybrid computer containing both analog and digital capability, such 
equipment was not available. The computer was a Mid-Century Instrumatic 
(5068), a Mid-Century MC-300 six-channel recorder was used to obtain traces of 


system input signal, response, error, etc. The computation scheme is shown in 
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Pag. 6-11. If in Fig. 6-8, one lets 


u sat] . 
is aes] (KO = 1 temporarily) (6-16) 
then 
I 
er ee + = =) 
Bu 2 tan es (6-17) 
Also, if 
=e 
vu =—ésSS (6-18) 
then 
Z 
Tl wn 
u,  st2Cw ee 
1 n 
and 
Vv, = i -2¢u $46 6 
] Ss n° oon" (onc), 
Finally, if one sets 
vo = KY] (6-21) 


then Fig. 6-11 is a represenation of the above relations. 

It should be recognized that as a result of performing the search on analog 
equipment, the step-response characteristics (t., Oo. ts) of the system were 
determined visually by reading the pen-recorder traces for v(t). Once a near- 
optimal solution was obtained, it was not possible to reduce the increment sizes 
(6) any further since any differences in the pen-recordings, for smaller 5 values, 
were not distinguishable. 

Generally, the difficulties in meeting the specifications are worse at the 
€=1.0 end of plant variability, i.e., worse for t and to therefore, the search 
is started at that end. The initial base point for the search is designated, bi, 

Fig. 6-12,the horizontal bar signifies the vector (multi-dimensional) nature of 
titsepoint, i.e., 
es VES, (6-22) 
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Fig.6-12 : "Direct Search" with ''Pattern Strategy" 


Only the first two iterations are shown 
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which is obtained from 


a a=2, B=0.5 (6-23) 


The step response of the system with a compensator which possesses such char- 


acteristics is 


(t_. Ot.) = (4. 8, 0, 9) (6-24) 
Therefore, 


Ml, 2,0.5) = 14. g-2)°+[0-1517+ [9-4] = 257.81 (6-25) 


This value for the artificial function, f, is now compared with the value obtained 
at the next observation. The next observation is performed at the point b, +8, 
where 6 is a vector with all components zero except the one representing the 


1 
independent variable, K Thus, 


b, +8) = (i222. Ss) (6-26) 
since %, | = 5) =0.2, from (6-15a). From these parameter values the response 
is found to be 

(t,o..t) = (3.5,.0762 8) (6-27) 
Therefore, 
onan e 2 
mee, 2, 0.5) = L3. 5-2)°+Llo-15]° +16. 8-4]° = 235. 05 (6-28) 


In general, when there is little or no intuitive feel for a potential solution and one 
direction offers as much hope for success as another during the initial exploratory 
moves,the response for b, -4, is also obtained to determine if it yields better 


] il 
1° In this case, since 


results than bao.) or b 


f(b, +5,) = f(b ,) (6-29) 


and, since it is recognized that a sluggish system (¢ =1.0) becomes even more 





Lod 


sluggish for reduced gain, the value of f produced by b, +8, 1s accepted as the 


smallest and b, +5 is designated the temporary base point, ty): 
indicate this 1s the first temporary head (in a vector sense) and only one variable 


The subscripts 


(Ko) has been perturbed. Thus, ae represents the temporary base point about 
which the next exploratory move for the second independent variable, 4, is made. 
To retain consistency of notation, let 


tig = by (6-30) 
To find to 1.e., the temporary head obtained from perturbing the second 
independent variable, the response for Ta +6, is first determined. Since b, 


Me@resentS an increment to Q of 0.1, the response for 


t,t = { 125. 250085) (6-31) 


must be found. The performance index for these parameter values yields 


(t oot) ivi Zae2eOe ss) (6-32) 
Therefore, 
mal 2, 2. 1, 0.5) = Eee iengis els.al¢ = 242.0 (6235) 
The response for 
ec = ( lye2s Oana) (6-34) 


is now found. The performance index becomes 


(t_.o_,t.) = (3, 0, 6. 8) (6-35) 
Therefore, 
Me 2. 1..9,.0/5), = 13-21° + (0-15 )¢ + Rememeieresens (6-36) 
Since 
i Sy) oul (chao) (6-37) 


the next temporary head is taken to be 
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To find, ty) 3: the third and last (in this example) independent variable, 8, is per- 


turbed by o.- First, the response for 


(t),+5,) = (1.2, 1.9, 0. 55) (6238) 


is found. 


The performance index yields 


(t O. t) = ( 3280 0G 5.28) (6-39) 
Therefore, 
mplee2, 1. 9, 0.55) = (3. 8-2) + Loss)“ +s. ee4)- Sean ae (6-40) 
The response for 
lope = (1. 2, 1.9, 0. 45) (6-41) 


is now found. The performance index yields 


CEE, 0 25) (6-42) 
Therefore, 
mie 2, 1.9, 0. 45) = L3,9-21°+(0-15]°4[5. 5-4] = 230. 86 (6-43) 
Since, 
£(t,,- 5) — f(t 1ot 5, ) < f(t, ,) (6-44) 


the next temporary head is taken to be 


t 2 =t 2105 See (6-45) 


All the independent variables have now been perturbed; the last temporary head, 
14 is taken to be the new, or second, base point, b,. 

With the establishment of the second base point, b,, the pattern strategy 
is brought to bear. The pattern move uses the information generated in the explora- 


tory moves which, initially, establish the pattern and later modify the pattern, as 
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the example is to demonstrate. The pattern moves are formulated upon the phi- 
losophy of ''nothing succeeds like success''. Based upon the success in the past, 
of motion to Be from by, the assumption is made that continued motion in the 
same direction remains the most desirable course of action. Thus, the new 
temporary head for the second set of exploratory moves is obtained by doubling 


the magnitude of (b5-b,} and proceeding inthe same direction. Hence, 


ts =b,+ 2(b,-b,) =b,+(b,- b,) = 2b, - by (6-46) 
Therefore, from (6-22) and (6-45) 
tog =2b,-b, = 2(1, 2,1. 9, 0, 45) - (15 2205) (6-47) 
= (2. 4,-3.0870. 9) a 220s) 
= (1-4 Ue Oe) (6-48) 
The performance index becomes 
(ti, Ocgnt a) Gag Cee) (6-49) 
cus 
Therefore, 
mo 4, 1.8, 0. 4) = 4-21" +L0-151° + lano-4l seemed (6-50) 
The response for 
ioe). = (12440. 2,1. 8, 0.4) = (io) Ween) (6-51) 


ZO, 


is now found. The performance index yields 


it, 0. taluasen0s. 6) (6-52) 
1g S S 
Therefore, , 
mile, 1.8, 0, 4) = 13-21 #lo-15) + hanea4l seen (6-53) 
feaeeround that 
t.2 0 = 6-54 
Mee ql <a (6-54) 


If one reasons as previously with regard to the gain, Ke the calculation for 


f= j t, be taken to be 
t56 by can be skipped and toy can be taken 
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—— 


tat > Soon (6-55) 
mite search for te now begins. To this end, the response for tog te is found, 
where _ o 

tat Op = Oy oe) (6-56) 
The performance index yields 
SMe ee So I 0h eb) (6-57) 
Therefore, 
mil, 1.9, 0.4) = L3.1-21°+ lo-15 f +04, 5-41° = eon ae (6-58) 
The response for 
eo =(1ei6el Oe) (6-59) 
yields 
(t O- t.) = 13. 000-.4..0) (6-60) 
Therefore, 
amen 704) = (e=2 1-4 oes ee Weal) 
Again, since 
f(t,)- 54) <f(t, + Sa)s f(t, 4) (6-62) 
the next temporary head is taken to be 
too aa 2 re 
io find, tos, the last independent variable, 8, is perturbed by First, 
the response is found for 
a 6. = 6-64 
Bor. (1.6, 1.7, 0. 45) ( ) 
The performance index yields 
(t Oust 0 = "(oro vom) (6-65) 
rs s 


Therefore, 
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lero, lq 245) eal o: g-21°+Lo-151°+6-41? = 232. 24 (6-66) 


The response is found for 


ee 


Bes = (Gul a Omen (6-67) 


The performance index yields 


(t .o4t.) = (26 024.5) (6-68) 
Therefore, 
: 2 2 2 
£(1.6.1.7,0.35) = U2, 8-2) aloe ed ar eee eee 
(6-69) 
Since 
ee ee = ; 
HES 3) f(t, 5) eaten: ,) (6-70) 
the next temporary head is taken to be 
jo] to = 
So ie lie mes erie 


Since all the independent variables have again been perturbed the last temporary 
head, to: is taken as the new, or third, base point, b,. The process is continued; 
from (6-45) and (6-71) is found 

ie 2b,-b, (6-72) 
e( 1.6; Vo 7703.35) — U2 
(3.2, 3.4, Oe) Clee. 1 oO) 
Sen Onl. 55:0..25) (6-73) 


1} 


It is noted that repeated success in the same direction, for all the independent 
variables, has caused the pattern to grow. Thus, although only two iterations 
mieemeeen Conducted, the gain value, KO is at 2.0 after starting at 1.0 and being 
constrained to basic step sizes 1 te Oran 


The performance index becomes for t 3 


(t_,0.,t.) =127 0 esec) (6-74) 
Stas 





Pie 


Therefore, 


£(2, 1.5, 0.25) = |2<2 oe lloe 75 le eee ele (6-75) 


It is noted that the set of values 
(K,4,8) = (270.71 Ole) (6-76) 


is the first triple for which the specifications are met by all three elements of 
the vector performance index (t O. t.). The results of the search to this point 
are encouraging. The artificial function f which has served the purpose of 
finding a surface of acceptable performance no longer provides a suitable com- 


parison test; the magnitudes of the individual differences 


It -T |, lo -O_], lt_- T_| (6-77) 
cst ss ss 
are monitored in an attempt to distinguish ''best'' performance. The response for 


tt 5, = (2.2, 1.5, 0.25) (6-78) 
yields a performance index 
({. 70,0 =— (le tere) (6-79) 
r' s'’ Ss 


Therefore, 
£(2.2, 1.5, 0.25) =L1. 7-2\°+ loots 14 ee gEAo 2 oneene (6-80) 


The response is still without overshoot, so the exploratory move (tf, 794) is 
again omitted and 


t..=t.+ 6 6-81 
tsi= Sot %y ( 


The choice has been made, not because 


f(t,ot 8) < f(t, 0) (6-82) 


although that is true, but because the response for yeu By is further inside 
specifications than that for ty 9: This is discussed in detail below. 


The search for ty, is begun. First, the response for 


6. = (6-83) 
t, 75, (2,2. \6n0825} 
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is found. The performance index gives 


(t.,0_,t.) = (2, 0, 4) (6-84) 
Therefore, 
Li 2 2 2 
£(2, 2. 1..6, 0: 25) = (222) -4 losis | ieee eee (6-85) 
For = = 
t.,-5, = (2, 2940. 25) ; (6-86) 
the performance index yields 
(t, oO. t.) = 2005 3, 5) (6-87) 


At this point, the investigator faces a quandary. With the boundary of 
acceptable performance pierced, and the artificial function f effectively dis- 
carded, how is the decision regarding the ''best'' temporary base point to be made? 


The three possibilities, sei). ae yytell fe oF yield performance indices 


"30° "30° 71 30 

(1.7, 0, 3.8), (2, 0, 4) and (2,0, 3.5), respectively. All are within specifications; 
which is best? While taot 5) yields a set just meeting requirements, and on that 
basis could be discarded, how can a decision be made between t29 and bor 6)? 
Obviously, what the one gains on rise time the other gains on settling time. Two 
possible techniques suggest themselves: a)a coin is tossed, i.e., one invokes the 
Principle of Insufficient Reason, b)the specifications are reduced; a new artificial 


function, f,, is formed, and another comparison measure of which set (K, a, 8B) 


| eg 
is ''best'' is thereby provided. 


To this end, a new set of requirements 
=1.5, O =15, T =3.0 (6-88) 
if S S 


are specified; it should be noted that rise time and settling time have been 


retained at the same ratio as previously. Under these circumstances, 


I 


f f,(2.2,1.5, 0,25) = C1. 7-1. 5)7+ (0-15) G. 8-3]° a e225neg 


1't3)) 
£y(t,;7 65) - ¢ 


(6-89) 


(252-7. oeaees) = (as | Calne ieee 226.25 
: (6-90) 


(2.2.1. 4, 0,25) = |2-1.5)-+ onistemaesaens 225.50 
- (6-91) 


It 


fi) (t,)- >) 





174 


Therefore, since 


f(t, ,-6,) < f(t.) = f,(t,,7 5 5) (6-92) 


the next temporary head is taken to be 


i325 ieee So's) 
Hro= find, ae the variable, 8, is augmented by 5.. First, the response 
for ence ‘ 
+ = 
t.5" 6. (2.22, 1. 4.03510) (6-94) 
is found. The performance index yields 
(t_. oe t) = (2, 0, 3.5) (6-95) 
Therefore, 
22a 150 0) = La. eae: Locale = 225. 50 
(6-96) 
The response for 
ta5-54 = (2.2, 1.4, 0. 2) (6-97) 
is found. The performance index yields 
(er ec ls eis Secs (6-98) 
r s s 
iieretore, 
2 c “eS 
1325 Lach. 2) = (1. 9-1. 5)]°# 0215) 413-4230 = 22a 
(6-99) 
Since 
t..- 6.) <£(t,,.) = £(t,,+ 5 6-100 
fon 3) i) 
the next temporary head is taken.to be 
t. = t ~& ——-} 6-101 
om se 3 D4 


Since all the independent variables have again been perturbed, the last 
temporary head, ty: is taken as the new, or fourth, base point, ive From (6-71) 
and (6-101) is calculated 
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es 


ago eb,-b, (6-102) 


(2. .2,.1,4,,0,2)— (1) oye Oe) 
(4.4.2. 830.08) - (loge? -0235) 


= hc oy ol 0.405) (6-103) 
The response for, or is determined and the characteristics yield 
(f. .O-t je = (ie 4 Oe) ; ‘6-104) 
r oS as 
Waarefore, 
f (2. gS. ) 20. 05) = wees Aye lioonG les Le. eee (6-105) 
The new function, fi: must now be discarded since the set of values for Ko Ozer 
(6-103), provides a response which again is within the specifications on te Ot. 


Thus, a second new comparison function is formed for the changed requirements 


Ties, OF = 00, eleaeaee (6-106) 
be S S 
Therefore, 
ie. 8, 1.1, 0.05) = GE Lane hake, 721° = MO. 72, (6-107) 
Miteeresponse for 
ao = -1 
Payot Ope woe Oo Me Nae) (6-108) 


is now found. The performance index yields 
((.2,90: 5°t = (le eae) (6-109) 
rs Ss 
Therefore, 
: 2 2. Jo 
f(s, 0) tee. 05) = Licd-1) ee lOel0e Ze s-2 1 loos 
(6-110) 
Since, the response is again without overshoot, the exploratory move (t35-5)) is 


again omitted. Also, 


— 


7 < t = el 
fa(tygt 64) < fatyo) (6-111) 


so 
Toe te (6=112) 
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aihe search for, aoe is begun. First, the response is found for 


t4,+5, = (3¢ lo2G0s) (6-113) 


The performance index is 


(t,o.,t.) = (2,0, 2.9) (Gaines) 
Therefore, 
£,(3, 1.2, 0.05) = 2-1 )-nkosloy al cea) mtomeen (Gamuey 
For 
tay7 o> F (3201200105) (6-116) 
the performance index gives 
(t,0.,t) = la? 0c. 4) (6-117) 
Therefore, 
= Sree =a 2 [ Z iva 
ayy )) = £,(3.1, 0,05) = 11. 7-1]°+lo-10]°+(L2. 4-2]*° = 100. 65 (6-118) 
Thus, since 
t. Pose o ne fe ts = 
Ey) ~ POM ee ere / oi), 
the next temporary head is taken to be 
tao = ty) (6-120) 
Motind, on the variable, 3, is increased by oe First, the response is found 
for - Z: 
ap = = 
ta ° 3 (320.312 aoe) (6-121) 
a hesperfoOrmance index yields 
tt 220 t Jia Gls 10 are) (6-122) 
1 Sies 


Therefore, 
Sewn al) = (Ly 811°" 10-101-+ [302221 2 oes ose 


t 12 3 (3. GO Ie 3 ) (6- ) 
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cannot be used due to practical limitations on the analog computer, i.e., due 
primarily to the layout used to solve the problem. An even more important 
justification for not reducing 8 to zero, however, is that it is not possible to 
realize such a compensation network, within the given class to which the com- 
pensator design is constrained, Furthermore, since the artificial function has 


been reduced twice already, the design is well within original specifications, for 


oe !.0. Thus, since ; 


f{t,5) < £,(t,,+ 53) (6-125) 


ae f,(t,,) = f(t, ,) = f(t 407 >) <f,(t,5+ 5) (6-126) 


FOO. 25> = 02-508 


the next temporary head is taken to be 


t 43 = tyo—rbe (6-127) 


Since all the independent variables have again been perturbed, the last temporary 


head, t 43° 
is calculated 


t.. = 2b_-b 6-128 
a6 2b, by ( ) 


2135. Vee O05 e222 lee ee) 
(6, 2.2, O- l) = (22, 1 ae) 
= (3. 8, 0. 8, -0. 1) ——#(3. 8, 0. 8, 0. 05) (6-129) 


is taken as the new, or fifth, basepoint, by. From (6-97) and (6-108) 


ch 
| 


That is, under the practical constraint that 8 not be less than 0.05, 3 is set 
equal to 0.05 whenever the process indicates B02 05. 


The response for 


€, = (3. 8, Os S90n 0} (6-130) 
50 
is determined. The performance index yields 
(t_.,0../ t= (1.47 08Ze Ss (6-131) 
rs s 
iiterefore, 


£,(3.8, 0. 8, 0. 05) = i. 4-1 oman [2 5-2 = 100,41 (ollees 








Lis 


The response is now found for 


teat o) = (4.0, 0. 8, 0. 05) 
The performance index yields 
(t On t) = (122,022.20) 


Therefore, : 


f,(4,0.8, 0.05) =[1, 2-1]°+lo-10}"+ lz-2}* = 100. 04 


meain, the exploratory move, tags Oa is omitted, and since 
+O daca ote 
falta gt Op) = f(t 9) 
Sie CS Onan 
The search for, t-., is begun. First, the response is found for 


BY 


t.1+5, 4 
te tos (A On 7 2) 10>) 


The performance index gives 
(t,o, t.)) =) 0a) 
r s Ss 


Therefore, 


f,(4, 0. 9, 0. 05) Pi tal oles ee 1 = eh. oI 


Por, 


tey7 ° 5 = (4, 02 70205} 


the performance index yields 


(tro. tira Ol oe 
Eos os 


Therefore, 


2 


Thus, since 


(40) 7 0805) = [1 3-11°+ fo-10 412. 5-216 = 1o0ns4 


(6-133) 


(6-134) 


(6-135) 


(6-136) 


(6-137) 


(6-138) 


(6-139) 


(6-140) 


(6-141) 


(6-142) 
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f(t, + 55) 2 fF tee < Gea) (6-143) 


the next temporary head is taken to be 


tao = tot 5, (6-144) 
To find, te the variable, 8, is increased by 6,. First, the response is found for 
teat, = (4,0.9, 0.1) (6-145) 


The performance index yields 


(t.,o.,t.) s(t oe O cae) (6-146) 


Therefore, 


£,(4, 0.9, 0.1) = C1, 2-1/+lo-10+l2,2221- = logs vemmniCenes 


Again, the response for t 6, cannot be found so the tot 6, response 


oe 
is matched against the tz5 response and it is found that 


f(t.) es f(t,,+5,) (6-148) 


4 


Therefore, the next temporary head is taken to be 


ta =t, 5b, (6-149) 


Since all the independent variables have again been perturbed, the last temporary 


head, is taken as the new, or sixth, base point, by. From (6-127) and (6-149) 


bc 3° 
is found 


te = 2b, -b, (6-150) 


2(4,,0,9, OF 05) 2 als.el L, 0, 05) 
(8, 1.8, 0.51) te 1, 0.05) 
= (5, 0. 707 05) (6-151) 


The response is determined for 


teo = (5, 0.770205) (6-152) 


The performance index yields 
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(to ie) = (lalvO52e1) 


Ss 
Therefore, 


(27047 0.105) = [1.9-1 14 orn ies De. Tail toon 


The response for 


trot Sy = (5.2409) O05) 


is now found. 


The performance index yields 
(t_,o.,t.) = (1.1,0,2. 2) 


Therefore, 


fy (242, 0-7, 0.05) c= leslie 121 ea oentodiao tua = 100. 05 


The response for 


ec =e, = (4. 8, 0.770. 05) 


60 
if found. The performance index yields 
(t ot) =alal, lO. geeees) xy 
Therefore, 
7 Z 2 f 
E456, 0. 7,005) 7 1) 1=1) “ lostolsal 2222 lecees 
Since 


(Go oo 


fa(ho) ~ falte gt 21) = folteg- oy 


the next temporary head (basepoint) is chosen to be 


St a0 


The search is started for te 93 the response is found for 


t, +A 


617 &2 = (5, 0. 6, 0, @5) 


The performance index measures 


(t,o; t..) ="( 1305, Ge, 
r’ s’ Ss 


(6-153) 


(6-154) 


6-155) 


(6-156) 


(6-157) 


(6-158) 


(6-159) 


(6-160) 


(6-161) 


(6-162) 


(6-163) 


(6-164) 
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Therefore, 
ie) 030) 05)e— bi. 05-1]°+Lo-10]°+L2-2]% = 100. 0025 (6-165) 
Bor, _ _ 
Pee SiGe 0s, 0.05, (6-166) 


the performance index yields 


(to, i =1( le OZaoy (6-167) 
Therefore, 
fie 9510.16, 0.05)! = Giant Wee ne ee [2-215 =s1oumol (6-168) 
Since, 2° an - 
f(t jt 5.) < fi(tei- 5 5) s f(t, ,) (6-169) 


the temporary head is taken to be 


— — — 


Uae Wonle ie Boo 0) 
To find, tea the third independent variable, 8, is increased by, 5, and the 
response is found for . 
{.3t ose 2 
too 3 (5,205.0. 1) (6-171) 
The performance index yields 
(t_,0_,t_) = (1, 0, 2) (6-172) 
ros Ss 
Therefore, 
£,(5, 0.8, 0.1) = (1-1]*+Lo-10]7+ L2-2]* = 100 G17) 


cannot be found so the i +6 response is 


Again, the response for Ue © or Ss 


3 
matched against the teo response and it is found that 


Ate) 6G. (6-174) 


Therefore, the next temporary head is taken to be 





to 3 > teat Sa-—aebe (6-175) 


Since all the independent variables have been perturbed, the last temporary head, 
tea is taken as the new, or seventh, basepoint, b Thus, from (6-149) and 


(6-175) 


7° 


tog = 2b2-b, (6-176) 
2(5, 0. 8, 0.1) - (4, 0.9, 0. 05) 
(10, 16; 0.2)" = (450, 90205) 


= (6,.057,,0. 15) (62177) 


If 


The response for 


70 


can be determined, etc. The process can be continued, thus providing much 


= (6,0. 7, 0-5) (6-178) 


additional information regarding the correlation between (K a, B) space and the 
space of (t oe t_). A solution providing satisfactory response across the range 
of varying € is (K a,3) =(5,1.0,0.05). The values of the performance index 

for various ¢ are displayed in Table 6-1. It is clear from Table 6-1 that, 

as Seo —-1.0,. and 2=0.05 forma satisfactory compensator which meets the 


Original specifications. 





TABLE 621 

‘ pe fe pepe 

< r S 
Mmm kc mimes] 
0n9. | 1095 On oS 
Oss) KOS ITO | eee 
nO a | Oe | CO 
oe eel 
(2 On oe 
MEM | Ube 
a 
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T2207 eed (6-14) 


for all€ inthe range, 0.01 <¢< 1.0. Indeed, this compensator meets the more 


stringent requirements of 
ie =1, O =10, T =2 (6-106) 


across the entire range specified for varying ¢. 

A review of the example indicates how quickly the search converges upon 
a sOlution; a digital computation program can, of course, provide the answer 
obtained with considerably more speed and accuracy. The addition of the delay 
time (pure lag), which is the major change to the example, introduced by the 
human operator,in now way modifies the search procedure used. The delay 
interval has its effect upon the system response characteristics and may reduce 
the prospect of success in finding a ''best-overall'' compensator of a given class 
but the direct search method with pattern strategy retains its general applicability. 

It should be noted again that, in the example just offered, the human- 
operator transfer function was set equal to unity and the variable-parameter- 
mapping plots used to estimate the approximate values that @ and 8 must have 
if compensation is to achieve the purpose prescribed. Thus, the three-variable 
mappings were allowed to degenerate to a two-variable representation and the 
values of 4 and 2 inthis instance, had a direct effect upon the location of the 
critical points governing stability. Normally, when the {) mappings are used to 
represent the human-operator adaptive range the influence of 4 and 8 values 
upon system stability is not so easily discernible. In such cases, the selection of 
values for @ and 3 influences the shape of the -L curve, which represents 
the system dynamics; the manner in which this curve passes throught the 0 
regions is, of course, what ultimately determines overall (with a man in the loop) 
system stability and performance. Thus, in an analysis with a human operator 
im the system, the effects of the Ko a and 3 values are only indirectly observable 
before they are selected. Upon selection, however, if the compensated -Lo is 
plotted the effects are clearly visible. This situation is no different than that 
routinely faced by the control engineer when dealing with less animate systems. 


Experience is quickly gained concerning the results of introducing particular com- 
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pensator configurations (picking values for a and 8); to this end, with a definite 
goal for the shape of the modified or curve the Bode’ diagram (more quickly 


plotted than a Nyquist curve) can be used to advantage. 


Ge 7 Summary 


A discussion Of scalar performance indices and some of their limitations 
is followed by a laconic survey of methods for solving the "'best''’ performance 
problem, especially for varying-parameter plants. The notion of a ''vector'' 
performance index is introduced and a technique for establishing one response as 
better'' than another is presented. This technique rests intimately upon the 
formulation of an "artificial" function which when reduced to zero identically 
satisfies the upper bounds on system requirements. These requirements are set 
forth in the ''vector'' index. The design for a ''best-overall'' compensator to 
meet specifications upon a varying-parameter plant is obtained by using a "direct 
search" method (especially suited to the electronic computer). This search 
method incorporates a "'pattern strategy'' to aid in finding the next set of trial 
values for the compensator; the strategy affords rapid convergence to a solution 
as is demonstrated by an example. 

Finally, the example itself, displays the versatility of the variable-para- 
meter mappings presented in Chapters Two-Five. By setting the human-operator 
transfer function equal to unity, it is shown that the mappings have use for 
“inanimate'' system design. 

Empirical justification of the conclusions from this theory (with a human 
operator in the loop) must wait upon additional research and the publication of 
experimental results issuing from deterministic inputs driving controlled elements 
other than pure gains. It is recognized that the bulk of the evidence justifying the 
form of human-operator transfer function used in this study has been accumulated 
for systems driven by stochastic input signals. Still, the initial range selection 
for the parameters KO si T, > T, and Ty included the effects of calculations pub- 
lished for deterministic inputs and modifications, if needed, are easily incor- 
porated. With a technique provided for comparing one multi-dimensional index 
against another, the methods of this chapter offer the prospect for more realistic 


performance measures than the standard ''figure of merit'' of scalar variety affords. 
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APPENDIXA 


fe Human Operator Zero Locus 


It is desired to establish (4-13) as the equation of a circle. The trans- 


Mermiation (A-1) is to be considered 


_atjb- (A-1) 
aa dae fi 


w=utjv= 


Of interest are the conditions x -0, y “wT, ; thus 


atjb _ (i +x)-jy 
(l+x)tjy (1 +x)-jy 


We = 


ee a(l+x)+by+jl b(1+x)-ay ] (A -2) 
Cac 
play, 


_ 


The real and imaginary parts of each side of (A-2) are set equal 


— ate (A -3) 
(letaxa) ey 


oe eiisey (A -4) 
Say 
emaatlarly, 


atjb ; u-jV 


lt+xtyy ee noe 


 autbvt+j(bu-av) 


ie tas 


Ptxrjy 
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therefore 


+ 
eee ee (A-5) 
Zee 
WW) ty 
Dien, 
OY =isenes (A-6) 
i 


For x=0, (A-5) gives the locus in the space of the mapping and (A-6) permits 


recovery of Tr fOr any given 7 


Ce 
u tv -au-bv = 0 (A-7) 


If the square in(A-7) is completed, one obtains the familiar equation of a circle 


Z Ze eee ts a 
(. 3] i(. a = Veena (A-8) 


ze Z eg 
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APPENDUSs 


mele the Parabolic Effect 


A typical plot of a face from (T, E » Ty) space, T,_ constant, appears in 


I L 
Fig. B-1l. Without loss of generality, a set of xy axes is shown so oriented as to 
simplify the following exposition. Of interest is the complete mapping corres- 
ponding to face ODCEO in Fig. B-2. This mapping can be found by allowing a 
typical element of this face such as ain Fig. B-2 to scan from line OD to line 
EC. The effect, as has been seen in Chapter Four, is a plot similar to that 
shown in Fig. B-!, except for translation of axes. When qis halfway down line 
DC, 68 is halfway down OE, Fig. B-2, and similarly in the mapping space, 

‘Fig. B-1, i.e., linearity holds through the mapping as is proved later in this 
Appendix. Thus, a’ is also halfway down the line DG’ of Fig. B-leecnianem 


halfway along O’E’. Each of these points is given arbitrary coordinate 


designation in Fig. B-1l andit is noted 


for this Bhoice of axes. The mapping of any element a8 is obtained by joining the 
points aa’ in Fig. B-!. As shown in Fig. B-3 the family of these lines extends 
outside the end point mapping of Fig. B-1|. It is shown that the envelope of this 
family of lines 4‘a' is a parabolic arc as suggested by the many dotted lines of 


fimeetamily drawn in Fig. B-3. For the geometry of Fig. B-1! 


ae) 
x 
/ 

oo Ye ee 
y 3. 5 3 
2 (oe eee 

0, > A ee 


=) = + 
a, (y, ¥4) V4 


5 / 
where the subscripts indicate the x and y components of the moving points q and 
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y 
(Ks Yo) 
c 
E} (x,y, } 
a 
B 
D ae 
fo) (X65. ) 


Pie, B-1? Typical Map itonm (it, Ty) Space 


L I 





aoeds (Ny 





I 


AL 6 


g 


L) ut eoeq e Butuuedg : 27-g “31 y 


EES 
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ENVELOPE YIELDS ARC 
OF PARABOLA 


Fig. B-3: The Parabolic Effect 
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g’ in the space of the mapping. From these designations may be written the slope 


of the moving line: 


@ -o hy, -Aly,-y,)-y 
d 3 2 
SS Se Gee eee ee (B-1) 
dx a’ -g’ Os 
x x Z 


Either by integration and evaluation of the constant of integration (when x =0, 
y ay) or through use of the slope intercept form y =mxtb the equation for the 


family of lines may be written 


y pe a y, 
(ie ae 4 
Ys e—  aeng (B-2) 
g ae 
From the theory of envelopes it is known that poe 
1 Sonie) =: (B-3) 


represents a family of curves, then differentiating with respect toi, holding x 


and y constant and cquating to zero, yields 
dF 
See ce) ee B-4 


The envelope is obtained by climinating ). between (B-3) and (B-4). 











ED na y: 
oe 8 4 
es 0) > ee ee x-Ay, = 0 
“—) ae) 
and 
wal 
= _ 7g) 
x 4A 
y.X =X 
i (B-5) 
(1) 


Ford, L.R., “Differential Equations”, McGraw Hill, New York, 1955. 
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If (B-4) is substituted into (B-2) one obtains, 


bee 4 = ere 
x5 (B-6) 


Equation (B-6) upon simplification becomes 











a — 
i fx (B-7) 
me 
z 
my a) 
a 
y = ae (B-8) 
i 
2 
A change in notation yields 
as 
Pee ooo 4 
= — and b= 
Z 
y -ax=b./x (B-9) 
Equation (B-9) squared becomes 
(S Zee 
y -2axyta x = bx (B-10) 


Equation (B-10) is a quadratic form whose parabolic nature is not 


obvious. Standard formulas for transformation of coordinates to principal axes 


(1) 


form are available. Briefly, a quadratic form 
Z 2 
Axatexy + Cy Dx thle (B-11) 


may have its axes rotated by an angle q so as to eliminate the cross-product 


(1) Thomas, G.B., Jr., "Calculus and Analytic Geometry", Addison- 
Wesley, 1956. 





eke, 


term. A general rotation applied to the axes of Eq. (B-11) results in a new 


quadratic form 


2 
A‘’x’ BY cy 4 Cty ee =a) (B-12) 


with coefficient relations given by 


z &, 
A’ =Acos go. + B €0S q sin aq 17 -simen 


2 Z 
= B(cos a-sin qa)+2(C-A)sin acos q 


B’ = 

/ ae . 2 

C =Asin o-Bsinacos q+tC cos q (B-13) 
D’=DcosatE sina 

BE’ ——Oesia 7 t+ Dees 

no = 


If (B-11) has coefficient B# 0 then an angle q can always be chosen such that 


- 


B’ =0; this angle may be found from 


oe 26 = (B-14) 


With the above information and some trigonometric manipulation one 


fmacdcethat since 








a 
_ti-a 
GOt ca. a 
then 
2 
cos 2q=-> 
a. th 


SO, 


sina = = (B-15) 
pA seal 
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and 
1 
€0s-00= =a (B-16) 
Sa@tl 
Use of (B-15) and (B-16) in Eqs. (B-13) produces 
Z Z 
2 a ab b 
(a“+1l)y/°+ ey’ - —— =x’ = 0 (B-17) 
aZ+l a“tl 
which has the form of a parabola. However, let 
2 Z 
WS and See 
{a +1) (a +1) 
then 
2 
y-hx’ +ky’ = 0 
and completing the square 
2 
» k& b 
4-3 —.-__—=. = 
(y 5) a ies 0 
2 Z 
; ik pas 
il s aees = aaa 
(y 5) h(x Tel 
Let 
Wp / res 
y ay. 2 
and ”) 
x” = x’ aL ah 
then 
me 
y = hx” (B-18) 


Equation (B-18) is the simple form of the equation of a parabola. 
It should be noted that the area added to the end point mapping to form 


the complete hull is very small. In the light of this smallness, its effect upon 





iOS 


the gross hull over the range of T and its location relative to the plots of 


is 
most physical systems, there is sufficient justification to call the mapping 


"almost complete" in the spirit of Section 4.3. 


; A p ; ; p 
Bb. 2 n Exploration of Tots Space 


The prospect of defining exactly the mapping of the rectangular 
parallelipiped of Fig. B-2 onto the 2) plane must be considered. Previous results 
were for the mapping of the boundaries of this volume only. As above for face 
ODCEO, the envelope of the region in) space corresponding to the mapping of 
every single point within the volume shown in Fig. B-2 must be obtained if both 
necessary and sufficient conditions for stability are to be provided for the dual 
Nyquist plot. 

The situation w=1 givenin Fig. B-4 (same as Fig. 4-18) is first 
considered. It should be noted that face 1243 of Fig. B-2 maps into the triangle 
plus short line segment 1'2°4°3" and face 5687 maps into the triangle and short 
line segment 5'6°8°7". Indeed, as face 1243 is moved along lines 1-5, 2-6, 3-7 
and 4-8 it eventually coincides with face 5687. Thus, every point of the volume 
al 


am ( E , Ty) space is scanned as face 1243 is moved along for varying T 


4 
goes from zero to 5.3. The effect in the complex plane plot is to 


i 1 
1.@€., as Ty 
have the triangular shaped mapping move from the position 1'2'4"3° to that of 

5'6'8'7'. This motion occurs under the constraint that point 2‘ moves along arc 


2‘6’, i.e., the motion of point 2 along the line 2-6 in (T,_, Tp Ty) space corres- 


Wy 
ponds to motion of 2’ along the circular arc 2’6° in the space of the mapping. 
Similar situations pertain to points 1’,3° and 4’ in the complex plane as they 
move toward coincidence with 5’,7° and 8’, respectively. 

It is desired to show that lines 2’-4’ and 4’-3’, as they move toward 


lines 6-8’ and 7'-8’, never contain points which lie outside the bounds of the 


hull plotted in Fig. B-4. Later, the results of the parabolic effects displayed 





Ty Burkreg joq033q : p-q “81g 


WVL3S0 33S 


WVL30 
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(.9).2 
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in the first part of this Appendix are considered. 
For the situation of Fig. B-4, i.e., w=1, the line E’C ‘(4’-2") has end 


points 1+j25 and -16.5+j25.7, therefore its equation is 


V=2o ry x-1 
25.7-25 -16.5-1 


(. 7) 








Me VG 
or 
= — -l)+ = 
yi 17.5 (x-1) +25 (B-19) 
i.e., the line has slope 
— iat 
aes 


The end points of the radius arm of arc F’E’ are 1 +j25 and the center 


l 
of arc ee’. Tee <.; 5 t jle. 5. The equation of the radius arm is therefore 





Via eo 
WoO aoe ole , 
g 
=12.5 
y-25 = a (x-1) 
“2 
or 
y = 25(x=))zs5 (B-20) 
i.e., the line has slope 
= 25 
m, 
Since 
-.7 ey a 
mom. = 2 Ses = -] 
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then line E’C’ is perpendicular to the radius arm OnE -E" or E’C’ (2 278 
is tangent to arc F’E’ (4'8’). The above is true for w=1. It must be shown 


Mae, in general, for any w. 


é ° 
For E‘’C*: The end points are 1+j25w and eile, Gio j¢5. 7w therefore 


line E’C’ has equation 


y-25W 2 x- 1 
25. 7w-25w -17. But 
es cet! een 
1 -17. 5ur iow 


For radius arm Onp-E’ in Fig. B-4 the end points are 1+j25w andthe 





ec 51 : : 
center One at ae the equation of Dag is then 
y-25w _x-!l 
12.5W-25w Le 
Z 
= =) 22 ow _ 
m, = aa = 25w 
a2 
and 
Fede A. ot Verl fae ot ae 
> > 17.50 (a rn 


indicating, again, that oT. 2 is perpendicular to line E‘C’, and therefore line 
EG is tangent to arc F'E’. 
A similar analysis answers the question, "Is ib) (Ge tangent to arc AD?" 
For D‘C’: The end points are 1+j.7” and 1-17. 5 w+ 25.7jW theretore 


the line D’C’ has equation 


2 
25. TP-. TW Pe ler 5 we 
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y-. /7w a x =] 
ae Pieest 
= 25 w 
and | - eninge) (B-21) 
-17.5w 


mor the radius arm of arc A’ D 00-0) 9. oO en eee Fig. B-4 the 








AD 0.7 AD 
end points are 1+j.7w and the center como j ea giving an equation for the 
radius arm of 
-./w —= x-] 
Tw oil 
2 aor 1 W 2 = l 
-. 35 
m, = = a = .7W (B-22) 
Equations (B-21) and (B-22) give 
m)m, = nt (a7 teal 
= lf 


and line- D'C’ is tangent to arc PSS IDY 

Thus, it has been proved that initially line 2°- 4° is tangent to the cir- 
cular arc F’E’ (48°) at point 4‘, Ifin (T,. Ty Ty) space, face 1243 moves 
some distance along the T, axis toward face 5687 then in the space of the map- 
ping, point a experiences some corresponding angular excursion along the arc 
2°6’ toward 6’. Point 4’ experiences the same angular excursion along the 
arc 4°8’ toward 8’, etc. 

Actually, it is sufficient to show that for a specified increment along Th 
equal angular increments occur on all the arc segments, i.e., it shall be proved 
Mite =x T,, O< x< l CT 


L L £ 
ment increments having equal angular measure. The proof is facilitated by the 


=max T,) corresponds, as x varies, to arc seg- 


passing of all circles through the origin. Thus, it is only necessary to show that 


a particular increase along qT) produces equivalent phase change on all the cir- 
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cles. 


For arc B’C’ (2'6’), since 


Ai esalle 5 w+ 25. 7jw 
i 1+T, jw 


then, for T. =0 


t Zz 
Oo SAC Se) commen 
and 
” = faveal eS 
}-17.5w 


fore |.=x T 


L ve 5 
q = Lal? 5H )+ 25. 7ju 
1+xT jw 
and 
@ = tan”) “2.18 tan”! (xT, w) 
1-17.5w 


Therefore, the extent of arc 2'6° is given by 


noe ee Tw (B-23) 


For arc_F’‘E’ (4°8°), since 


Ante eee ie 
\ a 1+T, ju 


then, for Ty =) 


a ycore 


and -] 
tan. 25 


R 
I 


For TL =xT, . 
eB 2 Seo 
L+xT)je 


and 


fs 


tan) 25 meee T, w 


Therefore, the extent of arc 4'8° is given by 
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CNC eg ea tan’ xT, (B-24) 


and similarly, along arc G’O’(3'7’) andarc A’D‘(1‘5‘), 

Since (B-23) equals (B-24) the statement is proved that corresponding 
angular motion occurs on all arc segments. Indeed, the fact follows from the 
properties of the bilinear transformation. (1) 

It can be mentioned at this point, also, that if Ty is set equal to xT, 
where T= max Ty and 0<x<1, similar to the above exposition, then the effect 
again is to move correspondingly along the mapped line. For example, Fig. B-5 
shows a plot, w=1, of the situation where 0< Ts 12.5, rather than 0< Ts 25%; 
as originally. The equations of lines D’-c’', O’-E’, A‘-B’ and G-F’ are 
the same as those in Fig. B-4, but the line segments cover only half of their 


previous extent. This property follows directly from the mapping 


- (l+jwT (1 +jwTQ) 


0 (amo 


For W=W), Tyy= constant, T) = constant 


= (1+j,T,) (C,+jC,) 
= (C,- C59 T t+ i(C,t C Ty) 


= [-C wt 5C, wy |T y+ Grace 


_ B-25 
Heal ) 


where K.=some complex constant. Equation (B-25) demonstrates that, irrespec- 
i 
tive of the value of Ty the mapping lies along the same straight line. 
To be considered now is the case where movement has occurred corres- 


/ / A “l / é 
ponding to 30° of arc along 2 6, shown as point 2, and 30° of arc along 4 8 , 


TY Churchill, R.V. , "Complex Variables and Applications", 2nd Ed., McGraw- 


Hill, New York, 1960. 
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OST, <12.5, OST) <0.7 AND 7, =0 
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3(0) 
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7(G) 


Fig. B-5: Effect of Reduced Range on T, 
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shown as 4, Fig. B-4. The line joining these new points corresponds in 

(T,. TY TH) space to an element of face 2486, Fig. B-2, which is parallel to line 
2-4, or equally well, line 6-8, In the complex plane the angle ek oe equals 15° 
(an inscribed angle equals one-half its intercepted arc). Also, the angle formed 
by the extended line 2‘-4' tangent to arc 4°8° and the line 4°-4” equals 15° (the 
angle between a tangent and a chord equals one-half its intercepted arc). These 
two angles are therefore "vertical angles'' and points 2’, 4°, 4” can be connected 
by a straight line. Thus, in the space of the mapping, as 2 -4° moves toward 
6‘-8°, the line connecting the moving points 2” -4” always passes through the 
original position of point 4°; thus, the line NEVER passes outside the "end 
point'' hull. Of course, this state of affairs pertains only until point 2° reaches 
4’, but after that it is obvious that the line 2-4" is inside the hull as in the case 
when 2 -4” has finally become 6°-8°. Similar conditions prevail in the pro- 
gression of line 3°-4° to coincidence with 7 =e Sie. ethene connecting any 
points a always passes through the original point 3°. Thus, except for the 
parabolic envelope effect mentioned at the beginning of this Appendix the end 


point mapping of the volume in (T al; » TH) space would be a complete one. It is 


granted that the area contributed by Be parabolic fillet adds little to the complete 
hull that was not already in the end point hull. This can be seen by envisioning 
the triangular shaped region le 43 sweeping down to 5'6'8'7', as Ty ranges 
from zero to 5.3, only keeping in mind during the sweep the presence of the 
parabolic fillet. If at any time the stability of a physical system is dependent 
upon points in the affected region they then can be taken into account. It is seen 
in Chapter Five that in those systems considered in this study the additional area 
may be disregarded. Due to the presence of this small effect the mapping is 


classified ''almost-complete. "' 
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APPEND Ee 


Gare A Iwo Parameter System 


The use of the three-parameter Variable-Parameter Mapping of Chapter 
Four for a two varying-parameter problem is considered in this Appendix. An 
interesting stability demonstration occurs in the following system; also, the 
example serves a useful purpose in the performance considerations of Chapter Six. 
The situation arises inthe quest for a best-overall compensator which is des- 
cribed in detail in Chapter 5ix; only that portion of the design task of tutorial 
pertinence is investigated. 

The system and compensator are as shown in Fig. C-l. The system is 


analyzed for values of ¢=0.1, w = 1.57. As in Chapter Four stability obtains if 
hsb IE Gl 10) (C-1) 


has no right half plane zeros, where 


ne 
12) 
n 


is eee 7 
F(s) = |e +i) | slot 2p ©) Wo) 


If (C-1)_is multiplied by 


So. 2 
=2 ae so), 


the fixed and varying portions for the plots are obtained. That is, 


2 
U 
so +l n _ C-4) 
sat] : s(s + 2y x) 
and 
vn - stl (C-5) 
VV 6©satl 


where V_ is again taken to be unity as in Section 4.2, i,e., @*= 8 = 0. 
n 


Also, PC 


L a ee (C-6) 


n s(st2ow_) 
n 





ras 






2 
Wry 


s(s+2Cw,) 






Fig. C-1: Second-Order System with Compensator 
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This results in the standard plotting structure 


= -L (cag 


Sac Lhe Degenerate Mapping 


The constraints on permissible Q@ and B aretakentobe 0<a< 5,3 
: poaiageet 

and 0<8 < 25. It is apparent that the =< mapping is obtained from the image 
of face 3487 shown in Fig. 4-16. This mapping is actually shown in Fig. 4-20, 
for w=1.0, i.e., it is represented by the bounds formed by line 3 ‘4°, arc 4‘=8%, 
line 8°7° and arc 7'-3'. The degenerate mapping just outlined is shown in Fig. 
C-2 and again in Fig. C-3, along with the plot of -Le given by the negative of 
(C-6). 

The primary purpose of this example is to demonstrate that sufficient 
caution must be exercised in interpreting the plots to preclude succumbing to the 
temptation of making stability decisions based upon limited frequency information. 
One must guard against using the thumb rule of Section 4. 2 superficially, by con- 
sidering only one or two frequency points. The remainder of the example clarifies 
this notion. Upon observation of Fig. C-3 one might infer that some set of a,B, 
such as a4=2, 8 =25, would result in a stable system. Thatis, for 8B =25, as a 


v 
varies, the image point of sa runs along arc segment 4°=8 , and for @=2 ifthe 


image point is "'to the left'' of -L_ (tracing -L_ as frequency increases), sucha 
n 
condition seems to imply a stable situation. The point corresponding to a=2, 


8 =25 and, of course, t=1.0 is shown as M(jl) in Fig. C-3 where 


7 742(25)+ 1 C-8 
Mo =e (yal (ae) 
That is, 
V 
M = ~F (a =2, 8 =25) 


and it is observed that M(jl) does fall to the left of -L. implying stability. 
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An actual stability check for these specified circumstances can be obtained 


from the Routh test. From (C-1) and the characteristics given one obtains 


(2. 46)(s2 + 1) 


PLASLESiecar (C-10) 
From (C-10) there follows 
3 Z 
s+. 354s +.20956s+.0985 =0 
Ss 1 . 20956 


S - 354 0965 


1 .0742-.0985 ~ (-) 
. 354 


>) 


S 7H) 


a 


The two sign changes in the array above signify two right-half-plane roots 
in the characteristic equation and, consequently, instability. The apparent in- 
consistency between the graphical indications and the analytical results can be 
cleared away by a more detailed plotting of the M(jw) locus at more frequencies, 
particularly those less than unity. 

The pertinent points of the locus for M(jw) are shown in Fig. C-4 along 
with -L mj. Itis noted that the ‘difference vector ‘drawn 110m -L, to M 
undergoes a clockwise rotation between M(0. 3) and M(0.7) thus indicating in- 
stability and consistency with the Routh test. The difference vector actually makes 
two complete clockwise rotations corresponding to the existence of the two right- 
half-plane roots in the characteristic equation. To observe the complete revolu- 
tions the negative frequency branches of -L_ and M must be considered; also, 
the effects of s=0 are handled by the standard technique of setting 
s = lim Re! : 

Reo 





e 
Oo 








eL0 


(ep) 


24 


22 
V 

A 7: 

y Ga 


20 


M (0.3) 
M(0.6) 


- ~ o M(0.7) 
5 w=0.5 = 
F y M(1.0) 


& M(2.0) 
\ 


elie 8 \ 
pa R(w) 
mw 2 4 é 8 io «2 


Two-Parameter Variation for 4=2, 8=25; also shows -L. 





21) 


- iT n —_ : : 
where 0 varies from x to 7) and the origin of the s-plane is circumvented by 


a small semi-circle of radius Bs 

Difficulties with ''encirclement" interpretation can normally be resolved 
by sketching the entire locus of -Lo for positive and negative frequencies 
(similarly for the varying portion); the excursions of the difference vector between 


these two loci may then be more readily observed. 


1) Del Toro, V. and S.R. Parker, Principles of Control Systems Engineering, 


foes30, McGraw-Hill, N. Y., 1960. 
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